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ABSTRACT 


This  thesis  develops  a  low-resolution  stochastic  simulation  model  to  assess  the 
impact  of  the  intelligence,  surveillance  and  reconnaissance  components  of  C^ISR,  and 
strike  capabilities  on  the  mission  success  of  a  United  States  carrier  battle  group  (CVBG). 
The  simulation  uses  a  stochastic  approach  to  model  a  two-day  conflict  between  a  CVBG 
and  a  land-based  enemy  which  incorporates  the  randomness  and  uncertainty  inherent  in 
warfare.  The  simulation  is  implemented  as  a  C-t-f-  computer  program  to  develop  a  tool  to 
analytically  exercise  a  prospective  new  system  in  order  to  predict  its  possible  effect  on 
combat  operations.  Experiments  were  run  which  simulated  a  two-day  battle  in  which  the 
United  States  CVBG  sensor  availability,  sensor  accuracy,  and  weapons  availability  were 
varied  to  study  their  affect  on  the  outcome  of  the  battle.  Statistical  analysis  techniques  are 
used  to  quantitatively  measure  the  results  of  the  battle  as  the  sensor  and  weapon 
parameters  change. 
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The  views  expressed  m  this  thesis  are  those  of  the  author  and  do  not  reflect  the 
official  policy  or  position  of  the  Department  of  Defense  or  the  United  States  Government 
The  reader  is  cautioned  that  computer  programs  developed  in  this  research  may 
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the  time  available,  to  ensure  that  the  programs  are  free  of  computational  and  logic  errors, 
they  cannot  be  considered  validated.  Any  application  of  these  programs  without 
additional  verification  is  at  the  risk  of  the  user. 
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EXECUTIVE  SUMMARY 


A.  PURPOSE 

The  purpose  of  this  thesis  is  to  develop  a  low-resolution  stochastic  simulation 
approach  to  model  a  naval  theater  conflict  against  a  land-based  enemy  which  incorporates 
the  aspects  of  the  randomness  and  uncertainty  inherent  in  warfare.  The  simulation  is 
designed  to  assess  the  impact  of  intelligence,  reconnaissance,  and  surveillance  information, 
derived  from  internal  and  external  sources,  and  strike  capabilities  (weapons  types, 
numbers,  and  characteristics),  specifically  and  illustratively  on  the  mission  success  of  a 
United  States  carrier  battle  group  (CVBG)  that  is  tasked  to  provide  support  for  a 
multinational  force  operating  in  a  hostile  region. 

B.  BACKGROUND 

Command,  control,  communications,  computer  systems,  intelligence,  surveillance, 
and  reconnaissance  (C'^ISR)  and  superior  strike  capability  are  two  keys  to  the  future 
success  of  United  States  carrier  battle  group  operations.  A  CVBG  is  frequently  tasked  to 
exert  the  policy  of  the  United  States  in  remote  areas  of  the  world  against  powerful  military 
nations.  The  leaders  that  make  the  decisions  on  the  type  of  C*ISR  and  weapons  assets  to 
develop  and  procure  need  tools  to  analytically  exercise  a  prospective  new  system  in  order 
to  predict  its  possible  effect  on  combat  operations. 

The  proposed  modeling  approach  examines  a  realistic  naval  engagement  in  which  a 
group  of  multinational  force  ships  is  initially  in  the  Persian  Gulf  and  wishes  to  transit  out  of  the 


XI 


area  through  the  Straits  of  Hormuz.  The  transit  is  opposed  by  Iranian  land-based  surface-to- 
surfece  missile  launchers  and  air  forces.  The  United  States  government  dispatches  a  CVBG 
to  protect  the  multinational  force  as  it  transits  the  straits  by  destroying  Iranian  mobile 
surface-to-surface  missile  launchers  and  defeating  Iranian  fighter,  attack  aircraft,  and 
reconnaissance  aircraft 

C.  RESULTS 

A  stochastic  simulation  model  to  assess  the  impact  of  information  and  strike 
capabilities  for  this  situation  has  been  programmed  in  the  C++  computer  language,  hiitial 
force  structures,  sensor,  and  weapons  capabilities  were  established  using  general 
unclassified  military  publications.  Computer  experiments  were  then  performed  which 
simulated  a  two-day  conflict  between  the  United  States  naval  carrier  battle  group  and  Iran. 
The  initial  conditions  of  the  computer  runs  were  varied  to  affect  the  sensor  and  weapons 
availability,  and  sensor  accuracy  for  the  United  States  and  multinational  forces.  Statistical 
analysis  techniques  were  used  to  examine  the  results  of  the  computer  simulation  runs  to 
determine  how  the  varied  sensor  and  weapon  capabilities  affected  the  outcome  of  the 
battle. 

The  analysis  of  tire  simulation  results  show  that  the  use  of  an  additional  aircraft 
carrier  combined  with  the  availability  of  satellite  intelligence  enables  the  U.  S.  carrier 
battle  group  to  complete  its  mission.  On  the  basis  of  this  initial  computer  experiment,  the 
proposed  simulation  modeling  approach  shows  promise  of  having  informative  predictive 
capabilities,  but  more  detail  in  the  models  is  required. 

xii 


L  INTRODUCTION 


A.  BACKGROUND 

Conmiand,  control,  communications,  computer  systems,  intelligence,  surveillance,  and  reconnaissance 
(C^'ISR)  is  presaitly  receiving  a  significant  amount  of  attention  in  the  Department  of  Defense.  The  United  States 
(U.  S.)  armed  services  are  concerned  with  acquiring  nearly  perfect  knowledge  of  enemy  units,  movements,  and 
plans  while  denying  the  enemy  the  same  informatioiL  C'iSR  includes  the  assignment  of  the  sensors  used  to 
obtain  data,  communication  links  used  to  transmit  and  receive  data,  the  capabilities  and  priorities  of  the 
machines  and  men  that  process  raw  data,  and  the  commander’s  perception  of  the  military  situation  on  the  basis 
of  the  analyzed  end  product.  The  leaders  that  make  decisions  on  the  ^pe  of  ClSR  assets  to  develop  and  procure 
need  tools  to  analytically  exercise  a  prospective  new  system  in  order  to  predict  its  possible  effect  on  combat 
operations.  Modeling  and  simulation  tools  are,  and  have  potential  to  be  used  extensively  the  Department  of 
Defense  to  influence  these  decisions.  This  thesis  used  a  ^pedfic  naval  conflict  to  construct  a  stochastic 
simulation  model  that  incorporates  some  of  the  randomness  and  uncertainty  inherent  in  warfare.  The  model  is 
programmed  as  a  computer  simulation  in  the  C++  computer  language  Simulation  trials  are  run  to  erplore  the 
effects  of  changing  conditions  that  relate  to  C’iSR  which  are  programmed  into  tire  model  on  the  outcome  of  tire 
naval  conflict 

The  scenario  that  the  model  was  designed  for  is  a  realistic  naval  scenario  that  may  be  of  interest  to  U. 
S.  military  planners.  The  scenario  envisions  a  groip  of  multinational  force  units  that  is  in  the  Persian  Gulf  and 
that  wishes  to  transit  out  of  the  area  through  the  Straits  of  Hormuz;  its  transit  is  opposed  by  Iranian  land-based 
surfece-to-surfece  missile  launchers  and  air  forces.  The  multinational  unit’s  exit  is  assisted  by  a  U.  S.  carrier 
battle  group  (CVBG).  The  interplay  of  C‘*ISR  and  strike  c^abilities  is  judged  by  use  of  the  ejploratory 
simulation  model  that  is  implemented  in  this  thesis. 

B.  THESIS  OBJECTIVE 

The  objective  of  this  thesis  is  to  develop  a  stochastic  modeling  ^proach  that  can  assess  the  impact  of 
intelligence,  reconnaissance,  and  surveillance  information,  derived  fi-om  internal  and  external  sources,  and  strike 
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capabilities  (wejqwns  types,  numbers,  and  characteristics),  on  the  mission  success  of  a  U.  S.  CVBG  in  the 
scenario  outlined  above.  Simulation  is  used  to  study  the  outcomes  of  the  stochastic  model. 

C  APPROACH 

This  thesis  develops  a  low-resolution  closed-loop  stochastic  simulation  model  of  a  specific  naval 
scenario.  The  element  of  chance  is  incoiporated  into  the  model  by  assigning  q)ecific  probabilities  for  the 
occurrence  of  certain  significant  events  (e.g.  detection  of  a  target,  weapon  hits  on  a  target,  operation  of  radar 
systems,  ^c.).  When  the  model  is  ta^ed  to  determine  the  outcome  of  a  significant  event,  the  outeome  is 
simulated  using  appropriate  probabilities.  The  model  design  attenq)ts  to  represent  the  impact  of  various  types  of 
informational  sensors  (own  unit  and  external),  the  battlefield  perception  drawn  fi-om  the  combined  sensor 
information,  and  the  subsequent  number  of  assets  utilized  and  losses  e^qteiienced  while  attempting  to  ccmqtlete 
an  assigned  missioa  This  modeling  ^proach  has  been  realized  as  a  time-stq)  simulation  programmed  in  the 
C++  computer  language  to  produce  quantitative  results.  The  units  that  comprise  the  Iranian  and  United  States 
forces,  and  the  sensor  and  weapons  capabilities  of  each  force  were  d^ned  for  this  model  to  allow  for  quantitative 
analysis  of  simulation  runs.  Analysis  of  the  simulation  results  using  a  specific  initial  force  stracture  and  unit 
characteristics  has  shown  that  this  model  does  not  contain  sufficient  detail  to  properly  assess  the  affect  of  C*ISR. 
This  model  presents  an  approach  that  has  promise  of  appraising  the  affea  of  C4ISR  on  U.  S.  naval  operations  if 
more  detail  is  programmed  into  the  coitq)uter  simulation. 
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n.  SITUATION 


A.  GENERAL  DESCRIPTION 

The  situation  considered  is  that  a  multinational  force,  consisting  of  a  few  U.  S.  combatants  escorting 
merchant  ships,  is  located  in  the  Persian  Gulf  Iran  decides  to  prevent  die  force  from  moving  through  the  Straits 
of  Hormuz,  and  threatens  aircraft  and  surfece-to-surfece  missile  (SSM)  attacks  on  the  vessels.  The  prospect  of  a 
complex  war  in  this  region  precludes  the  U.  S.  from  obtaining  any  friendly  air  or  littoral  support  The  U.  S. 
sends  a  CVBG  to  protect  the  force  whenever  it  chooses  to  transit  through  the  straits  to  the  Arabian  Sea.  The 
situation  is  shown  in  Figure  1. 


Figure  1.  Conflict  Region  and  Approximate  Locations  of  Opposing  Forces 


The  goal  of  the  multinational  force  is  to  transit  the  straits  into  the  Arabian  Sea,  but  to  do  so  while 
experiencing  minimal  damage.  The  CVBG  is  tasted  to  protect  the  multinational  force  destroying  Iranian 
SSM  weapon  sites  around  the  Straits  while  minimizing  losses  of  aircraft  and  Hamags  to  the  hattip  group.  The 
goal  of  the  Iranian  forces  is  to  prevent  both  the  multinational  force  from  complying  its  transit,  and  to  Tnavimiyp: 
damage  to  the  U.  S.  naval  forces  assigned  to  protect  the  mnltinafinnal  force. 
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B.  ORDER  OF  BATTLE 


The  multinational  force  consists  of  1  DDG  (Burke  class  with  vertical  launch  ^^stem  (VLS)X  1  DDG 
(Spruanoe  class  with  VLS),  and  2  Merchants.  The  U.S.  CVBG  consists  of  1  CVN  (Nimitz)  and  an  augmented 
aircraft  contingent  which  includes  46  F/A-18  strike  aircraft,  30  F-14  fighter,  4  E-2C  early  warning  aircraft,  and 
4  reconnaissance  aircraft,  2  CGs  (Ticonderoga  class  wATS),  3  DDGs  (Spruance/Bufice  w/  VLS),  and  1  AOR 
The  mission  of  the  U.S.  carrier  is  to  deter  Iranian  air  strikes  by  destroying  their  attack  aircraft,  attaching  their  air 
bases,  destroying  their  SSM  launchers,  and  eliminating  their  primary  sensors  (reconnaissance  aircraft  and  long 
range  radar  sites);  the  warships  are  tasked  to  provide  defense  for  the  aircraft  carrier  and  strike  land-based  targets 
using  Tomahawk  cruise  missiles.  The  United  States  force  can  also  have  two  satellites  available  for  overhead 
visual  and  electronic  intelhgence  gathering. 

The  Iranian  force  consists  of  land-based  attack,  fighter,  and  reconnaissance  aircraft  squadrons,  SSM, 
and  SAM  sites.  The  Iranian’s  primary  sensors  are  land-based  long  range  radar,  and  reconnaissance  aircraft. 

C  STUDIED  SITUATION 

The  assumptions  of  the  model  for  the  scenario  are  as  follows.  In  the  corrqruter  implementation^  an 
external  file  is  used  to  define  the  initial  force  strocture  for  both  sides.  This  external  database  defines  each  unit  by 
type,  weapons  systems  (type,  capacity,  range),  and  movemait  patterns.  The  characteristics,  tactics,  and  locations 
of  all  units  and  we^ns  incorporated  in  this  mockl  are  derived  fiom  unclassified  general  pubhcations  and 
information.  The  initial  database  is  designed  to  force  a  two-day  battle  between  United  States  and  Iranian  forces. 
As  the  battle  commences,  the  multinational  force  is  located  in  the  southern  Persian  Gulf  transiting  towards  the 
Straits  of  Hormuz,  and  the  CYBG  is  transiting  towards  a  rendezvous  point  southeast  of  the  straits.  The  tracks 
that  all  ships  will  follow  during  the  two-day  battle  are  pie-determined,  and  ontyf  their  rate  of  movement  varies 
around  an  average  transit  speed.  The  multinational  force  travels  through  the  Straits  of  Hormuz,  joins  the 
CVBG,  and  transits  towards  the  Arabian  Sea.  Both  the  CVBG  and  the  multinational  force  follow  tracks  that 
maintain  them  at  the  maximum  range  fi’om  Iran  while  remaining  in  suffidentty  deep  water.  The  Iranian  forces 
commence  hostilities  as  soon  as  they  detect  any  naval  units,  and  the  U.  S.  forces  launch  pre-emptive  strikes  when 
they  detect  the  Iranians.  The  battle  concludes  after  48  hoins,  and  on  the  average  the  CVBG  and  the 
multinational  force  (assuming  the  multinational  force  survived)  have  reached  the  entrance  to  the  Arabian  Sea. 
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When  the  United  States  has  satellite  assets  available,  the  satellites  examine  regions  of  the  Iranian 
coxmttyside  where  it  is  suspected  that  operating  air  bases  and  mobile  SSM  launchers  are  located.  These  pre¬ 
determined  flight  patterns  result  in  possible  detections.  Ihe  satellite  Cc^rability  is  represented  by  two  sensors. 
One  of  the  sensors  performs  a  three  hour  scan  of  a  r^on  of  Iran  where  threat  contacts  are  located  and  then  takes 
nine  hours  to  orbit  the  earth  before  returning  to  conduct  another  scan  of  the  same  area  The  other  sensor 
performs  a  2.5  hour  scan  over  a  different  area  where  threat  contacts  are  located  and  then  takes  six  hours  to  orbit 
the  earth  before  returning  to  conduct  another  scan  of  the  same  area.  All  other  sensors  are  located  on  moving  or 
stationary  objects  and  the  area  where  one  of  these  sensors  has  the  possibility  of  detecting  an  enemy  unit  depends 
on  both  the  location  and  range  of  a  given  sensor,  which  is  defined  for  each  sensor  in  the  initial  data  base. 

D.  TACTICAL  ASSUMPTIONS 

1.  This  scenario  excludes  any  combat  between  naval  units,  other  than  between  carrier  aircraft  and 
land-based  strike  and  recomiaissance  aircraft. 

2.  Intelligence  sources  are  restricted  to  the  sensors  in  the  simulation,  and  Human  Tnt<»lligftnf«» 
(HUMINT)  is  not  built  into  the  model. 

3  The  naval  umts  attached  to  the  CVBG  or  multinational  task  force  move  in  Jtormatioti  There  is  no 
independent  steaming  of  ships,  and  they  proceed  along  pre-deteimined  courses.  The  movement  track  of  the 
CVBG  ensures  that  all  units  remain  out  of  range  of  Iranian  SSM  launchers.  The  CVBG  units  therefore  are  onty 
susceptible  to  attacks  by  aircraft. 

4.  The  naval  units  are  considered  to  be  independent  entities  for  detection  purposes.  Therefore,  the 
detection  of  one  unit  does  not  reveal  the  location  of  another  unit.  In  practice,  identi^g  the  location  of  one  ship 
of  a  naval  task  force  implies  that  other  ships  are  nearby,  but  this  is  ignored  in  this  model. 

5.  The  dissemination  of  information  between  fiiendly  units  is  assumed  to  be  perfect  and  effectively 
instantaneous.  This  irt^jlies  that  information  available  to  one  unit  is  available  with  only  a  small  ne^gible  time 
delay  to  all  units,  so  a  common  intelligence  picture  can  be  formulated. 

6.  The  Iranian  forces  utilize  continuous  radar  coverage  to  detect  an  incoming  U.S.  strike.  The  United 
States  Navy  operates  using  the  Aegis  radar  system,  which  is  resident  on  Ticonderoga  class  cruisers,  as  the 
primary  air  defense  asset  protecting  a  carrier  battle  groi^p.  Use  of  the  Aegis  radar  ^stem  exposes  the  location  of 
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the  ship  it  is  located  on  to  the  enemy  electronic  intelligence  sensors.  The  U.  S.  believes  that  the  benefits  of  using 
this  ^stem  outweigh  ejqwsing  the  location  of  naval  assets  to  enemy  electronic  intelligence  sensors  (EUNI). 

7.  SAM  sites  are  modeled  as  if  they  are  collocated  to  an  air  base,  SSM,  or  command  and  control 
center  and  are  represented  as  a  close-in  weapons  q^stem  (defensive  asset).  ThQ^  are  assumed  to  only  engage  an 
enemy  unit  that  is  attacking  the  facility  that  is  collocated  with  the  SAM  site. 

8.  Once  a  missile  launch  occurs,  it  may  only  be  defeated  by  deception,  jamming,  or  a  close-in  weapons 
^em.  Airaaft  are  not  given  the  capability  to  shoot  missiles  out  of  the  sky. 

9.  Iranian  sui^ly  dqpots  that  reload  mobile  SSM  launchers  with  missiles  are  disbibuted  throughout  the 
model  of  the  Iranian  country.  They  cannot  be  detected  sensor  assets  (i.e.  are  effectively  undetectable),  and 
therefore  are  never  targeted  or  destroyed.  This  assumption  is  not  necessarily  realistic,  and  can  be  changed  in 
future  revisions. 

10.  The  unit  that  has  the  least  distance  between  itself  and  an  enenty  target,  and  that  has  a  weapon 
onboard  that  can  attack  the  enemy,  will  fire  at  the  target 

Note:  Various  of  the  restrictive  assun^rtions  made  above  can  be  altered  in  later  revisions  of  the  software. 
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in.  MODEL 


A.  GEOGRAPHIC  REPRESENTATION 

The  region  is  represented  by  a  700  X  400  mile  m^  consisting  of 20  X  20  mite  grids  (700  grids,  20 
rows  by  35  columns).  The  20  mile  grids  maintain  a  ship  within  the  same  grid  for  at  least  30  minutes  rturing  a 
transit  (the  typical  transit  speed  of  a  merchant  is  15  knots,  and  a  naval  warship  is  20  knots).  Each  grid  is 
numbered  sequentially  as  shown  in  Figure  2.  A  grid  contains  information  identifying  the  grid  boundaries,  the 
type  of  grid  (land  or  water),  the  general  weather  condition  over  the  grid,  separate  hsts  of  United  States  and 
Iranian  umts  that  are  located  on  the  grid,  and  separate  lists  of  the  Iranian  units  that  the  U.  S.  forces  have  detected 
within  the  grid,  and  U.  S.  units  that  the  Ranian  forces  have  detected. 
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Figure  2.  Geographic  grid  and  numbering  system 

The  database  to  build  this  region  is  located  in  ^pendix  A  and  is  diqrlayed  gr^hically  in  Figure  3.  F^rh  list  of 
units  that  are  located  on  the  grid  identify^  individual  units  Ity  a  unit  number  and  Cartesian  location,  i.e.  they 
spedfyf  the  current  ground  truth.  Both  the  U.  S.  and  haman  perception  of  the  location  of  enemy  units,  which  is 
realistically  imperfect,  is  maintained  using  two  additional  lists  that  reflect  the  perceived  (not  actual)  inhahitantg 
of  this  grid  itystem.  Movement  of  naval  units  is  restricted  to  sea  grids,  and  ground  units  to  ground  grids,  while 
aircraft  can  move  freefy  among  both  types  of  grids. 
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Figure  3.  Geographic  Data  Base  as  Stored  in  the  Computer  Simulation 
B.  COMBATUNTTS 
1.  Overview 

Status  of  the  militaiy  units  is  maintained  in  a  40-element  air^.  Each  individual  platform  (U.  S.  ship, 
U,  S.  reconnaissance  aircraft,  U.  S.  satellite,  Iranian  air  base,  Iranian  SSM  launcher,  Iranian  reconnaissance 
aircraft,  or  Iraman  radar  station)  is  identified  by  a  unique  unit  number  and  the  information  describing  the 
operational  parameters  of  each  unit;  individual  parameters  are  modified  whenever  an  event  occurs  that  changes 
a  parameter  during  the  course  of  the  simulation.  The  information  that  is  stored  for  each  unit  is  as  follows: 

a.  Unit  Characteristics 

1)  unit  number 

2)  type  of  unit 

3)  current  electronic  emission  status 

4)  number  of  hits  to  kill 

5)  number  of  hits  against  unit 

6)  unit  is  targeted  (incoming  weapon) 

b.  Location/Movement 

1)  grid  location 

2)  X  “  coordinate 

3)  y  -  coordinate 

4)  motion  status 
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5)  next  waypoint 

6)  time  at  current  waypoint 

7)  course 

8)  speed 

9)  nominal  operating  speed 

10)  waypoints  (six  different  waypoints  are  stored  for  each  unit) 

1 1)  loiter  time  at  a  wj^^int  (six  different  loiter  times  are  stored  for  each  unit) 

Note:  a  w^point  is  a  position  defined  by  X  and  Y  coordinates,  that  a  unit  is  to  pass  through 
at  some  time  during  the  conduct  of  the  simulatioa  Th^  are  sequentially  numbered  from  one  to  six,  with  the 
niunber  indicating  the  order  in  which  a  unit  will  pass  through  each  waypoint 

c.  Weapons  Capability 

1)  number  of  missiles  or  attack  aircraft  carried 

2)  effective  strike  range  of  missile  or  attack  aircraft 

3)  number  of  fighter  aircraft 

4)  effective  range  of  fighter  aiiciaft 

Note:  effective  range  of  an  aircraft  is  the  nuiriber  of  miles  the  aircraft  can  travel  from  its 
home  base  (assumes  normal  operation  speed)  before  it  must  commence  its  return  journey  (i.e.  sufficient  fuel). 

d.  Sensor  Capability 

1)  sur&ce  sensor  maximum  detection 

2)  air  sensor  maximum  detection  range 

3)  ESM  sensor  maximum  detection  range 

The  database,  which  specifies  the  parameters  of  the  characteristics  of  each  unit,  constructed  for  use  in 
the  performance  of  computer  trials  and  analysis  in  this  thesis  is  included  in  Appendix  B.  The  base  operational 
parameters  of  each  unit  (speed,  we^ns  load  out,  weapons  range,  number  of  aircraft,  sensor  ranges)  are 
estimated  fiom  unclassified  publications  which  describe  the  military  assets  of  the  United  States  and  Iranian 
governments;  (Sharpe,  1996)  and  (Gunston,  1980).  Specific  irtitial  parameters  are  changed  during  t^sring  and 
exercise  of  the  simulation  for  analysis.  The  total  numbers  of  weapons  fired  during  the  conflict,  and  the 
subsequent  losses,  are  stored  for  analysis  purposes  as  separate  variables. 

2.  Weapons 

The  U.  S.  weapon  systems  that  are  modeled  are  ship-launched  land-attack  cruise  missiles 
(Tomahawk  conventional  missiles),  carrier-based  fighter  aircraft  (F-14)  and  carrier-based  attack  aircraft 
(F/A-18).  The  Iranian  weapon  systems  modeled  are  land-based  fighter  (F-14)  and  attack  aircraft,  and 
surface-to-surface  missiles  (Exocet)  launched  from  mobile  land-based  platforms.  Each  unit  has  a  limited 
supply  of  weapons  that  is  carried  onboard.  Aircraft  carriers  and  air  bases  are  given  an  initial  number  of 
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jBghter  and  attack  aircraft  as  its  weapons,  and  U.  S.  naval  warships  cany  Tomahawk  missiles.  The 
Iranian  SSM  launchers  are  the  only  units  that  can  replenish  their  weapons  supply  during  the  simulation. 
Each  SSM  launcher  begins  the  simulation  with  an  initial  loadout  of  A/ missiles,  and,  when  they  are 
expended,  the  SSM  launcher  transits  to  a  supply  depot  and  receives  new  missiles.  Other  types  of  weapons 
that  are  not  represented  in  an  effort  to  simplify  the  model  could  be  easily  added  later  as  an  improvement 
of  the  model. 

Each  unit  maintains  counters  that  keep  track  of  the  number  of  weapons  it  has  available  at  a  given 
time,  and  variables  which  state  the  effective  range  of  each  weapon  carried.  After  a  unit  assigns  a  weapon 
to  engage  a  target,  the  number  of  weapons  available  for  the  unit  is  correspondingly  decreased.  The  model 
also  maintains  counters  that  track  the  total  numbers  of  each  weapon  type  fired  by  the  U,  S.  and  Iranian 
forces,  and  the  number  of  weapons  destroyed  by  the  defensive  mechanisms  of  their  intended  targets. 

3.  Sensors 

The  different  types  of  intelligence  gathering  ^sterns  that  are  rqjresented  in  this  model  include  surface 
search  and  air  search  radar  (ship  and  land  based),  aiibome  radar,  visual  detection  by  pilots,  electronic  support 
measures  (ESM),  satelhte  imagery,  and  satellite  electronic  intelligence  (ELINT).  The  sensors  are  modeled 
using  a  cookie-cutter  approach.  In  the  model  each  unit  is  given  three  different  types  of  sensors  with  their 
associated  effective  detection  ranges,  and  probabilities  of  detecting  a  contact  given  that  it  is  within  the  sensor’s 
effective  detection  range.  The  three  sensors  are  designated  surfece  search,  air  search,  and  ESM  search.  For 
instance,  the  surface  search  sensor  represents  a  surfece  search  radar  for  the  carrier,  while  it  represents  a  pilot’s 
eyes  for  a  reconnaissance  aircraft.  All  units  with  the  exception  of  the  mobile  SSM  launchers  and  reconnaissance 
aircraft  are  assumed  to  be  constantly  emitting  (operating  radar).  On  every  hour  and  half-hour  the  SSM 
launcher  decides  whether  it  operates  its  radar  for  the  next  half-hour.  A  BemoulH  trial  (probability  equal  to  0.5) 
determines  whether  the  launcher  operates  its  radar  continuously  during  the  next  half-hour.  If  the  SSM  launhcer 
is  moving,  this  trial  is  not  conducted  and  its  radar  remains  off.  Using  this  approach  it  is  possible  that  the  SSM 
launcher  continuously  operates  its  radar  for  long  periods  of  time.  A  more  complex  model  should  be  used  to 
represent  the  SSM  launcher  radar  operations  which  takes  into  account  the  location  of  enenry  surface  units  and 
reconnaissance  aircraft  which  can  detect  the  SSM  launcher.  When  the  SSM  launchers  are  transiting,  or  are  at 
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the  missile  reload  site,  thqr  do  not  operate  their  radar.  The  U.  S.  and  Iranian  reconnaissance  aircraft  operate 
their  radar  whenever  they  are  in  fhght,  and  secure  their  radar  the  moment  they  land. 

C  DECISION  MODULES 
1.  Movement 

a.  United  States  Naval  Units 

Each  ship  has  a  set  of  ^  (IT  =  6  for  the  simulation  runs  conducted  in  this  thesis)  waypoints 
designated  prior  to  running  the  simulation.  The  ships  will  transit  in  straight  line  paths  between  waypoints  at 
speed  V^p,  where  is  distributed  as  a  normal  random  variable  with  mean  pvsh^  (nominal  operating  speed) 
and  variance  cTvship^  conditioned  to  be  positive.  To  simulate  this  random  variable,  the  simulation  re-computes  the 
value  of  V^ip  imtil  a  value  greater  than  zero  is  obtained.  The  variance  accounts  for  deviations  in  course  and 
speed  attributed  to  variations  in  wind  and  current  conditions,  persormel  errors,  and  equipment  operational 
variability.  The  equations  for  confuting  the  course  and  sampling  for  the  speed  of  the  ship  are  shown  in 
equations  (1)  and  (2).  All  courses  in  this  simulation  are  corrqiuted  and  stored  using  the  radian  as  the  unit  of 
measurement. 

Ax  =  0,  Ay  >  0 

Course  =  -^2  Ax  =  0,  Ay  <  0 

arctan(Ax/Ay)  Ax  0  (1) 

(Pvsh^^vsh^)  (2) 

where  A^fp^avshy)  represents  the  distribution  of  a  normal  random  variable  with  mean  p.,ship 
and  variance  vsh^that  is  conditioned  to  be  positive 

Once  a  ship  reaches  a  waypoint,  it  may  loiter  at  the  wtypoint  throughout  a  pre-determined  time  interval  or  it 
may  immediately  commence  its  transit  to  the  next  wtypoint.  If  the  ship  is  programmed  to  remain  at  a  waypoint, 
its  motion  will  be  stopped  until  it  reaches  the  time  for  it  to  commence  the  transit  to  the  next  waypoint  The 
computed  track  of  each  ship  will  not  aUow  it  to  traverse  land  grids  and  maintains  each  ship  within  flie 
geographic  boundaries.  The  algorithm  that  determines  ship  motion  is  described  in  Appendix  C. 

Improvements  to  this  model  should  allow  both  the  carrier  battle  group  and  multinational  force  ships 
to  react  to  incoming  strike  aircraft,  missiles,  and  reconnaissance  aircraft  by  changing  their  course  and  speed  to 
avoid  them.  The  enertty  perception  module  that  is  described  later  does  not  compute  an  enemy  target  track.  In 
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the  current  model  the  motion  of  the  multinational  force  commences  at  the  beginning  of  the  simulation;  an 
improved  model  could  assess  the  number  of  Iranian  assets  destroyed  by  the  CVBG  strikes,  and  start  the 
multinational  force  movement  through  the  straits  after  a  set  amount  of  Iranian  weapons  capability  has  been 
destroyed. 

b.  Iranian  Mobile  Surface-to^urf^e  MissUe  Units  (SSM) 

Mobile  SSM  sites  will  move  from  their  present  location  to  a  randomfy  determined  position 
(waypoint),  or  to  the  missile  reload  site.  On  the  hour,  and  half-hour,  the  SSM  unit  will  use  a  Bernoulli  trial 
(probability  defined  in  the  initial  data  base)  to  decide  whether  it  should  maintain  its  jx)sition,  or  conGpute  a  new 
positioa  New  positions  are  computed  using  equations  (3)  and  (4).  The  new  X  position  (XssiUnew )  is  computed 
by  drawing  an  independent  random  variable  fi'om  a  normal  distribution  with  mean  zero  and  standard  deviation 
o^x_diange  adding  this  number  to  the  current  X  position.  The  new  Y  position  (Y ssiUnew  )  is  computed  similarly 
using  a  normal  distribution  with  mean  zero  and  standard  deviation  ay^djangc-  The  ax  diange  and  ay  diange  values 
represent  a  typical  movement  distance  for  the  SSM  launcher  (one  standard  deviation)  and  is  defined  by  the  user 
of  the  model,  the  value  of  10  miles  was  arbitrarily  selected  for  quantitative  analysis.  The  course  to  travel  to  this 
new  position  is  computed  using  equation  (1). 

XSSninew  -^(0,  Ox^char^)  Xcutrent_position  (^X_diai:^ge  ”10)  (3) 

YSSnijjeT^  ^Y_diange )  YcyijQj^_poj^^  (cry_chaqge  10)  (4) 

where  A^(0,  ax  d^nge)  represents  the  distribution  of  a  normal  random  variable  with  mean  0  and  ^ 
variance  ax  diange,  and  A^(0,  ay  diange )  represents  the  distribution  of  a  normal  random  variable  with  mean  0  and 
variance  ay^diange- 

The  new  position  may  be  toward  the  coastline  (closer  to  enemy  naval  ships)  or  inland,  and  the 
direction  is  purely  by  chance.  On  the  average  the  mobile  SSM  launchers  will  remain  in  the  same  general  area, 
which  maintains  the  launchers  in  range  of  ships  transiting  near  the  Iranian  coast,  but  confounds  the  United 
States  intelligence  picture  because  of  random  changes  in  locatioa 

Once  a  new  position  is  computed,  and  before  the  SSM  starts  moving,  two  checks  are  performed  to 
ensure  the  new  position  is  acceptable.  The  first  check  ensures  the  new  position  is  within  the  boimdaries  of  the 
geographic  map.  If  the  position  is  off  the  map,  equations  (5)  and  (6)  compute  a  new  position  that  is  within  the 
geographic  boundaries.  This  check  first  determines  if  either  the  new  X  or  Y  coordinate  (Xssninew  Yssrrw)  is 
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negative;  if  a  coordinate  is  negative  it  is  set  to  zero  (Xleftjboundary  =  Ylowerjboundaiy  =  0),  which  is  an 
acceptable  X  or  Y  coordinate.  It  then  determines  if  the  new  X  coordinate  is  greater  than  the  maxitmim  X 
coordinate  ofthe  map  (XjdghtJboimdary= 700);  if  it  is  greater,  it  is  set  equal  to  this  maximum  value.  Ifthenew 
Y  coordinate  is  greater  than  the  maximum  Y  coordinate  (Y iqrperjboundary  =  400),  it  is  similarly  set  equal  to  its 
maximum  possible  value.  This  boundary  check  is  appbed  for  any  unit  that  has  just  computed  a  new  position  to 
transit  toward. 

^Qeft_boundaiy  if  ^^SSTUnew  ^  -^Qeft_boundaiy 

XSSmjiew  ^Qi^_boundaiy  ifXsSUlnew  ^  ^Ci^tJxjundaiy  (5) 

XsSninew  if  ^eftjboundary  '^XsSUlnew  Xri^Jxjundaiy 


Yssnw  = 


^  lowcr_boundaiy 
Y^)per_bouncfary 
YSSUlnew 


ifYSSninew  ^  Yiowtr_boundaiy 
ifYSSUloew  ^  Y^,pCT_bcnmdaty 
if  Yioiwer_.boundary  ^  YsSUlnew  Yuppeii^boundaiy 


(6) 


The  second  check  that  is  performed  ensures  the  new  position  is  located  on  a  nuiribered  grid  that  is  designated 
land.  The  first  step  of  this  check  determines  the  numbered  grid  (Gridnew)  where  the  unit  is  to  be  located  once  its 
complete  its  journey  to  its  new  position  CXisiiinew  Yssmiw).  If  Gridnew  is  a  land  grid,  this  new  position  is 
acceptable  and  the  check  is  complete.  If  the  Gridnew  is  a  sea  grid,  the  position  is  changed  using  equations  (7)  and 
(8)  which  compute  a  new  position  that  has  a  better  chance  of  being  within  a  land  grid.  In  equation  (7)  and  (8) 
Gridpr^ent  is  the  grid  where  the  unit  is  presently  located.  Xcomsction  and  Ycoirectionare  numbers  used  to  shift  Xssrrw 
and  Yssnincw  so  that  the  new  SSM  position  is  located  on  a  numbered  grid  which  is  a  land  grid.  The  values  of 
XcorT«:tion  Bud  Yconectioo  ^  sct  cqual  to  ouc  quaTtcr  of  the  length  and  width  of  a  grid  For  this  simulation,  the  grids 
are  20  miles  by  20  miles  and  the  values  of  the  Xcoirectbn  and  Ycoi«ction  are  both  five  miles.  This  check  is  repeatedly 
performed  until  the  Gridnew  is  a  land  grid 

XsSmnew  “  XSSrrinew  ^correction  if  Gridnew(COl)  -  Gridp(csent(COl)  0  (7) 

^^^SStUnew  "  Xcorrection  if  Gridnew  (OOl)  *■  Gtldpi^esentC^l)  ^  ^ 


YSSItlnew  YSSninew  Ycom0ction 
Yssmnew  “  Ycoirection 


if  Gridnew  (row)  Gridpresait(col)  <  0  (8) 
if  Gridnew  (row)  -  Grid^^(col)  >  0 


Once  an  acceptable  new  position  (waypoint)  is  identified,  the  SSM  launcher  commences  its  journey  after  a  one 
minute  time  delay  at  speed  where  is  distributed  as  an  independent  normal  random  variable  with  mean 
jrvssm  (nominal  operating  speed)  and  variance  cfvssm^  conditioned  to  be  positive.  V^sm  is  computed  using  equation 
(9).  The  simulation  re-computes  the  value  of  Vssm  until  a  value  greater  than  zero  is  obtained  to  produce  a 
replication  of  the  conditioned  normal  random  variable. 
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^ssm  (^vssnb^wam) 


(9) 


If  the  SSM  is  moving  on  the  hour  or  half-hour,  it  will  not  stop  its  journey  to  compute  a  new  random 
position,  but  will  continue  on  its  present  coiuse  and  speed  until  it  reaches  the  waypoint  If  the  SSM  has  fired  all 
of  its  missiles,  it  will  commence  a  transit  to  its  assigned  missile  depot  to  replenish  its  missile  battery.  The  course 
and  speed  to  reach  the  depot  are  calculated  using  equations  (1)  and  (9).  Once  the  SSM  has  reached  the  missile 
depot,  it  will  loiter  for  R  mimrtes  (R  =  60),  during  which  time  period  the  missile  battery  will  be  loaded  with  its 
full  complement  of  missiles.  Once  the  reload  process  is  corr^lete,  a  new  waypoint  will  be  chosen  using 
equations  (3)  and  (4).  The  original  intent  was  for  the  SSM  laimcher  to  return  near  its  initial  position  and  then 
commence  its  normal  movemoit  pattern  The  computer  simulation  maintains  the  laimcher  near  the  reload  site 
for  the  remainder  of  the  simulatioa  This  is  an  incorrect  implementation  of  the  simulation  and  should  be 
corrected  in  fiiture  revisions.  The  algorithm  that  determines  SSM  movement  is  described  in  Appendix  C. 
c  United  States  and  Iranian  Reconnaissance  Aircraft 

The  recormaissance  aircraft  fly  in  straight  line  paths  between  three  waypoints  at  speed  Voir,  where 
Voir  is  a  normal  random  variable  with  mean  pvai  (nominal  operating  speed),  and  variance  conditioned  to  be 

positive.  The  random  qreed  accounts  for  deviations  in  course  and  speed  resulting  fi'om  operational  and 
equipment  variations.  An  independent  replication  of  is  drawn  each  time  the  reconnaissance  aircraft  departs 
a  waypoint  Every  recormaissance  aircraft  is  assigned  a  unique  search  pattern  in  the  data  base  which  is 
computed  with  respect  to  the  location  of  the  air  base  or  current  location  of  the  aircraft  carrier  on  which  it  is 
located.  Every  unit  has  a  position  which  is  defined  by  an  X  and  Y  coordinate.  A  new  position  can  be  dsfinmt  in 
terms  of  the  direction  and  distance  fi:om  a  current  X  coordinate  and  the  direction  and  distance  firom  a  current  Y 
coordinate.  The  following  procedure  defines  the  east  (Xdirection  =  1)  or  west  (Xdirection  =  -1)  direction  a  new 
X  coordinate  is  located  with  respect  to  the  current  X  coordinate  and  stores  it  in  the  variable  Xfocdon-  A  new  X 
coordinate  is  flien  computed  by  multipfying  the  Xfcecticn  by  a  specific  X  distance,  denoted  Xiavei_diaance,  and  adding 
this  product  to  the  current  X  coordinate.  A  similar  procedure  is  then  used  to  compute  the  Y  coordinate  using  the 
variables  Ydirection  (north  =  -1,  south  =  1)  and  Ytravel_distance.  The  travel  distances  are  the  distances  in  miles 
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that  an  aircraft  will  travel  in  the  x  and  y  directions  to  reach  its  next  waypoint.  These  distances  are  constant  and 
set  prior  to  the  commencement  of  a  simulation  run  in  the  initial  data  base 

A  reconnaissance  aircraft  flies  from  its  aircraft  carrier  or  air  base  to  its  first  waypoint  (X[l],  Y[l]) 
which  is  computed  using  the  carrier  or  air  base  as  the  current  position  and  the  variables  Ydiwikn, 

X[l]travd_distaiice,  and  Y[l]tavei_distance-  In  the  cusB  of  an  Iranian  reconnaissance  aircraft,  aU  variables  are  defined  in 
the  simulation  or  the  initial  data  base  and  do  not  change  during  the  simulation  For  U.  S.  reconnaissance 
aircraft,  the  only  variable  that  is  not  permanently  defined  is  the  Xdiecdoto  which  is  computed  to  be  the  same  east  or 
west  direction  as  the  carrier’s  present  course.  After  the  first  wtypoint  is  corr^uted,  the  second  waypoint  pf[2], 
Y(2])  is  computed  with  the  first  waypoint(X[l],Y[l])  defined  as  the  current  position,  and  the  variables  XdhKjicm, 
Ydipgctioib  ^^PltraBrei_dist3iicc?  and  The  same  ^CdirectiGn  and  ^^dirccdcn  Variables  are  used  ^br  this 

computation,  but  for  an  Iranian  aircraft,  the  product  of  the  and  X[2]ta,ei_distance  is  subtracted  finm  the 

current  X  coordinate  instead  of  added  as  done  in  computing  die  first  waypoint.  The  third  w^rpoint  is  not 
defined  until  the  aircraft  has  arrived  at  the  second  waypoint  The  third  waypoint(X[3],Y[3])  is  alw^s  set  equal 
to  the  X  and  Y  coordinates  of  the  air  base  or  aircraft  carrier  from  which  the  aircraft  started  its  fligbi  Since  an 
aircraft  carrier  is  constantly  changing  position,  a  U.  S  recormaissance  aircraft  is  required  to  compute  a  new  third 
w^point  (X[3],Y[3])  every  minute  until  it  reaches  the  location  of  the  carrier  and  completes  its  flight  Figure  4 
shows  an  exanqile  Iranian  reconnaissance  aircraft  fli^t  pattern,  which  is  indicated  by  the  solid  1iTip<;  and 
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detennined  using  equations  (10)  through  (15). 


Xdireaicn  =  1  (fly  east)  or  -1  (fly  west)  {set  fi'om  initial  data  base}) 

Ydrction  =  1  (ffy  south) 


Y[l]  Ypresentjocation  Ydirection  •  Y[l]tracvdl_daiaiicc 

(10) 

X[l]  Xpresent  location  Xdircction  •  X[l]trayel  distance 

(11) 

Y[2]  =  Y[1]  ■ 

(12) 

X[2]  =  X[l]  -  Xdirecdon  •  X[2]travcl_dist2nce 

(13) 

X[3]  =  Air  Base^q^rescntjoatm 

(14) 

Y[3]  Air  BSSC  ypresent_location 

(15) 

V^r  ^  (Mvairj^vair) 

(16) 

The  Iranian  aircraft  departs  fi-om  the  location  of  the  air  base  and  travels  in  a  istraight  line  to  its  first 
waypoint(X[l],Y[l]).  The  aircraft  then  computes  a  new  course  and  speed  to  reach  the  second  waypoint 
(X[2],Y[2])  and  commences  travel..  The  aircraft  then  flies  directly  back  to  the  air  base(X[3],Y[3]).  The  course 
and  speed  for  each  leg  of  an  Iranian  recoimaissance  aircraft  flight  pattern  are  computed  using  equations  (1)  and 
(16). 

The  two  waypoints  a  U.  S.  reconnaissance  aircraft  flies  through  are  computed  prior  to  launch,  and 
depend  on  the  present  location  of  the  aircraft  carrier.  Figure  5  provides  an  example  of  a  U.  S.  reconnaissance 
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aircraft  search  pattern,  which  is  indicated  by  the  sohd  lines  and  determined  using  equations  (10)  through  (12), 
and  (17)  through  (21).  The  definition  of  the  Xaitcdon  variable  determined  in  Equation  (17)  causes  a  U.  S. 
reoormaissance  aircraft  to  fly  a  pattern  away  from  the  carrier,  but  in  the  same  direction  of  the  carrier’s  motion 
east  or  west  This  allows  the  aircraft  to  detect  threats  that  are  located  in  the  area  of  the  carrier’s  intended 
movement 

A  U.  S.  aircraft  departs  fi'om  the  currait  location  of  the  carrier  and  travels  in  a  strai^t  line  to  its 
first  waypoint  (X[1],Y[1]).  The  aircraft  then  determines  a  new  course  and  speed  and  commences  travel  to  the 
second  w^point(X[2],Y[2]).  The  aircraft  then  returns  to  the  aircraft  carrier(X[3],Y[3]).  During  its  return 
flight  to  the  carrier,  the  aircraft  computes  a  new  course  and  speed  every  minute  to  account  for  the  movement  of 
the  aircraft  carrier.  The  course  and  speed  for  each  leg  of  a  U.  S.  aircraft  flight  ate  computed  using  equations  (1) 
and  (16). 


N 
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X[2I  travel_distance 
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’^ni-Y[ll] 


Yin  travel 
distance 


Figure  5.  United  States  Reconnaissance  Aircraft  Search  Pattern 
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Xdirecdon  =  1  (fly  east) 

0  <  CV„^  <71/2 

(17) 

-1  (fly  west) 

Ydirection  =  -1  (fly  nordi) 

7l/2<CVoouree  <51 

(18) 

Y[l]  ”  Yppes(ntj[ocation  Ydjfection  ® 

^[l]tiavel_distancc 

(10) 

X[l]  Xprcsent  location  Xdirgction  • 

^[l]travel_distancc 

(11) 

Y[2]  =  Y[1]  ” 

(12) 

X[2]-X[l]+Xar^^.X[2U, 

[^distance 

(19) 

X[3]  —  CVjgjTCsmtJocation 

(20) 

Y[3]  “  CVypresent^Iocation 

(21) 

Once  the  aircraft  reaches  the  carrier  or  air  base,  it  remains  on  the  surfece  for  a  certain  time  interval 
(loiter  time),  which  represents  the  refueling  and  maintenance  of  the  aircraft,  and  rest  for  the  crew.  The  loiter 
time  for  the  aircraft  is  tracked  %  the  same  method  used  for  naval  ships.  When  the  loiter  time  has  elapsed,  new 
wa5points  are  computed  based  on  the  location  of  the  air  base  or  present  location  of  the  aircraft  carrier,  and  a 
reconnaissance  aircraft  commences  flying  its  search  pattern. 

In  practice,  reconnaissance  airaafl  tend  to  follow  courses  that  are  judged  to  be  potentially 
informative.  The  specific  waypoints  defined  for  the  anafysis  performed  in  this  thesis  ensure  that  the  courses  are 
potentially  informative.  In  a  realistic  environment,  aircraft  cueing,  e.g.  by  an  external  (satellite)  or  shipboard 
sensor,  is  used  to  vector  aircraft  to  a  specific  locatiom  This  feature  is  not  represented  in  the  model,  but  is  a 
desired  improvement  A  reconnaissance  aircraft  also  needs  the  ability  to  alter  its  flight  path  if  a  hostile  aircraft 
or  missile  is  detected  near  the  plane.  In  this  model  the  recormaissance  aircraft  are  susceptible  to  attack  Ity  figtitpr 
aircraft  that  defend  the  CVBG  and  the  air  base.  Fighter  and  attack  aircraft  are  considered  we^x>ns  systems  for 
the  purpose  of  this  model  and  are  discussed  later  in  this  section.  The  algorithm  that  determines  reconnaissance 
aircraft  motion  is  described  in  Appendix  C. 

d.  Satellites 

Each  satellite  moves  between  two  waypoints  that  are  defined  in  the  initial  data  base.  The  course 
that  a  satellite  follows  is  computed  using  equation  (1),  and  the  speed  is  constant  (no  variance)  and  is  also  dpfinptd 
in  the  initial  data  base.  Each  satellite  travels  fi^om  its  initial  waypoint  to  the  other  waypoint.  Once  it  reaches  the 
new  location,  it  loiters  a  set  amount  of  time,  which  is  intended  to  represent  the  time  the  satellite  is  p  ^ing  over 
other  regions  of  the  earth  during  its  orbit.  Once  the  time  is  reached  for  the  satellite  to  commence  a  pass  over 
Iran,  the  satellite  travels  from  its  present  location  to  the  other  waypoint  This  .sitnnlates  a  back  and  forth  pass 
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over  the  same  region  which  is  not  correct,  but  is  done  as  a  programming  simplification.  For  the  purpose  of  this 
stmfy,  the  wi^points  are  defined  so  that  the  satdlites  fly  over  areas  of  Iran  where  SSM  launchers  and  air  bases 
are  located.  The  initial  data  base  provides  the  U.  S.  forces  with  a  pair  of  satellites.  One  of  the  satellites  flies  over 
Iran  south  to  north,  and  one  flies  east  to  west 

2.  Search  and  Detection 

a.  Assumptions 

1 .  The  ESM  and  ELINT  sensors  are  only  able  to  detect  and  localize  targets  that  are  emitting. 

2.  The  air  search  sensor  only  detects  airplanes. 

3.  The  sut&ce  search  sensor  detects  ships,  SSM  launchers,  air  bases,  and  radar  sites. 

4.  This  model  does  not  take  into  account  the  altitude  of  the  aircraft. 

5.  The  U.  S.  satdlites  are  undetectable  by  the  Iranian  sensors. 

6.  All  detection  information  is  immediately  shared  among  all  friendly  units  (i.e.  communications 
links  work  perfectly:  there  is  no  congestion  dd^). 

7.  The  U.  S.  assets  always  emit  (radar  in  operation),  and  therefore  are  susceptible  to  detection  by 
the  Iranian  ESM  sensors.  Iranian  air  bases  and  long  range  radar  sites  always  emit,  and  therefore  are  susceptible 
to  detection  by  U.  S.  ESM  sensors. 

All  of  these  assurtqrtions  are  candidates  for  modification  in  later  revisions  of  the  model. 

b.  Search  and  Detection  Theoiy 

Each  sensor  system  has  a  maximum  detection  range  expressed  in  the  depth  of  surrounding  grids 
throughout  which  there  exists  a  positive  probabihly  of  detecting  an  enen^f  unit  That  is,  if  the  range  is  one  grid, 
detections  occur  within  the  grid  at  which  the  sensor  is  located,  and  all  adjacent  grids,  (assuming  a  target 
susceptible  to  detection  is  present).  The  detection  range  assumes  that  the  sensor  is  located  in  the  center  of  the 
grid. 

Each  sensor  has  a  probability  of  detection  given  a  unit  is  within  its  detection  range.  Upon  detection 
of  an  enenty  unit,  the  perceived  location  of  the  eneny  unit  is  reported  by  the  sensor,  which  is  the  actual  location 
(ground  truth)  ofi&et  by  errors  in  the  x  and  y  direction.  This  perceived  location  is  distributed  as  a  circular  normal 
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distribution  with  parameters  p*,  py,  and  ajenscr-  The  values  of  p*  and  py  are  the  ground  truth  x  and  y  coordinates 
of  the  target,  and  asnsor  which  is  a  sensor  error  term  accounting  for  the  accuracy  of  the  physical  sensor  and 
human  operator  of  the  systeim  The  information  derived  from  each  sensor  is  assimilated  with  previous  detections 
of  the  same  unit  to  refine  the  inteUigence  estimate  of  enemy  capabilities. 

c.  Search  and  Detection  Lnplementation 

On  every  time  stq)  (one  minute),  each  unit  searches  for  enemy  units  using  all  three  of  its  sensor 
systems  in  order:  surface,  air,  and  ESM  For  each  sensor,  the  model  determines  the  grids  surrounding  the  unit 
where  detections  occur  given  a  target  is  present,  and  then  sequential^  searches  each  grid.  Within  each 
grid,  the  model  individually  examines  each  enemy  unit  located  on  the  grid’s  list  of  eneirq^  units  and  determines  if 
it  is  eligible  for  detection  the  sensor  (i.e.  aircraft  by  air  search  sensor,  unit  that  is  emitting  by  ESM  sensor, 
etc.).  Given  it  is  a  candidate,  an  indqjendent  uniform  random  variable  with  range  [0,1]  is  drawn,  fr  the  random 
variable,  which  is  denoted  Haetecuon  ,  is  less  than  the  probabihty  that  the  sensor  detects  a  target  given  it  is  within 
its  search  area,  a  detection  occurs,  ff  a  target  has  been  detected,  the  sensor  reports  the  perceived  location  of  the 
target  as  described  previously.  Equations  (22)  and  (23)  illustrate  this  approach. 

^detection  ^  1^(0, 1) 

Detection  occurs  if  ^  Pdettsemor] 

No  detection  ^  ^detection  ^  Pdet[semor]  (^2) 

Perceived  Location  ~  A^(p*,  py,  ct^)  (23) 

This  simulation  assumes  that  all  sensors  perfectly  determine  the  identity  of  a  target  (unit 
identification  number)  and  there  are  no  felse  detections, 
d.  Intelligence  Estimate 

As  discussed  earlier,  enemy  unit  detections  are  maintained  on  lists  that  are  resident  on  the  grid 
system.  Each  force  has  700  fists  of  detected  enenty  targets  (a  list  for  every  grid).  Some  of  the  fists  mntain 
numerous  targets,  while  others  remain  empty  the  entire  simulatioa  These  fists  identify  the  units  by  a  specific 
unit  number,  the  current  estimate  of  its  Cartesian  location,  the  number  of  detections  of  this  enemy  unit,  and  the 
time  of  the  latest  detection.  Each  time  a  sensor  makes  a  new  contact,  i.e.  achieves  a  detection,  it  HatfirminRs  the 
grid  location  corresponding  to  the  perceived  location  reported  by  the  sensor.  The  fist  located  on  that  grid  is  then 
scanned  to  locate  previous  observations  on  this  particular  contact.  If  there  is  a  contact  history,  an  average 
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updating  process  is  used  to  refine  the  contact  location,  otherwise  the  new  contact  is  added  to  the  list  as  a  different 
entity.  The  perceived  location  of  a  target  may  be  on  a  grid  where  the  target  is  not  actually  located.  If  two 
separate  detections  locate  the  same  unit,  but  the  perceived  locations  are  significantly  different  and  place  the  units 
on  two  different  grids,  there  will  be  an  entry  for  the  unit  on  each  grid.  In  a  real  scenario,  a  person  analyzing  this 
information  may  believe  there  are  two  separate  contacts.  However,  the  weapon  targeting  model  of  this 
simulation  does  not  allow  the  firing  of  two  separate  sorties  of  weapons  at  the  same  target  based  on  entries  on  two 
different  grid  hsts. 

Equations  (24)  through  (26)  describe  the  information  fijsion  process.  The  variables 

Yperc€ption_oid,  ^ud  uumbeT  of  observations  are  defined  for  each  unit  that  is  located  on  a  grid’s  list  If  a  unit  is 

located  on  two  sq)arate  grid  lists,  it  has  a  set  of  three  variables  for  each  list.  The  information  contained  within 

those  two  sets  of  variables  which  are  defined  on  different  numbered  grids  is  never  combined  to  provide  a  more 

precise  estimate  of  a  unit’s  actual  location.  The  'Xj,accpikm_oid  and  variables  are  the  average  X  and  Y 

perceived  locations  over  all  detections  of  a  unit  which  place  it  on  a  particular  grid,  but  does  not  include  the  most 

recent  detection.  The  number  of  observations  is  the  number  of  times  the  unit  was  reported  on  this  grid  prior  to 

the  most  recent  detection.  The  Xpe^^  new  and  Ypeiception^ncw  variables  are  the  average  perceived  X  and  Y 

location  of  the  unit  after  th^  have  been  averaged  with  the  most  recent  detectioa  The  nurriber  of  observations  is 

incremented  to  indicate  incorporation  of  this  new  detection  in  computing  the  unit’s  average  perceived  location. 

^^jerceptiorLjriew  ”  (Observations  •  ^^pcn:;qptkm_old  1) 

Yperoqptioiujiew  (Observatioiis  •  Ypcjcqjtion_oid  ^pcrcqjtkxujiew)/ (Observations  1)  (25) 

Observations  =  Observations  +  1  (26) 

The  model  deletes  a  contact  position  on  a  specific  grid  if  the  contact  has  not  been  reported  as 

detected  on  that  grid  by  at  least  one  sensor  among  the  entire  force,  during  the  latest  T  ipdate  cycles  (  T  =  15 

minutes).  This  is  done  to  eliminate  contacts  on  the  targets  that  have  changed  grids,  were  reported  incorrectly  as 

on  a  grid,  or  have  been  killed.  The  mcxlel  does  not  allow  a  sensor  to  detect  a  target  that  has  been  killed. 

3.  Weapons  Engagements 
a.  Engagement  Theory 
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Each  weapon  system  is  assigned  probability  of  defeat  by  the  target’s  defensive  mechanisms 
(Pdefeat)  and  3  probability  of  hit  (phit)  against  a  particular  target  type  by  the  user  prior  to  a  simulation  run. 
The  defensive  capability  of  each  unit  against  a  type  of  weapon  is  aggregated  into  a  single  number  (pdefeat ) 
that  represents  the  probability  that  the  unit  will  defeat  an  incoming  weapon  prior  to  the  weapon  releasing 
its  ordnance.  Each  unit  has  is  assigned  a  different  probability  of  defeating  an  incoming  weapon  for  each 
type  of  weapon  that  it  may  encoimter.  Defensive  weapons  systems  such  as  Iranian  land-based  surface-to- 
air  missile  (SAM)  sites,  which  protect  high-value  targets  (air  bases,  radar  sites)  are  included  in  this 
aggregation.  The  only  time  an  incoming  weapon  is  susceptible  to  destruction  is  when  it  reaches  its  target 

In  a  real  conflict,  the  probability  that  a  unit  can  defend  itself  is  affected  by  previous  damage 
incurred  by  the  unit.  The  computer  implementation  of  this  simulation  model  assumes  a  unit’s  defensive 
capabilities  are  fully  operational  until  it  is  destroyed.  This  model  could  be  improved  by  determining  the 
strength  of  a  unit’s  defensive  capabilities  based  on  the  number  of  hits  it  has  previously  received  and 
decreasing  the  probability  that  its  defensive  mechanisms  defeat  an  incoming  weapon  (parfea). 

Additionally,  in  the  case  of  naval  ships  operating  as  a  task  group,  the  damage  to  the  group  as  a  combined 
entity,  affects  the  ability  of  each  unit  to  defeat  incoming  weapons.  Under  wartime  conditions,  the 
Ticonderoga  class  cruisers  are  responsible  for  the  air  defense  of  the  carrier  battle  group.  If  the  cruiser  is 
damaged  then  the  ability  of  all  other  units  to  defend  themselves  against  air  attacks  is  diminished.  This 
feature  should  also  be  addressed  in  future  improvements  of  the  model. 

The  determination  of  whether  a  weapon  scores  a  hit  against  a  target  depends  on  its 
probability  of  hit  (p«,)  against  that  target  type  and  the  accuracy  of  the  targeting  information  at  the  time  of 
fire.  In  this  model  the  probability  of  hit  represents  all  weapon  delivery  error  (wind,  targeting  computer, 
etc.).  The  probability  of  hit  is  reduced  due  to  target  location  errors.  A  targeting  model,  which  is  discussed 
later,  designates  a  sortie,  which  is  a  specific  weapon  type  and  the  number  of  weapons  of  this  type  to  fire  at 
a  target.  The  sortie  is  targeted  at  the  location  where  the  aggressor  believes  the  enemy  is  positioned 
(XpeK;eivcd,Yp„ceived)-  The  (Ustance  between  this  position  and  the  ground  truth  location  of  the  target  (Xactuai, 
Yactuai)  is  used  to  compute  a  multiplier  (£>«,),  which  degrades  the  probability  that  a  weapon  scores  a 
hit  against  the  target  (Phit)Phit.  This  simple  approach,  which  is  shown  in  equations  (27)  and  (28), 
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describes  how  the  accuracy  of  the  perceived  target’s  position  influences  the  probability  that  the  weapon 
scores  a  hit  against  a  target.  Since  it  computes  the  aimpoint  error  at  the  time  of  fire,  it  does  not  enable  the 
weapon  to  receive  mid-course  guidance,  or  model  weapons  that  can  manually  adjust  their  course,  such  as 
piloted  aircraft  and  radar-seeking  seeking  missiles. 

■>j(Xactual  -  Xperceived)^  +  (Yactual  -  Yperceived)^ 

AimptError  = - -  (27) 


_  —AoTpEnvr 

LMt  —  c  if  Aimpoint  Error  <  10  miles 

.2  if  Aimpoint  Error  >  10  miles  (28) 

Until  the  sortie  reaches  the  target  and  the  result  of  the  strike  is  determined,  no  other  sorties 
can  be  fired  at  the  target.  This  does  not  allow  for  a  force  to  utilize  the  tactic  of  saturating  the  enemy  with 
incoming  weapons.  When  employed,  this  tactic  increases  the  probability  of  some  of  the  weapons  evading 
a  target’s  defensive  mechanisms  and  striking  the  target.  The  weapon  selection  decision  module  does  not 
select  a  mixed  group  of  weapon  types  to  fire  at  a  single  target. 

The  simulation  time  at  which  the  sortie  is  to  arrive  at  the  target’s  location  is  computed  using 
the  distance  between  the  sortie’s  point  of  origin(Xfcing^u„it,  YSring^umt)  and  the  ground  truth  location  of  the 
target(Xactuai,  Yactuai)>  the  nominal  (operational)  speed  of  the  weapon  type  (Weapon  speed)  which  comprises 
the  sortie,  and  a  preparation  time.  The  preparation  time  represents  the  time  required  to  compute  a  firing 
solution  and  prepare  weapons.  This  time  is  referred  to  as  the  on  target  time  (OTT)  and  is  computed  using 
the  equation  (29). 


OTT  =  'Present  Time  +  Prep  Time  + 


^(Xfiring_ 


unit  -  Xactual)^  +  (Jfiring _unit  -  Yactual)^ 


(29) 


Weapon  speed 

This  assumes  that  the  target  does  not  change  location  diuing  the  sortie  transit  time,  which  is  not  realistic 
but  is  done  for  simplicity  of  modeling.  The  flight  path  of  missiles  and  fighter  and  attack  aircraft  are  not 
placed  on  the  grid  system.  This  is  a  simplification  of  the  model  and  ignores  the  intelligence  gained  by 
placing  it  on  the  grid.  An  attacking  aircraft  provides  additional  reconnaissance  during  its  flight,  and  the 
direction  of  an  incoming  strike  gives  the  target  a  rough  idea  of  the  location  of  the  platform  that  launched 
the  strike. 
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Once  the  simulation  time  is  reached  that  indicates  a  sortie  has  arrived  at  the  target’s  location, 
an  adjudication  of  the  attack  is  conducted.  During  conflict  adjudication,  the  model  first  determines  how 
many  of  the  weapons  in  the  sortie  evaded  the  target’s  defensive  mechanisms  During  a  real  strike, 

the  direction  and  time  that  multiple  weapons  arrive  at  a  target  are  not  always  the  same.  For  example,  a 
group  of  attacking  aircraft  may  separate  on  their  final  approach  to  a  target  and  attack  from  different 
directions  and  in  waves  to  confuse  and  complicate  the  defensive  response  of  the  target.  As  a  result,  in  this 
simulation,  the  determination  of  whether  a  single  weapon  of  a  sortie  evaded  the  target’s  defensive 
mechanisms  is  considered  independently  of  the  sortie’s  other  weapons.  Hence,  the  nmnber  of  weapons 
that  evade  the  target’s  defensive  mechanisms  has  a  binomial  distribution  with  probability  of  success  equal 
to  p defeat  and  number  of  trials  equal  to  the  number  of  weapons  in  the  sortie.  The  number  of  hits  scored 
against  the  target  by  those  weapons  that  get  through  the  defensive  mechanisms  is  determined.  The 
number  of  hits  has  a  binomial  distribution  with  probability  of  success  equal  to  Phit«Dhitand  number  of 
trials  equal  to  the  number  of  weapons  that  survived  the  target’s  defensive  mechanigms 

Each  unit  has  a  counter  that  keeps  track  of  the  total  number  of  hits  it  has  sustained  during  the 
course  of  the  simulation.  The  total  number  of  hits  the  target  received  during  this  strike  is  added  to  the 
counter,  and  a  check  is  performed  to  ensure  the  target  has  not  exceeded  its  maximum  number  of  hits.  A 
target  that  has  sustained  at  least  its  maximum  number  of  hits  is  considered  destroyed.  If  the  target  is 
destroyed,  its  type  is  changed  to  indicate  it  no  longer  exists.  This  means  that  its  sensor,  weapon  and 
propulsion  (movement)  ^sterns  are  rendered  non-fimctional  and  provide  no  input  for  its  respective  force 
the  remainder  of  the  simulation.  If  the  attacking  sortie  consists  of  aircraft,  the  number  of  the  aircraft  that 
are  not  destroyed  by  the  target’s  defensive  mechanisms  are  returned  as  available  assets  onboard  the 
carrier  or  air  base  after  a  time  delay  which  accounts  for  their  return  flight, 
b.  Engagement  Implementation 

(1)  Overview.  The  model  maintains  two  engagement  lists  for  conflicts  between  the  two  forces, 
one  for  the  United  States  force  striking  Iranian  units,  and  another  for  an  Iranian  force  striking  U.  S.  units. 
During  every  update  cycle,  the  targeting  module  examines  every  grid’s  inteUigence  estimate  list  sequentially. 
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starting  at  grid  one  and  finishing  at  grid  700.  Three  steps  are  performed  in  this  examination  phase  for  each 
enemy  unit  on  a  grid  list  that  is  not  currently  targeted  (i.e.  a  weapon  is  not  enroute  to  the  target)  or  dead,  which 
determines  the  action  to  be  taken  against  an  enemy  unit  The  first  step  classifies  the  threat  level  of  the  target, 
and  determines  the  closest  missile-capable  and  aircraft-capable  fiiendfy  unit  that  has  weapons  available  to 
engage  the  target  In  the  next  step,  one  of  the  two  units  is  assigned  to  prosecute  the  target  using  a  sortie.  The 
decision  of  which  unit  to  assign  is  based  on  the  target  type  and  is  described  by  the  targeting  algorithms  included 
in  Appendix  C.  Finally,  the  time  for  the  sortie  to  reach  the  enemy  unit  (time  of  conflict  adjudication)  is 
computed,  and  the  conflict  placed  on  the  engagement  list  It  is  irrqxjrtant  to  note  that  these  three  steps  are 
performed  on  each  target  individually  before  addressing  the  next  target  Also  all  contacts  on  the  grid  one  list  are 
examined  and  decisions  made  on  we^ns  to  fire  at  all  of  the  targets  listed  as  on  grid  one  prior  to  examining 
targets  that  are  listed  as  on  grid  two,  and  so  forth.  A  result  of  this  process  all  targets  are  prosecuted,  whether 
high  or  low  threat  until  we^ns  supplies  are  dqjleted  starting  at  the  lowest  numbered  grid.  If  a  high  threat 
target  is  located  on  a  grid  with  a  high  grid  nuniber,  the  possibility  exists  that  weapons  will  not  be  available  to 
prosecute  the  target  In  future  revisions  of  this  model  it  is  recommended  that  a  mote  realistic  weapons  allocation 
procedure  be  implemented.  Once  all  targets  which  exist  in  the  entire  grid  system  have  been  placed  on  this  hst,  a 
determination  on  which  targets  to  prosecute  can  be  made  on  the  basis  of  threat  level. 

(2)  Threat  Determination  and  Closest-Unit  Selection.  During  this  stq>,  the  enemy  imit  is 
classified  as  a  high,  medium,  or  low  threat,  m  accordance  with  the  algorithm  located  in  Appendix  A.  The 
algorithm  first  computes  the  distances  between  the  perceived  enemy  position  (Xpo„ived,  Ypocoveci)  and  every  one  of 
the  fiiendly  units.  If  any  of  the  distances  is  smaller  than  the  effective  we^n  range  (of  a  we^n  that  can  attack 
a  given  fiiendly  unit)  of  the  enemy  unit,  it  is  classified  as  a  high  threat;  otherwise  it  is  a  low  threat.  There  are 
two  exceptions;  an  Iranian  long-range  radar  site  is  classified  as  a  medium  threat,  and  an  Iranian  air  base  which 
is  located  where  its  aircraft  are  not  within  the  range  that  they  can  strike  any  mnltinatinnal  force  or  U.  S.  naval 
ships  is  classified  a  medium  threat  instead  of  a  low  threat.  During  this  step  the  distance  information  is  used  to 
identify  the  closest  fiiendly  missile-capable  unit  and  aircraft-c^able  unit  that  has  both  weapons  available 
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(#weapoiis  >  imnimum  reserve  number  of  a  we^n  that  must  remain  onboard  a  unit  at  all  times),  and  is  located 
within  effective  weapon’s  strike  range  of  the  enemy  unit  for  possible  target  prosecution. 

(3)  Assign  Unit  and  AUocate  Weapons.  There  is  a  sqjarate  decision  module  for  each  force 
(U.  S.  and  Iranian)  which  determines  which  umt  prosecutes  a  given  enemy  unit,  and  the  numbers  and  type  of 
we^n  to  utilize.  Both  of  these  decision  modules  are  somewhat  arbitrary  and  do  not  follow  any  known  mihtaiy 
guidance.  Alternative  decision  processes  can  be  reasonabfy  contemplated  and  studied  within  the  structure  of 
this  simulation.  The  decision  module  bases  its  unit  selection  on  the  level  of  threat,  and  the  type  of  enemy  unit 
For  example,  a  high  threat  Iranian  air  base  will  be  assigned  to  an  aircraft  carrier  using  numerous  attack  aircraft, 
whereas,  a  low-threat  Iranian  air  base  would  be  assigned  to  the  closer  of  a  Tomahawk-c^able  warship  or  a 
carrier  dispatching  a  smaller  number  of  attack  aircraft.  The  two  algorithms  that  determine  U.  S.  and  Iranian 
weapons  allocation  are  located  in  Appendix  C.  Once  the  unit  and  weapon  assets  are  determined,  the  counters 
that  track  umt  weapon  availability  and  the  force  weapon  use  are  updated  accordingly.  In  the  case  of  attack  or 
fighter  aircraft,  the  carrier  or  Iranian  air  base  will  not  gain  the  availability  of  those  assets  until  the  strike  is 
completed  (conflict  adjudicated)  and  aircraft  £ty  back  to  their  point  of  origination. 

(4)  Engagement  List  There  is  a  separate  engagement  list  for  each  force  (U.  S.  attacking 
Iranian  and  its  counterpart).  The  time  the  weapon  is  on  target  and  airrpoint  error  degradation  {Dhit)  is  computed 
as  stated  previously  in  the  engagement  theory  discussion.  The  engagement  is  then  added  to  the  engagement  list. 
Each  item  on  the  list  identifies  the  time  the  weapon  reaches  the  target,  the  unit  that  fired  the  weapon,  the  enerity 
umt  that  is  under  attack  from  this  we^n,  the  type  of  weapon,  the  number  of  weapons  fired,  and  the  aimpoint 
error.  Once  a  conflict  is  adjudicated  it  is  removed  from  its  engagement  list.  This  entire  targeting  process  may  be 
improved  by  developing  more  detailed  threat  prioritization  rules,  and  refining  the  rules  to  assign  weapons  to  a 
target  based  on  current  militaiy  doctrine. 
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c.  Conflict  Adjudication 

(1)  Defensive  Weapons.  When  a  sortie  reaches  its  intended  target,  independent  Bernoulli 
trials  are  simulated  for  each  weapon  of  the  sortie  in  accordance  with  equation  (30)  to  determine  if  the 
number  of  weapons  that  survive  the  target’s  defensive  mechanisms.  Hoi^eat  is  a  draw  from  a  ntiifnnu 
distribution  over  [0,1], 

H,efea,  ~  U(0,l) 

Weapon  survived  if  Hje/eat  >  pdcfeat 

Weapon  destroyed  if/f*f«<,r  ^Pdefeat  (30) 

The  independent  assumption  is  used  since  weapons  fired  at  a  single  target  are  often  fired  at 
different  times  and  may  approach  the  target  from  varied  directions.  If  the  weapon  is  not  destroyed,  it  then 
has  a  probability  of  scoring  a  hit  against  the  target. 

This  representation  of  the  defensive  weapon  systems  is  simplistic  in  nature,  and  further 
improvements  should  model  significant  defensive  systems  as  separate  entities.  Another  improvement  to 
the  model  would  allow  the  engagement  of  incoming  weapons  prior  to  reaching  their  target.  In  a  red 
conflict,  U.  S.  strike  aircraft  would  have  to  avoid  Iranian  SAM  sites  located  just  inside  Iran’s  borders,  and 
Iranian  strike  aircraft  could  be  engaged  by  U.  S.  combat  air  patrol  aircraft  and  surface-to-air  missiles  from 
U.  S.  warships  long  before  they  reach  their  intended  target. 

(2)  Offensive  Weapons.  Once  the  number  of  weapons  that  survive  the  target’s  defense 
mechanisms  is  computed,  the  munber  of  hits  scored  against  the  target  is  determined  in  accordance  with 
equation  (3 1).  Independent  Bernoulli  trials  are  conducted  for  each  weapon  striking  the  target  to 
determine  how  many  hits  the  target  sustained.  Hku  is  a  draw  from  a  uniform  distribution  over  [0,1]. 

Hut  ~  U(0,1) 

Weapon  hit  target  if/f«t  <Phit*I>«/ 

Weapon  missed  target  if  Hut  >  Phit  •  Dut  (31). 

Whenever  an  aircraft  carrier  or  air  base  is  a  target,  it  is  assumed  that  both  the  area  where 
planes  are  launched  (runway)  and  the  location  of  planes  on  the  surface  are  targeted.  This  model  does  not 
possess  the  detail  required  to  adequately  account  for  the  attrition  of  planes  under  these  circumstances. 
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however,  a  simple  approach  is  implemented  which  recognizes  the  attrition  of  planes.  A  number  is 
defined,  Acy,capon  type,,  which  attempts  to  represent  both  the  location  of  aircraft  throughout  an  air  field,  or 
aircraft  carrier,  and  the  lethal  area  of  an  incoming  missile  or  bomb.  These  munbers  for  U.  S.  weapons 
against  an  Iranian  air  base  and  Iranian  weapons  against  a  carrier  are  defined  at  the  beginning  of  the 
simulation.  This  niunber  does  not  address  airplanes  that  are  grouped  together  in  the  same  area,  so  an 
additional  random  number  is  computed  called  the  effectiveness  (E)  of  the  strike  to  identify  this  situation; 
£■  is  a  uniform  random  variable  over  [0,2].  The  range  of  the  uniform  random  variable  was  arbitrarily  set 
between  zero  and  two  for  analysis  purposes.  The  total  number  of  aircraft  killed  is  computed  using 
equation  (33).  This  equation  uses  the  number  of  aircraft  that  are  located  on  the  air  base  or  carrier  at  the 
time  the  conflict  is  adjudicated.  Since  both  attack  and  fighter  aircraft  may  be  destroyed  at  a  location,  the 
model  divides  the  number  of  aircraft  killed  in  proportion  to  the  number  of  each  aircraft  that  is  located  at 
the  base  or  on  the  carrier.  If  the  number  of  fighter  or  attack  aircraft  killed  exceeds  the  number  of  aircraft 
at  a  location,  the  maximiun  number  of  that  kind  of  aircraft  is  destroyed.  This  approach  is  illustrated 
using  equations  (32)  through  (35). 


E  ~  £7(0,2) 

(32) 

#  Aircraft  Killed  =ifhits  •  Acweapontype  •  #  Aircraft  •  E 

(33) 

r-  .  re  ^ FighterAircraft^ 

Fighter  Killed  =  - 

^HAircraftKilled 

(34) 

V  #  Aircraft  ) 

f  H  AttackAircraft^ 
Attack  Aircraft  =  - 

\#AircraftKilled 

(35) 

V  §  Aircraft  J 
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IV.  SIMULATION 


A.  OVERVIEW 

The  simulation  consists  of  the  model,  variables  that  record  aggregate  weapons  and  unit  inthrmatirtn 
and  decision  modules  that  force  the  units  within  the  model  to  interact  over  the  course  of  a  multiple-day  battle. 
The  simulation  occurs  in  cyclic  ipJates;  it  is  performed  at  one-minute  intervals  of  cycle  lengths,  but  is  easily 
modified  to  change  the  length  of  cycle.  The  simulation  is  programmed  in  the  C++  languagp:  to  perform 
quantitative  analysis. 

There  are  three  different  sources  of  mformation  which  affect  the  results  of  the  simnlatinn  The  three 
sources  are  the  initial  database  (geogr^hic  and  force  stractures),  the  decision  modules  that  govern  force 
interactions,  and  the  decision  parameters  utilized  m  the  decision  modules.  AH  three  of  these  are  fixed  prior  to 
conducting  a  simulation  run.  This  thesis  concentrates  on  sturfying  the  effects  of  adjusting  the  force  structure  and 
characteristics  of  individual  units  in  the  initial  dataha.«ie 

The  output  of  the  simulation  consists  of  three  files  that  track  specific  parameters  over  the  course  of  the 
simulatiotL  The  output  files  are  used  to  keep  track  of  the  locations  of  all  units  during  the  simulation, 
damage  incurred  by  individual  units,  the  total  number  of  weapons  fired,  weapons  destroyed,  and 
casualties  incurred  by  all  forces  at  the  end  of  the  two-day  battle. 

B.  UPDATE  CYCLE 

The  update  cycle  is  diagrammed  in  Figure  6.  The  initialization  phase,  which  loads  the 
geographic,  U.  S.,  and  Iranian  force  structure  database,  and  initializes  all  counters,  is  conducted  at  the 
beginning  of  a  simulation  run.  During  a  typical  cycle,  each  unit  loiters  or  continues  its  transit  towards  its  next 
waypoint  Each  unit  then  uses  all  ofits  available  sensors  to  detect  enemy  contacts.  After  each  force  processes  all 
the  information  received  from  its  sensors,  a  refined  intelligence  estimate  is  drawn.  Farh  force  then  Hpfprminpt 
the  enen^r  units  that  it  considers  threats,  and  commits  units  to  use  their  weapon  systems  to  strike  enen^^  units.  If 
any  weapons  have  reached  their  intended  target,  the  resulting  conflict  is  adjudicated.  If  the  time  is  reached  for 
an  aircraft  to  return  to  its  point  of  origin,  the  number  of  available  air  assets  for  its  carrier  or  air  base  is  updated  to 
reflect  its  return.  The  clock  is  then  incremented  and  the  process  continues. 
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figure  6.  Simulation  Update  Cycle 

1.  Initialization 

The  initialization  phase  of  the  simulation  loads  the  grid  infonnation  onto  the  map  from  a  data  file.  The 
simulation  then  reads  the  list  of  units  from  a  data  file  into  the  array.  The  simulation  then  examines  the  initial 
location  of  each  unit  to  determine  the  grid  where  each  unit  is  located,  and  adds  each  unit  onto  the  list  for  its 
respective  grid.  Initially,  all  units  are  stationary  and  will  change  position  according  to  the  movement  algorithms. 
There  is  luesently  no  initial  intelligence  estimate.  A  separate  data  file  (resulting  from  random  initial  sightings 
of  enemy  units)  may  be  incorporated  later  as  an  initial  intelligence  estimate. 

2.  Qock 

The  simulation  performs  an  update  cycle  at  one-minute  intervals  of  real  time.  This  implies  that  each 
unit’s  position  will  change  using  the  distance  it  covered  over  the  last  minute  of  time.  The  one-minute  interval  is 
chosen  to  accommodate  aircraft  and  missiles  which  move  at  high  speed. 

3.  Movement 

This  phase  of  the  simulation  determines  whether  a  unit  is  in  motion,  or  should  start  movement  If  it  is 
in  motion  it  computes  its  present  location  and  checks  to  see  if  it  moved  to  a  new  grid.  If  it  has  moved  to  a  new 
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grid  it  updates  the  associated  grids  to  reflect  this  change.  If  motion  starts  on  this  time  step,  it  updates  the  unit  to 
show  that  it  is  in  motion  and  determines  the  time  at  which  the  unit  will  reach  its  next  ws^point  All  motion  is 
conducted  according  to  the  algorithms  located  in  Appendix  A.  The  movement  of  all  units  is  updated  during 
every  cycle  which  equates  to  one  minute  of  scenario  time. 

4.  Sensor  Update 

Every  sensor  in  the  simulation  examines  all  the  grids  within  its  maximum  detection  rangft  during  each 
time  step.  The  intelhgence  estimate  is  then  updated  accordingly,  ffthere  is  no  detection  ofa  contact  for  the  last 
T  minutes  of  simulation  time,  the  contact  will  be  dropped  from  the  intelligence  estimate.  The  satpillitp  assets, 
which  are  an  important  &ctor  in  this  simulation,  have  a  behavioral  artifidahty:  each  scans  over  a  collection  of 
grids  (a  region  of  the  Iranian  country),  but  during  the  next  pass  of  the  satellite  over  the  same  region,  it  q~.ane  the 
grids  in  reverse  order.  This  is  a  computer  artificiality  because  of  the  method  used  to  simulate  traveling  between 
waypoints. 

5.  Enemy  Targeting 

The  model  generates  the  two  engagement  lists  and  rqxlates  global  variables  which  track  total  force 
we^n  use.  Whenever  an  enenty  unit  is  targeted,  the  targeted  variable  of  the  unit  is  flagged  to  indtraf^  a  weapon 
is  enroute,  which  prevents  multiple  weapons  fired  at  the  unit  at  the  same  time.  The  simulatinn  does  not  model 
the  situation  in  which  commanders  desire  to  saturate  one  target  with  multiple  types  and  numbers  of  weapons. 

6.  Conflict  Adjudication 

During  this  step,  the  simulation  examines  each  conflict  on  the  enpgement  list  and  determines  if  the 
we^n  is  now  over  the  target.  All  the  engagements  that  meet  this  criterion  are  removed  fiom  the  ftngagPTnpnt 
hst  and  adjudicated.  After  the  conflict  is  adjudicated,  the  model  ipdates  variables  to  indicate  weapons  and  units 
destroyed,  and  places  any  surviving  attacking  aircraft  on  the  time  delay  list  If  a  unit  survived  the  its 

targeted  variable  is  unflagged,  and  it  may  now  be  fired  at  again. 

After  all  the  conflrcts  for  a  given  time  stq)  have  been  adjudicated,  the  simulation  examines  the  hst  of 
aircraft  that  are  flying  back  to  aircraft  carriers  and  air  bases.  All  attack  and  fighter  aircraft  that  have  completed 
their  return  flight  during  this  time  step  are  added  as  available  assets  to  th^r  respective  carrier  or  air  base. 
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C  SOURCES  OF  INFORMATION 


As  previously  discussed,  the  simulation  is  influenced  Ity  three  diJBFerent  sources.  The  first  is  the 
algorithms  and  their  decision  parameters  (nurnbers  of  asset  to  use  in  a  strike  attack,  the  number  of  weapons  a 
unit  must  maintain  in  reserve,  etc.)  that  define  the  decision  modules  utilized  in  the  simulation  They  specify  the 
doctrine  under  which  the  battle  is  conducted,  wiiich  remains  unchanged  for  the  anafysis  of  this  thesis.  The 
algorithms  that  were  used  to  realize  the  simulation  as  a  conqjuter  program  are  included  in  Appendix  C.  The 
second  influence  is  the  list  of  constants  related  to  sensor  and  weapons  effectiveness  probabilities  that  are  used  in 
the  decision  modules.  The  constants  specify  probabilities  (?<!««»,  Pdeta,  and  Phs),  sensor  accuracy,  weapons 
speed,  and  lethality  of  weapons.  These  parameters  are  constant  for  a  given  set  of  simulation  runs  and  are  listed 
Appendix  D.  The  third  source,  is  the  initial  characteristics  of  each  unit,  and  the  numbers  of  units  that  are  loaded 
into  the  simuladorL  The  initial  characteristics  include  weapon  fypes,  nunibers,  and  ranges,  unit  operating 
speeds;  intended  track;  number  of  hits  a  unit  can  withstand  and  remain  fimctional;  sensor  range;  and  omiagiong 
control  status.  This  information  is  maintained  in  the  initial  force  database  and  is  the  easiest  source  of 
information  to  vary. 

D.  COMPUTER  MODEL 

The  simulation  is  coded  using  an  object-orientated  approach  in  the  programming  language  C-h-  which 
is  included  as  Appendix  E.  Mar^^  of  the  programming  techniques  and  data  structures  were  drawn  from  C++ 
How  to  Program  (Deitel,  1994).  The  data  structure  which  stores  the  characteristic  of  each  grid  and  unit  defined 
by  the  initial  database,  and  the  decision  module  algorithms  are  hard  coded  into  the  simulation  The  parameters 
which  are  utilized  by  the  algorithms  are  also  hard  coded  into  the  simulation,  but  are  grouped  together  at  the 
begirming  of  the  mam  source  code  file  and  are  easily  modified  by  anyone  &miliar  with  general  programming 
practices.  The  geographic  and  force  structure  databases  were  built  in  one  file  using  a  spreadsheet.  This  file  is 
easily  modified  by  a  user  ^miliar  with  the  operation  of  spreadsheets,  and  requires  no  moHifioatinn  of  the 
computer  code. 
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E.  OUTPUT  FILES 


1.  Unit  movement  (Motion.dat) 

This  file  is  updated  during  every  time  step  and  records  the  Cartesian  location  of  each  unit  along 
with  its  unit  identification  number.  This  file  is  used  to  verify  the  units  are  moving  correctly  in  accordance 
with  both  the  initial  database  parameters  and  movement  algorithms. 

2.  Unit  losses  (Loss.dat) 

This  file  is  updated  each  time  a  imit  is  killed.  The  simulation  records  the  unit  identification 
number  and  the  simulation  time  at  which  it  was  killed. 

3.  Unit  status  (Unitdat) 

This  file  is  the  primary  data  source  for  measures  of  effectiveness  (MOE)  calculations  and 
evaluation  of  the  battle.  This  file  contains  a  complete  listing  of  the  characteristics  of  each  unit  at  the 
begiiming  of  the  simulation  (t  =  0),  and  the  completion  of  the  simulation.  The  information  includes  the 
number  of  hits  suffered  by  each  unit,  and  number  of  weapons  left  over  after  the  battle.  The  first  list  is  used 
in  verifying  that  the  simulation  correctly  input  the  data  fi-om  the  initial  database.  The  second  list  provides 
information  for  measure  of  effectiveness  (MOE)  calculations.  At  the  end  of  the  simulation,  the  program 
also  outputs  to  this  file  the  variables  that  keep  track  of  the  number  missiles  fired,  attack  and  fighter 
aircraft  launched,  luiits  killed,  and  aircraft  lost  on  strike  or  interdiction  missinns 
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V.  EXPERIMENTATION  AND  ANALYSIS 


A  OVERVIEW 

The  output  of  this  simulation  is  a  function  of  numerous  input  parameters.  The  goal  of  this  thesis  is  to 
construct  and  exercise  a  proposed  simulation  model  that  has  the  cap^ility  to  examine  the  effect  of  varying 
conditions  that  relate  to  the  sensors  and  we^ns  available  to  the  units  that  comprise  a  United  States  carrier 
battle  group.  A  fectorial  experiment  was  used  to  stu<fy  one  instance  of  the  conflirt  described  in  this  thesis  to 
examine  how  changing  three  ofthe  initial  conditions  affect  the  outcome  of  the  battle.  ThQ  Minitab®  Release  11 

for  Windows  statistical  software  package,  whose  capabilities  are  described  in  the  Minitab  ®  Reference  Manual 

Release  11  for  Windows  (Minitab,  1996),  is  used  to  demonstrate  techniques  for  analyzing  the  ouqjut  from  this 
simulation 

The  result  of  the  conflict  is  sununarized  using  seven  measures  of  effectiveness  (MOE),  which  judge  the 
outcome  of  simulation  trials  with  respect  to  the  success  of  the  United  States  and  mnltinatinnal  forces.  Output  is 
gathered  by  performing  computer  trials  of  the  simulation;  this  information  is  used  to  calculate  the  seven  MOE. 
A  2^  fifll  &ctorial  erqreriment  using  ten  replications  of  the  computer  simulation  for  eadi  &ctor  level  was 
conducted  and  analysis  performed  on  the  output  (specific  parameters  of  the  initial  data  base  and  program  are 
modified  as  required  to  vary  the  conditions  to  conduct  the  experiment). 

B.  MEASURES  OF  EFFECTIVENESS 

The  measures  of  effectiveness  are  defined  to  represent  the  ability  of  the  U.  S.  force  to  accon^Ush  its 
mission  while  minimizing  damage  to  its  own  forces,  and  to  describe  the  results  in  terms  of  the  United  States 
combat  effectiveness  and  survivability.  The  measures  of  effectiveness  do  not  measure  the  attrition  rate  of 
attacking  weapons  since  they  depend  only  on  the  probability  of  defeat  (p^jg^)  by  a  target’s  defensive  mechar^ms 
in  the  computer  implemontation  of  the  simulation  The  following  MOE  are  defined  for  the  model: 

1.  Total  number  of  hits  sustained  among  all  ships  comprising  the  carrier  battle  group  over  the  two-d^ 

battle. 
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4.  Number  of  multinational  force  ships  remaining  at  the  conclusion  of  the  two-day  battle. 

5.  Combined  number  of  U.  S.  missiles  and  attack  aircraft  launched  against  Iranian  targets  over  the 
two-day  battle. 

6.  Total  Niunber  of  hits  scored  against  all  Iranian  land-based  targets  over  the  two-day  battle. 

7.  Total  number  of  Iranian  attack  aircraft  launched  against  U.  S.  sh^  over  the  two-d^  battle. 

C.  EXPERIMENTATION 

The  2^  factorial  experiment  varied  three  conditions  which  are  the:  (1)  The  availabihty  of  sateUite 
intelligence  to  the  U.  S.  forces;  (2)  U.  S.  sensor  accuracy;  and  (3)  rtuniber  of  U.  S.  aircraft  carriers.  The  United 
States  either  receives  satellite  information  from  one  pair  of  satellites,  or  it  receives  no  satellite  information  the 
entire  simulation,  ^pon  detection,  a  sensor  reports  the  perceived  location  of  the  target  which  is  the  actual 
location  offset  by  a  normal  error.  For  one  level  of  the  sensor  accuracy  fector,  the  normal  errors  for  all  sensors 
have  mean  zero  and  variance  of  16  miles  for  surfoce  and  ESM  sensors  and  .25  mile  for  air  sensors.  For  the  other 
level  of  sensor  accuracy,  the  normal  errors  of  the  United  States  and  multinational  force  sensors  errors  have  mean 
zero  and  variance  zero  and  so  the  target  is  located  perfectly,  Iranian  sensor  errors  have  mean  zero  and  variance 
of  16  miles  for  surfoce  and  ESM  sensors  and  .25  mile  for  air  sensors.  Hence,  the  sensor  accuracy  foctor 
eliminates  the  sensor  error  and  causes  all  U.  S.  and  midtinational  force  sensors  to  provide  perfect  information  on 
the  location  of  a  target.  The  U.  S.  operates  with  one  or  two  aircraft  carriers  in  the  battle  group.  The  addition  of  a 
second  aircraft  carrier  provides  a  second  air  wing  with  its  attadc,  fighter,  and  reconnaissance  aircraft.  The 
design  of  the  ejqperiment  and  the  analysis  techniques  that  are  applied  to  devdop  conclusions  are  described  in 
Statistics  for  Experimenters,  An  Introduction  to  Design,  Data  Analysis,  and  Model  Building  (Box,  1978, 
pp.  306  -  321).  The  experimental  runs  are  arranged  m  the  standard  order  as  shown  in  Tables  1  and  2. 
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Table  1.  Factor  Levels  for  2^  Factorial  E^eriment 


Variable 

- 

+ 

1.  Satdlites; 

0 

2 

2.  Sensor  accuracy: 

Normal  errors 

Perfect  information 

3.  Aircraft  Carriers : 

1 

2 

Table  2.  2^  Factorial  Design,  Trials  Arranged  in  the  Standard  Order 

Trial  fector  in  simulation 

Number  Satellites  Sensor  accuracy _ Carriers 


1 

2 

3 

4 

5 

6 

7 

8 


0  Normal  errors 

2  Normal  eirors 

0  Perfect  information 

2  Perfect  information 

0  Normal  errors 

2  Normal  errors 

0  Perfect  infoimadon 

2  Perfect  information 


D.  ANALYSIS 

1.  Overview 


1 

1 

1 

1 

2 

2 

2 

2 


The  full  "t  &ctorial  experiment  was  performed  with  ten  replications  performed  for  each  factor  level. 
Each  of  the  MOE  computed  from  the  output  of  the  simulation  runs  is  examined  using  statistical  analysis  tools. 
The  experiment  results  are  included  as  Appendix  F. 

2.  Analysis  Techniques 

The  first  step  in  exploring  how  the  initial  conditions  affected  the  outcome  was  to  graphically  explore 
the  results  of  the  e?q}erimentation.  The  individual  data  points  were  plotted  and  both  these  gr^hs  and  the  data 
are  included  as  Appendix  F.  The  plots  were  examined  to  note  significant  trends  and  determine  what  analysis 
techniques  to  ^ly  for  each  MOE. 

If  the  plot  of  the  observations  of  an  MOE  leads  one  to  believe  that  the  observational  distribution  is 
normal,  a  fractional  factorial  anafysis  fit  was  computed  for  the  MOE  using  the  Minitab^  statistical  software 
package.  This  analysis  first  computes  a  fractional  fectorial  estimate  of  how  each  fector  individually  (main 
effect),  and  when  combined  with  other  factor(s)  (multiple-factor  interaction  effect)  affects  the  MOE.  The  result 
of  this  analysis  provides  an  estimate  of  how  each  effect  changes  the  MOE  and  a  probability  value  (p-value). 
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which  is  the  piobabilily  of  getting  a  test  statistic  value  at  least  as  extreme  as  flie  one  actually  obtained  under  the 
null  hypothesis  that  there  is  no  effect.  If  the  p-value  is  less  than  a  desired  significance  level  a,  the  effect  be 

significant  and  requires  interpretatioa  The  main  effect  of  a  variable  diould  be  individualfy  interpreted  onty  if 
there  is  no  evidence  that  the  variable  interacts  with  other  variables.  When  there  is  evidence  of  one  or  more  such 
interaction  effects,  tte  interacting  variables  should  be  considered  jointly,  and  other  anatysis  techniques  used  to 
further  explore  the  main  effects.  In  using  this  analysis  technique,  first  look  at  the  three-feaor  interaction  and 
decide  if  it  is  significant  Ifthethree-toctor  interaction  is  significant,  interpretation  (ff  main  effects  and  two- 
fector  interaction  effects  are  not  valid  on  the  basis  of  this  anafysis  tedmique  alone.  If  there  is  no  three-factor 
interaction,  but  one  or  more  two-factor  interactimis,  only  the  two-fector  interactions  are  interpreted  and  no 
conclusions  can  be  drawn  about  the  main  effects.  ^Box,  1978,  r).  3 17-3 18). 

The  software  performs  a  three-fector  analysis  of  variance  (ANO  VA),  which  is  a  technique  that 
determines  whether  the  different  effect  estimates  differ  firom  zero  more  than  could  reasonably  be  expected  under 
the  null  Itypothesis  of  no  effect.  The  ANOVA  technique  tests  the  null  hypothesis  that  the  means  of  the  effects  for 
each  level  of  interaction  are  equal  to  zero,  against  the  alternative  hypothesis  that  there  is  at  least  one  that  is  not 
zero.  The  software  corrqrutes  an  F-ratio  and  a  p-value  for  each  level  of  interaction  (main,  two-way,  three-way). 
This  p-vahie  is  the  probabihty  of  obtaining  a  realization  from  the  F-distnfrution  under  the  null  hypothesis  of  no 
effect  at  least  as  large  as  the  observed  F-ratio.  The  theory  behind  the  computation  and  use  of  the  F-ratio  and  p- 
value  in  the  ANOVA  analysis  can  be  found  in  Probability  and  Statistics  for  Engineering  and  the  Sciences 
(Devore,  1995,  pp.  444-452). 

A  second  procedure,  the  Kruskal-Wallis  test,  which  tests  whether  an  arbitrary  number  of  X:  independent 
populations  are  identical  is  also  used,  and  is  applicable  to  all  the  MOE  results  regardless  of  their  distributioa 
This  test  is  similar  in  spirit  to  one-w^  ANOVA  analysis  and  tests  the  null  hypothesis  that  the  k  pqjulations  are 
identical,  against  the  alternative  hypothesis  that  there  are  at  least  two  of  the  peculations  that  are  Hifferfint  The 
Kruskal-Wallis  test  computes  an  /f-value  and  compares  it  to  a  distribution  to  compute  a  p-value.  If  the 

p-value  is  so  small  that  the  null  hypothesis  is  not  supported,  it  is  concluded  that  the  compared  populations 
are  not  all  the  same.  A  theoretical  development  and  explanation  of  this  test  is  found  in  Nonparametric 
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Methods  for  ^mtitcUive  Amlysis  {GMam,  1985,  h).  173-179).  TheMratoZ»®Kniskal-Wallis  test  that  is 
adjusted  for  ties  was  used  to  obtain  an  //-value  and  a  p-vahie  for  each  &ctoi. 

The  A///3/7ah®  statistical  software  package  was  then  used  to  develop  cube  plots  of  the  means  for  each 
MOE  and  conq)ute  basic  summary  statistics  (mean,  standard  error  of  the  mean,  and  standard  deviation)  of  each 
individual  trial  for  all  MOE. 

3.  MOE  1 

Since  the  data  points  for  MOE  1  appear  normally  distributed,  as  shown  in  the  plot  included  in 
Appendix  F,  a  factorial  analysis  was  computed  and  is  shown  in  Table  3.  For  any  effect  with  a  p-value  less 
than  .05  (a  =  .05  significance  level),  the  null  hypothesis  that  the  mean  of  the  effect  is  equal  to  zero  is 
rejected  and  the  alternative  hypothesis  that  the  mean  is  different  from  zero  is  accepted;  hence,  the  effect  is 
considered  significant.  The  use  of  this  significance  level  is  consistent  throughout  all  analyses  that  are 
discussed  in  this  thesis.  The  ANOVA  table  indicates  that  there  is  at  least  one  strong  two-way  interaction, 
with  a  conesponding  p-value  of  .03.  Due  to  the  strong  two-way  interaction  which  occurs  between  the  use 
Table  3.  Analysis  of  Variance  Table  for  MOE  1  with  Estimated  Effects 


Source 

df 

SS 

MS 

F-ratio 

p-value 

Main  effects 

3 

5145.40 

1715.10 

12.37 

0.00 

Two-way  interactions 

3 

1343.30 

447.80 

3.23 

0.03 

Three-way  interaction 

1 

515.10 

515.10 

3.71 

0.06 

Residual  error 

72 

9986.50 

138.70 

Total 

79 

16990.40 

Factor 

Effect 

p-value 

mean 

28.30 

- 

satellite 

-0.43 

0.87 

sensor 

-10.38 

0.00 

carrier 

12.23 

0.00 

satellite  &  sensor 

3.83 

0.15 

satellite  &  carrier 

5.93 

0.03 

sensor  &  carrier 

4.18 

0.12 

all 

5.08 

0.06 

Note:  The  value  for  the  mean  effect  is  the  overall  mean  value  of  MOE  1  for  the  eighty 
simulation  runs.  The  values  for  the  other  effects  represent  the  estimated  change  in  MOE  1  when  an 
individual  factor  or  combination  of  factors  are  at  their  high  level  (+)  as  compared  to  the  when  aU  of  the 
factors  are  at  their  low  levels  (-).The  p-value  of  the  estimated  effect  is  the  probability  that  the  value  for  the 
effect  is  actually  zero. 
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of  satellite  intelligence  and  an  additional  carrier  (p-value  =  .03),  the  large  main  effects  of  perfect  sensor 
information  or  adding  an  additional  carrier  cannot  be  interpreted  based  on  this  analysis  technique.  Before 
accepting  the  results  of  the  factorial  analysis,  a  normal  probability  plot  of  the  residuals  is  performed  to 
verify  the  normality  assumption  used  to  enter  into  this  analysis.  The  plot  of  the  residuals  in  Figure  7, 
shows  that  the  normality  assumption  is  probably  reasonable,  and  the  results  of  this  analysis  technique  are 
valid. 


Normal  Probability  Plot  of  the  Residuals 

(response  Is  moe  1) 


Figure  7.  Normal  Probability  Plot  of  the  Residuals  for  the  Fitted  Factorial 
Model  of  MOE  1 

The  most  information  is  obtained  about  this  MOE  by  examining  a  cube  plot  of  the  means,  which 
is  shown  in  Figure  8.  A  comparison  of  the  means  for  the  two  levels  of  sensor  accuracy  clearly  shows  that 
when  perfect  sensor  information  is  used,  the  mean  number  of  hits  sustained  decreases  in  all  cases  except 
when  an  additional  carrier  is  added  to  the  battle  group.  If  the  means  of  the  two  levels  for  the  number  of 
carriers  in  theater  are  compared,  it  is  noted  that  the  mean  number  of  hits  increases  with  the  addition  of  a 
second  carrier,  regardless  of  the  changes  in  the  other  two  factors  in  all  but  one  case  The  exception  is 
when  there  are  two  carriers,  perfect  sensor  information,  and  no  satellite  intelligence.  No  reasonable 
ejqplanation  is  offered  for  this  exception,  and  it  reqmres  further  study.  The  only  conclusion  that  can  be 
drawn  is  that  using  perfect  sensor  information  when  operating  with  one  carrier  causes  the  mean  number 
of  hits  sustained  by  the  ships  among  the  carrier  battle  group  to  decrease. 
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-  Means  for  MOE  1 


Cube  Plot 


24.1  38.S 


satell 


satellites 

sensor  accuracy 

aircraft  carriers 

mean 

se  mean 

std  dev 

0 

normal  errors 

1 

32 

3.5 

11.06 

2 

normal  errors 

1 

26.9 

4.19 

13.25 

0 

perfect  information 

1 

18.7 

4.37 

13.82 

2 

perfect  information 

1 

11.1 

4.33 

13.68 

0 

normal  errors 

2 

39.2 

4.29 

13.56 

2 

normal  errors 

2 

35.8 

2.62 

8.28 

0 

perfect  information 

2 

24.1 

3.61 

11.42 

2 

perfect  information 

2 

38.5 

2.26 

7.14 

How  to  Read:  T  he  main  effect  of  each  factor  is  seen  to  be  a  difference  between  two  averages, 
half  of  the  eight  results  being  included  in  one  average  and  half  in  the  other.  The  main  effects  may  be 
viewed  as  a  contrast  between  observations  on  parallel  faces  of  the  cube  plot.  Similarly,  each  two-factor 
interaction  is  viewed  as  a  contrast  bet^veen  results  on  two  diagonal  planes.  The  mean  values  displayed  on 
the  cube  are  also  listed  in  tabular  form  along  with  their  standard  error  of  the  mean  and  standard  deviation. 

Figure  8.  Means  of  the  Total  Number  of  Hits  Sustained  Among  All  Ships  Comprising  the  CVBG 
Over  the  Two-Day  Battle  and  Summary  Statistics. 


4.  MOE  2 

The  number  of  hits  sustained  by  the  ships  comprising  the  multinational  force  varies  between  two 
numbers  and  does  not  appear  normally  distributed,  so  an  ANOVA  analysis  was  not  attempted.  Each  ship 
comprising  this  group  has  a  maximum  number  of  hits  it  can  sustain  before  it  is  considered  killed.  The 
combined  total  hits  for  the  four  ships  is  36,  and  in  27  out  of  80  simulation  runs  all  ships  are  killed.  Since 
a  sortie  may  contain  numerous  weapons,  it  may  score  more  hits  than  needed  to  destroy  a  ship.  This 
means  the  number  of  hits  sustained  may  be  greater  than  the  maximum  required  to  kill  a  ship.  In  two 
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cases,  the  mean  number  of  ships  killed  is  greater  than  36.  When  satellite  intelligence  and  an  additional 
carrier  is  in  theater,  the  mean  number  of  hits  sustained  by  the  multinational  force  is  estimated  as  20  from 
the  data  listed  in  Figure  9.  This  suggests  that  some  of  the  ships  are  not  killed. 


A  Kruskal-Wallis  test  is  performed  for  each  individual  factor.  Each  test  compared  the  40 
observations  of  MOE  2  for  a  factor  at  its  initial  level  (-)  with  the  40  observations  at  its  optimal  (+)  level. 
The  results  of  this  analysis  are  summarized  in  Table  4.  Since  the  p-values  for  using  satellite  intelligence 
and  an  additional  carrier  are  less  than  0.05  they  are  considered  significant.  This  suggests  that  using 
satellite  intelligence  or  an  additional  carrier  causes  the  number  of  hits  sustained  by  the  multinational  force 
to  decrease. 


Table  4.  Kruskal-Wallis  Analysis  for  MOE  2 


factor 

H-statistic 

p-value 

satellite 

15.52 

0.00 

sensor 

0.19 

0.66 

carrier 

13.12 

0.00 

Comparing  the  mean  number  of  hits  sustained  for  each  case,  it  is  determined  that  only  when  both  an 
additional  carrier  and  satellite  intelligence  is  available  does  the  mean  number  of  hits  sustained  by  the 
multinational  force  decrease,  regardless  of  sensor  accuracy.  Since  the  ANOVA  table  indicates  that  a  two- 
way  interaction  appears  to  be  present,  the  results  of  the  Kruskal-Wallis  test  should  be  viewed  with  caution. 


Cube  Plot  -  Means  for  MOE  2 
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Figure  9.  Means  of  the  Total  Number  of  Hits  Sustained  Among  All  Ships  Comprising  the 
Multinational  Force  over  the  Two-Day  Battle  and  Summary  Statistics. 

The  only  conclusion  that  can  be  stated  is  that  the  mean  number  of  hits  decrease  when  the  U.  S. 
forces  utilize  both  satellite  intelligence  and  an  additional  carrier. 

s.  MOE3 

Only  a  few  ships  among  ttie  seven  comprising  the  carrier  battle  group  are  defeated  during  the 
two-day  battle,  and  the  most  ships  lost  during  any  of  the  80  Riniiilatinn  runs  is  two.  This  data  is  not 
normal  and  only  the  Kniskal-Wallis  test  is  appropriate  and  the  results  are  shown  in  Table  5.  The  results 
of  the  Kruskal-Wallis  tests  indicate  that  the  use  of  perfect  sensor  information  is  significant.  Examining 
Figure  10  shows  that  the  mean  fiaction  of  ships  that  survive  increases  when  perfect  sensor  information  is 
used  in  all  cases  except  when  there  are  two  carriers.  These  data  points  are  so  clustered,  that  any  difference 
in  the  means  may  be  due  to  normal  variations  as  indicated  by  the  summary  statistics.  No  conclusions  can 
be  drawn  for  this  MOE  from  this  analysis. 


Cube  Plot  -  Means  for  MOE  3 
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Figure  10.  Means  of  the  Fraction  of  CVBG  Ships  Remaining  at  the  Conclusion  of  the  Two-Day 
Battle  and  Summary  Statistics. 

Table  5.  Kruskal-Wallis  Analysis  for  MO£  3 


factor  H-statistic  p-value 

satellite  2.45  0.12 

sensor  5.89  0.02 

carrier  2.02  0.16 


6.  MOE4 

The  number  of  multinational  force  ships  that  stirvive  the  battle  is  directly  related  to  the  number 
of  hits  that  the  ships  sustain,  and  there  is  a  strong  correlation  between  MOE  2  and  MOE  4.  Again 
ANOVA  analysis  is  not  appropriate,  and  Kruskal-Wallis  tests  are  performed.  The  tests,  which  are 
summarized  in  Table  6,  indicate  that  both  the  use  of  satellite  intelligence  and  an  additional  aircraft  are 
significant.  The  cube  plot  of  the  means  clearly  shows  similar  findings  to  those  for  MOE  2;  when  satellite 
intelligence  is  available  along  with  a  second  carrier  the  mean  number  of  multinational  force  ships  that 

Table  6.  Kruskal-Wallis  Analysis  for  MOE  4 


factor 

H-statistic 

p-value 

satellite 

30.55 

0.00 

sensor 

0.39 

0.53 

carrier 

13.49 

0.00 

survive  the  two-day  battle  increases.  This  is  consistent  with  MOE  2  where  the  use  of  satellite  intelligence 
with  two  carriers  caused  the  mean  number  of  hits  sustained  by  the  multinational  force  to  decrease.  Since  a 
two-way  interaction  appears  to  be  present,  the  results  of  the  Kmskal-Wallis  test  should  be  viewed  with 
caution.  The  only  conclusion  that  can  be  stated  is  that  the  mean  number  of  multinational  force  ships  that 
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survive  the  two-day  battles  increases  when  the  U.  S.  forces  utilize  both  satellite  intelligence  and  an 
additional  carrier. 


Cube  Plot  -  Means  for  MOE  4 
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Figure  11.  Means  of  the  Number  of  Multinational  Force  Ships  Remaining  at  the  Conclusion  of  the 


Two-Day  Battle  and  Summary  Statistics. 

7.  MOE  5 

The  data  points  for  MOE  5  appear  normally  distributed  as  shown  in  the  plot  included  in 
Appendix  F.  A  factorial  fit  was  computed  and  a  plot  of  the  residuals  was  generated  to  check  the  validity 
of  the  fit.  The  residual  plot,  which  is  displayed  in  Figure  12,  refutes  the  normality  assumption  and  the 
factorial  fit  was  discarded.  The  Kruskal-Wallis  tests,  which  are  summarized  in  Table  7,  once  again 
indicate  that 


Table  7.  Kruskal-Wallis  Analysis  for  MOE  5 


factor 

H-statistic 

p-value 

satellite 

9.25 

0.00 

sensor 

0.01 

0.91 

carrier 

44.88 

0.00 
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Normal  Probability  Plot  of  the  Residuals 

(response  is  moe  5) 


Figure  12.  Normal  Probability  Plot  of  the  Residuals  for  the  Fitted  Factorial 
Model  of  MOE  5 

both  the  addition  of  satellite  intelligence  and  a  second  carrier  are  significant.  It  can  be  seen  in  Figure  13 
that  when  a  second  carrier  is  added,  the  mean  number  of  U.  S.  missiles  and  attack  aircraft  launched 
against  Iranian  targets  significantly  increases.  The  mean  number  of  missiles  and  attack  aircraft  launched 
also  increases  whenever  satellite  information  is  available,  but  in  the  case  when  only  one  carrier  is  in  the 
battle  group;  this  perceived  increase  is  within  the  bounds  of  the  standard  errors  of  the  means.  This  leads 
one  to  believe  that  the  second  carrier,  which  brings  an  additional  air  wing  to  both  attack  the  Iranian 
targets,  and  provides  a  second  set  of  reconnaissance  aircraft,  is  the  reason  the  mean  number  of  missile  and 
attack  aircraft  launches  increases. 
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Figure  13.  Means  of  the  Combined  Number  of  U.  S.  Missiles  and  Attack  Aircraft  Launched  Against 
Iranian  Targets  over  the  Two-Day  Battle  and  Summai^  Statistics. 


8.  MOE6 

Since  the  data  points  for  MOE  6  appear  normally  distributed,  as  shown  in  the  plot  included  in 
Appendix  F,  an  ANOVA  table  was  computed.  The  results  of  the  ANOVA  appear  in  Table  8  and  indicate 
that  there  is  at  least  one  strong  two-way  interaction,  with  a  p-value  of  .00.  Due  to  the  strong  two-way 
interaction  between  an  additional  carrier  and  satellite  intelligence  (p-value  =00),  the  two  must  be 
considered  jointly  and  no  conclusions  about  the  large  main  effects  of  the  availability  of  satellite 
intelligence  or  adding  an  additional  carrier  can  be  stated  based  on  this  analysis  technique.  A  normal 
probability  plot  of  the  residuals  shown  in  Figure  14  supports  the  normality  assumption  which  is  required 
for  this  analysis. 
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Normal  Probability  Plot  of  the  Residuals 

(response  is  moe  6) 


I 

I 


Figure  14.  Normal  Probability  Plot  of  the  Residuals  for  the  Fitted  Factorial 
Model  of  MOE  6 


The  estimate  of  the  two-way  interaction  of  satellites  and  carriers  suggests  that  the  combined  effect  of  using 
satellite  intelligence  and  an  additional  carrier  suggests  causes  the  mean  number  of  U.  S.  missile  and  attack 
aircraft  launches  to  increase  by  about  7. 

Table  8.  Analysis  of  Variance  Table  for  MOE  6  with  Estimated  Effects 


Source 

df 

SS 

MS 

F-ratio 

p-value 

Main  effects 

3 

8120.90 

2706.98 

63.55 

0.00 

Two-way  interactions 

3 

2167.20 

722.41 

16.96 

0.00 

Three-way  interaction 

1 

56.11 

56.11 

1.32 

0.26 

Residual  error 

72 

3067.10 

42.60 

Total 

79 

13411.40 

Factor 

Effect 

p-vaiue 

mean 

104.79 

- 

satellite 

17.13 

0.00 

sensor 

2.03 

0.17 

carrier 

10.43 

0.00 

satellite  &  sensor 

0.58 

0.70 

satellite  &  carrier 

10.38 

0.00 

sensor  &  carrier 

-0.63 

0.67 

all 

-1.68 

0.26 

Note:  The  value  for  the  mean  effect  is  the  overall  mean  value  of  MOE  6  for  the  eighty 
simulation  runs.  The  values  for  the  other  effects  represent  the  estimated  change  in  MOE  6  when  an 
individual  factor  or  combination  of  factors  are  at  their  high  level  (+)  as  compared  to  the  when  all  of  the 
factors  are  at  their  low  levels  (-).The  p-value  of  the  estimated  effect  is  the  probability  that  the  value  for  the 
effect  is  actually  zero. 
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Further  examination  of  the  cube  plot  of  the  means  in  Figure  13  shows  that  the  same  eiBfect  is  occurring  as 
in  MOE  5.  The  addition  of  a  carrier  causes  the  mean  number  of  hits  scored  against  Iranian  targets  to 
increase  as  more  attack  aircraft  are  launched  due  to  the  additional  airwing  located  on  the  second  carrier. 
The  increase  when  there  is  one  carrier  and  satellite  intelligence  available  falls  within  the  standard 
deviations  of  the  compared  means.  A  definitive  statement  about  the  one  way  effect  of  satellite  intelligence 
is  questionable  since  it  may  be  accounted  for  by  the  standard  errors  of  the  means. 

Cube  Plot  -  Means  for  MOE  6 
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Figure  15,  Means  of  the  Combined  Number  of  U.  S.  Missiles  and  Attack  Aircraft  Launched  Against 
Iranian  Targets  over  the  Two-Day  Battle  and  Summary  Statistics 

9.  MOE  7 

Since  the  data  points  for  MOE  7  appear  normally  distributed,  as  shown  in  the  plot  in  Appendix 
F,  a  fractional  factorial  fit  was  computed.  A  residual  plot,  which  is  shown  in  Figure  16,  verifies  that  the 


assumption  of  normality  is  reasonable.  The  results  of  the  factorial  analysis  performed  for  MOE  7  are 
summarized  in  Table  9.  Since  there  are  no  significant  two-way  and  three-way  interactions,  the  use  of 
satellite  intelligence  (p-value  =  0.00)  and  perfect  sensor  information  (p-value  =  0.00)  are  identified  as 
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significant  main  effects.  The  estimate  of  the  main  effect  of  satellites  indicates  that  use  of  satellite 
intelligence  causes  the  mean  number  of  attack  aircraft  launched  against  U.  S.  and  multinational  force 
ships  to  decrease  by  about  36.  The  estimate  of  the  main  effect  of  sensor  accuracy  indicates  that  when 
perfect  sensor  information  is  used  the  mean  number  of  attack  aircraft  launched  decrease  by  about  19. 


Normal  Probability  Plot  of  the  Residuals 

(response  is  moe  7) 
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Figure  16.  Normal  Probability  Plot  of  the  Residuals  for  the  Fitted  Factorial  Model  of  MOE  7 


Table  9.  Analysis  of  Variance  Table  for  MOE  7  with  Estimated  Effects 


Source 

df 

ss 

MS 

F-ratio 

p-value 

Main  effects 

3 

32790.00 

10930.00 

20.17 

0.00 

Two-way  interactions 

3 

4370.60 

1456.90 

2.69 

0.05 

Three-way  interaction 

1 

583.20 

583.20 

1.08 

0.30 

Residual  error 

72 

39008.00 

541.80 

Total 

79 

76751.80 

Factor 

Effect 

p-value 

mean 

150.05 

- 

satellite 

-35.90 

0.00 

sensor 

-18.70 

0.00 

carrier 

-1.00 

0.85 

satellite  &  sensor 

9.10 

0.09 

satellite  &  earner 

-9.60 

0.07 

sensor  &  carrier 

6.60 

0.21 

all 

5.40 

0.30 

The  Kruskal-Wallis  tests  and  an  examination  of  the  cube  plot  in  Figure  17  support  these  conclusions. 


Table  10.  Kruskal-Wallis  Analysis  for  MOE  7 


factor 

H-statistic 

p-value 

satellite 

27.62 

0.00 

sensor 

6.40 

0.01 

carrier 

.07 

0.79 

Cube  Plot  -  Means  for  MOE  7 
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Figure  17.  Means  of  the  Total  Number  of  Iranian  Attack  Aircraft  Launched  Against  U.  S.  ships 
over  the  Two-Day  Battle  and  Summary  Statistics 


E.  EXPERIMENT  CONCLUSIONS 

The  track  which  the  ships  in  the  U.  S.  battle  group  follow  provides  protection  keeping  them  outside 
of  the  maximum  weapon  range  of  most  Iranian  weapons  for  significant  portions  of  time  over  the  two-day  battle. 
Hence,  MOE  1  and  MOE  3  which  measure  the  survivability  of  the  ships  among  the  CVBG  do  not  show  large 
reqxtnses  to  changes  in  each  of  the  three  conditions.  The  mean  number  of  hits  against  these  ships  does  decrease 
when  perfect  sensor  information  is  used  for  a  CVBG  consisting  of  one  carrier.  The  only  other  MOE  which  is 
affected  by  using  perfect  sensor  information  is  MOE  7  which  measures  the  number  of  attack  aircraft  sorties  that 
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Iran  launches  at  U.  S.  and  multinatiQnal  force  ships.  It  appears  that  the  use  of  perfect  infonnation  the  U.  S. 
increases  the  probability  that  a  United  States  missile  or  aircraft  strike  scores  a  hit  against  an  Iranian  target, 
therelty  reducing  number  of  wejqwns  that  Iran  can  fire  at  the  U.  S.  ships.  It  is  concluded  that  the  use  of  pafect 
sensor  information  increases  the  probability  that  ships  among  the  carrier  battle  group  survive  over  the  two-day 
battle. 

The  trade  of  the  multinational  force  places  these  ships  within  Iranian  weapons  range  early  in  the  battle. 
Since  reconnaissance  aircraft  are  only  flown  from  a  carrier,  the  ability  to  locate  and  target  Iranian  SSM 
launchers  and  air  bases  which  are  located  near  the  multinational  force  trade  does  not  occur  until  the  CVBG 
nears  the  Straits  of  Hormuz.  As  a  result,  the  specific  scenario  used  in  the  simulation  makes  the  Hftgtnicfion  of  all 
multinational  force  ships  by  Iran  likely.  The  survivability  of  the  multinational  force  is  measured  hy  MOE  2  and 
MOE  4.  These  measures  show  that  flie  only  condition  that  enables  the  some  of  the  multinational  force  shipg  to 
successfully  complete  their  transit  occur  when  satellite  intelligence  is  available  and  a  second  aircraft  carrier  is 
among  the  battle  group.  The  availability  of  satellite  intelligence  provides  the  U.  S.  forces  a  method  of  detecting 
Iranian  targets  that  are  either  not  detected  Ity  shipboard  sensors  and  reconnaissance  aircraft,  or  are  detected  after 
they  have  launched  a  significant  amount  of  weapons  against  the  multinational  force.  Not  onfy  must  the  U.  S. 
detect  these  targets,  but  additional  strike  assets  in  the  form  of  the  second  carrier  air  wing  must  be  available  to 
attack  these  targets.  Hence,  for  the  scenario  of  the  simulation,  it  is  concluded  that  for  the  CVBG  to  successfully 
acconplish  its  mission,  which  is  to  ensure  the  multinational  force  reaches  the  Arabian  Sea,  satellite  intelligence 
and  a  second  aircraft  carrier  are  required. 

The  mean  number  of  U.  S.  missiles  and  attack  aircraft  launched  at  Iranian  targets  (MOE  5)  and  the 
mean  number  of  hits  scored  against  Iranian  targets  (MOE  6)  both  inaease  when  an  aHditinnal  carrier  is  artHwi  to 
the  battle  group.  This  is  not  a  surprising  conclusion.  It  ^jpears  that  the  availability  of  satellite  intelligence 
have  also  caused  these  MOE  to  increase,  but  the  analysis  shows  that  the  change  due  to  this  availability  may  be 
due  to  normal  variatioa 
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VL  CONCLUSIONS  AND  FURTHER  STUDY 


A.  CONCLUSIONS 

The  analysis  of  the  results  obtained  from  the  frctorial  e)q}eiiment  demonstrate  that  the  outcome  of  this 
simulation  model  is  reqwnsive  to  the  changes  related  to  sensors  and  weapons  which  were  varied  during  the 
ejqjeiiment  The  simulation  results  for  this  one  specific  conflict  show  that  the  use  of  intelligence  infotmation 
obtained  from  satellite  assets  combined  with  the  assistance  of  a  second  aircraft  carrier  enables  the  CVBG  to 
accomplish  its  missioa  When  these  conditions  are  not  met,  all  multinational  force  ships  are  usually  destroyed. 
The  use  of  sensors  which  provide  ground  truth  position  reports  for  all  detections  does  not  significantly  improve 
the  mission  success  of  the  CVBG  which  implies  that  the  simulation  is  not  as  saisitive  to  changes  in  sensor 
accuracy  as  desired. 

The  simulation  model  approach  presented  shows  promise  of  having  infinmative  predictive  raphilittpg 
but  it  present^  does  not  appear  to  possess  sufficient  detail  to  properly  measure  all  aspects  of  CiSR  A  more 
detailed  stu^  of  this  model  is  needed  to  prqjerly  evaluate  its  present  capability  to  assess  the  impact  of  cbangfli!  m 
C^'lSR  and  inqrrove  on  its  shortcomings.  After  conducting  this  initial  erq)eiiment  and  anatysis,  it  is  dear  that  for 
this  simulation  model  to  be  a  useful  predictive  tool  to  a  dedsion  maker,  it  needs  to  incorporate  more  realistic  data 
fusion  and  decision  models  (rules)  under  which  the  battle  is  to  be  conducted. 

B.  FURTHER  STUDY 

This  is  an  initial  approach  to  a  stochastic  simulation  model  which  nti1i7e.R  very  simplistir  decision 
algorithms.  The  first  step  in  irtqrroving  fliis  simulation  model  so  that  it  is  a  more  powerful  predictive  tool  is  to 
refine  the  decision  modules  that  address  targeting,  detection,  conflict  adjudication,  and  sensor  employment  The 
algorithms  developed  do  not  follow  any  formal  U.  S.  or  Iranian  doctrine.  This  model  can  be  aHaptpH  to 
incorporate  more  detailed  decision  algorithms  using  U.  S.  naval  and  Iranian  tactics  and  dnrtrinp 

Another  area  where  further  study  is  recommended  is  in  the  detection  model.  The  sensors  are  given 
generic  probabihties  of  detection  which  are  not  conditioned  on  a  specific  type  of  target.  Also  the  inflnpnrp  of 
weather,  jamming,  and  decoys  are  not  addressed  and  when  incorporated  may  improve  the  realism  of  the  model. 
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The  ejqjerimentation  performed  using  this  simulation  model  was  veiy  limited  in  scx)pe.  The  response 
of  this  simulation  can  be  further  studied  by  changing  other  parameters,  such  as  detection  probabilities,  sensor 
detection  ranges,  we^ns  loadout,  or  weapons  range. 
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APPENDIX  A.  GEOGRAPHIC  DATA  BASE 


XLeft  X  Right  Y  Upper 


0  20  0 

20  40  0 

40  60  0 

60  80  0 

80  100  0 

100  120  0 

120  140  0 

140  160  0 

160  180  0 

180  200  0 

200  220  0 

220  240  0 

240  260  0 

260  280  0 

280  300  0 

300  320  0 

320  340  0 

340  360  0 

360  380  0 

380  400  0 

400  420  0 

420  440  0 

440  460  0 

460  480  0 

480  500  0 

500  520  0 

520  540  0 

540  ,560  0 

560  580  0 

580  600  0 

600  620  0 

620  640  0 

640  660  0 

660  680  0 

680  700  0 

0  20  20 

20  40  20 

40  60  20 

.  60  80  20 

80  100  20 

100  120  20 
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140  160  20 

160  180  20 

180  200  20 
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140 

160 

220 

240 

2 

1 

160 

180 

220 

240 

2 

1 

180 

200 

220 

240 

2 

1 

200 

220 

220 

240 

2 

1 

220 

240 

220 

240 

2 

1 

240 

260 

220 

240 

2 

1 

260 

280 

220 

240 

2 

1 

280 

300 

220 

240 

2 

1 

300 

320 

220 

240 

2 

1 

320 

340 

220 

240 

2 

1 

340 

360 

220 

240 

2 

1 

360 

380 

220 

240 

2 

1 

380 

400 

220 

240 

2 

1 

400 

420 

220 

240 

2 

1 

64 


65 


40 

60 

260 

280 

1 

1 

60 

80 

260 

280 

1 

1 

80 

100 

260 

280 

1 

1 

100 

120 

260 

280 

2 

1 

120 

140 

260 

280 

1 

1 

140 

160 

260 

280 

1 

1 

160 

180 

260 

280 

1 

1 

180 

200 

260 

280 

2 

1 

200 

220 

260 

280 

2 

1 

220 

240 

260 

280 

2 

1 

240 

260 

260 

280 

2 

1 

260 

280 

260 

280 

2 

1 

280 

300 

260 

280 

2 

1 

300 

320 

260 

280 

2 

1 

320 

340 

260 

280 

2 

1 

340 

360 

260 

280 

2 

1 

360 

380 

260 

280 

2 

1 

380 

400 

260 

280 

1 

1 

400 

420 

260 

280 

1 

1 

420 

440 

260 

280 

2 

1 

440 

460 

260 

280 

2 

1 

460 

480 

260 

280 

2 

1 

480 

500 

260 

280 

1 

1 

500 

520 

260 

280 

1 

1 

520 

540 

260 

280 

1 

1 

540 

560 

260 

280 

1 

1 

560 

580 

260 

280 

1 

1 

580 

600 

260 

280 

1 

1 

600 

620 

260 

280 

1 

1 

620 

640 

260 

280 

1 

1 

640 

660 

260 

280 

1 

1 

660 

680 

260 

280 

1 

1 

680 

700 

260 

280 

1 

1 

0 

20 

280 

300 

1 

1 

20 

40 

280 

300 

1 

1 

40 

60 

280 

300 

1 

1 

60 

80 

280 

300 

1 

1 

80 

100 

280 

300 

1 

1 

100 

120 

280 

300 

1 

1 

120 

140 

280 

300 

1 

1 

140 

160 

280 

300 

1 

1 

160 

180 

280 

300 

2 

1 

180 

200 

280 

300 

2 

1 

200 

220 

280 

300 

2 

1 

220 

240 

280 

300 

2 

1 

240 

260 

280 

300 

2 

1 

260 

280 

280 

300 

2 

1 

280 

300 

280 

300 

2 

1 

300 

320 

280 

300 

2 

1 

320 

340 

280 

300 

2 

1 

340 

360 

280 

300 

2 

1 

66 


67 


680 

700 

300 

320 

2 

1 

0 

20 

320 

340 

1 

1 

20 

40 

320 

340 

1 

1 

40 

60 

320 

340 

1 

1 

60 

80 

320 

340 

1 

1 

80 

100 

320 

340 

1 

1 

100 

120 

320 

340 

1 

1 

120 

140 

320 

340 

1 

1 

140 

160 

320 

340 

1 

1 

160 

180 

320 

340 

1 

1 

180 

200 

320 

340 

2 

1 

200 

220 

320 

340 

2 

1 

220 

240 

320 

340 

2 

1 

240 

260 

320 

340 

2 

1 

260 

280 

320 

340 

2 

1 

280 

300 

320 

340 

2 

1 

300 

320 

320 

340 

2 

1 

320 

340 

320 

340 

2 

1 

340 

360 

320 

340 

1 

1 

360 

380 

320 

340 

1 

1 

380 

400 

320 

340 

1 

1 

400 

420 

320 

340 

1 

1 

420 

440 

320 

340 

2 

1 

440 

460 

320 

340 

2 

1 

460 

480 

320 

340 

2 

1 

480 

500 

320 

340 

2 

1 

500 

520 

320 

340 

2 

1 

520 

540 

320 

340 

2 

1 

540 

560 

320 

340 

2 

1 

560 

580 

320 

340 

2 

1 

580 

600 

320 

340 

2 

1 

600 

620 

320 

340 

2 

1 

620 

640 

320 

340 

2 

1 

640 

660 

320 

340 

2 

1 

660 

680 

320 

340 

2 

1 

680 

700 

320 

340 

2 

1 

0 

20 

340 

360 

1 

1 

20 

40 

340 

360 

1 

1 

40 

60 

340 

360 

1 

1 

60 

80 

340 

360 

1 

1 

80 

100 

340 

360 

1 

1 

100 

120 

340 

360 

1 

1 

120 

140 

340 

360 

1 

1 

140 

160 

340 

360 

1 

1 

160 

180 

340 

360 

1 

1 

180 

200 

340 

360 

2 

1 

200 

220 

340 

360 

2 

1 

220 

240 

340 

360 

2 

1 

240 

260 

340 

360 

2 

1 

260 

280 

340 

360 

2 

1 

280 

300 

340 

360 

2 

1 

68 


300 

320 

340 

360 

2 

1 

320 

340 

340 

360 

2 

1 

340 

360 

340 

360 

1 

1 

360 

380 

340 

360 

1 

1 

380 

400 

340 

360 

1 

1 

400 

420 

340 

360 

1 

1 

420 

440 

340 

360 

1 

1 

440 

460 

340 

360 

2 

1 

460 

480 

340 

360 

2 

1 

480 

500 

340 

360 

2 

1 

500 

520 

340 

360 

2 

1 

520 

540 

340 

360 

2 

1 

540 

560 

340 

360 

2 

1 

560 

580 

340 

360 

2 

1 

580 

600 

340 

360 

2 

1 

600 

620 

340 

360 

2 

1 

620 

640 

340 

360 

2 

1 

640 

660 

340 

360 

2 

1 

660 

680 

340 

360 

2 

1 

680 

700 

340 

360 

2 

1 

0 

20 

360 

380 

1 

1 

20 

40 

360 

380 

1 

1 

40 

60 

360 

380 

1 

1 

60 

80 

360 

380 

1 

1 

80 

100 

360 

380 

1 

1 

100 

120 

360 

380 

1 

1 

120 

140 

360 

380 

1 

1 

140 

160 

360 

380 

1 

1 

160 

180 

360 

380 

1 

1 

180 

200 

360 

380 

2 

1 

200 

220 

360 

380 

2 

1 

220 

240 

360 

380 

2 

1 

240 

260 

360 

380 

2 

1 

260 

280 

360 

380 

2 

1 

280 

300 

360 

380 

2 

1 

300 

320 

360 

380 

2 

1 

320 

340 

360 

380 

2 

1 

340 

360 

360 

380 

1 

1 

360 

380 

360 

380 

1 

1 

380 

400 

360 

380 

1 

1 

400 

420 

360 

380 

1 

1 

420 

440 

360 

380 

1 

1 

440 

460 

360 

380 

2 

1 

460 

480 

360 

380 

2 

1 

480 

500 

360 

380 

2 

1 

500 

520 

360 

380 

2 

1 

520 

540 

360 

380 

2 

1 

540 

560 

360 

380 

2 

1 

560 

580 

360 

380 

2 

1 

580 

600 

360 

380 

2 

1 

600 

620 

360 

380 

2 

1 

69 


620 

640 

360 

380 

2 

1 

640 

660 

360 

380 

2 

1 

660 

680 

360 

380 

2 

1 

680 

700 

360 

380 

2 

1 

0 

20 

380 

400 

1 

1 

20 

40 

380 

400 

1 

1 

40 

60 

380 

400 

1 

1 

60 

80 

380 

400 

1 

1 

80 

100 

380 

400 

1 

1 

100 

120 

380 

400 

1 

1 

120 

140 

380 

400 

1 

1 

140 

160 

380 

400 

1 

1 

160 

180 

380 

400 

1 

1 

180 

200 

380 

400 

2 

1 

200 

220 

380 

400 

2 

1 

220 

240 

380 

400 

2 

1 

240 

260 

380 

400 

2 

1 

260 

280 

380 

400 

2 

1 

280 

300 

380 

400 

2 

1 

300 

320 

380 

400 

2 

1 

320 

340 

380 

400 

2 

1 

340 

360 

380 

400 

1 

1 

360 

380 

380 

400 

1 

1 

380 

400 

380 

400 

1 

1 

400 

420 

380 

400 

1 

1 

420 

440 

380 

400 

1 

1 

440 

460 

380 

400 

1 

1 

460 

480 

380 

400 

2 

1 

480 

500 

380 

400 

2 

1 

500 

520 

380 

400 

2 

1 

520 

540 

380 

400 

2 

1 

540 

560 

380 

400 

2 

1 

560 

580 

380 

400 

2 

1 

580 

600 

380 

400 

2 

1 

600 

620 

380 

400 

2 

1 

620 

640 

380 

400 

2 

1 

640 

660 

380 

400 

2 

1 

660 

680 

380 

400 

2 

1 

680 

700 

380 

400 

2 

1 

70 
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lit# 

Description 

HitstoKiU 

Unit  Type 

Op  Speed 

Not  Used 

#  Missiles  or  A/C 

Eff.  Range 

1 

U.S.  CV 

40 

1 

20 

0 

46 

550 

2 

U.S.  warship 

10 

2 

20 

0 

40 

600 

3 

U.S.  warship 

10 

2 

20 

0 

40 

600 

4 

U.S.  warship 

10 

2 

20 

0 

40 

600 

5 

U.S.  warship 

10 

2 

20 

0 

40 

600 

6 

U.S.  warship 

10 

2 

20 

0 

40 

600 

7 

AOR 

8 

3 

20 

0 

0 

0 

8 

U.S.  warship 

10 

2 

15 

0 

40 

600 

9 

U.S.  warship 

10 

2 

15 

0 

40 

600 

10 

merchant 

8 

3 

15 

0 

0 

0 

11 

merchant 

8 

3 

15 

0 

0 

0 

12 

E-2C 

4 

12 

100 

0 

0 

0 

13 

U.S.  recon  a/c 

2 

13 

325 

0 

0 

0 

14 

U.S.  recon  a/c 

2 

13 

325 

0 

0 

0 

15 

U.S.  recon  a/c 

2 

13 

200 

0 

0 

0 

16 

Satellite 

400 

19 

150 

0 

0 

0 

17 

Satellite 

400 

19 

150 

0 

0 

0 

18 

Satellite 

400 

19 

50 

0 

0 

0 

19 

Satellite 

400 

19 

50 

0 

0 

0 

20 

U.S.  CV 

40 

1 

20 

0 

46 

550 

21 

Iranian  air  base 

60 

21 

0 

0 

50 

230 

22 

Iranian  air  base 

60 

21 

0 

0 

50 

230 

23 

Iranian  radar 

8 

22 

0 

0 

0 

0 

24 

Iranian  radar 

8 

22 

0 

0 

0 

0 

25 

Iranian  radar 

8 

22 

0 

0 

0 

0 

26 

Iranian  radar 

8 

22 

0 

0 

0 

0 

27 

Iranian  SSM 

1 

24 

15 

0 

16 

50 

28 

Iranian  SSM 

1 

24 

15 

0 

16 

50 

29 

Iranian  SSM 

1 

24 

15 

0 

16 

50 

30 

Iranian  SSM 

1 

24 

15 

0 

16 

50 

31 

Iranian  SSM 

1 

24 

15 

0 

16 

50 

32 

Iranian  SSM 

1 

24 

15 

0 

16 

50 

33 

Iranian  SSM 

1 

24 

15 

0 

16 

50 

34 

Iranian  SSM 

1 

24 

15 

0 

16 

50 

35 

Iranian  SSM 

1 

24 

15 

0 

16 

50 

36 

Iranian  SSM 

1 

24 

15 

0 

16 

50 

37 

Iranian  SSM 

1 

24 

15 

0 

16 

50 

38 

Iranian  recon  a/c 

6 

33 

600 

0 

0 

0 

39 

Iranian  recon  a/c 

6 

33 

600 

0 

0 

0 

40 

Iranian  recon  a/c 

6 

33 

600 

0 

0 

0 

Unit  Type:  1  =  aircraft  carrier,  2  =  U.S.  warship,  3  =  merchant,  12  =  E2-C  aircraft,  13  = 

U.S.reconaissance  aircraft,  19  =  U.S.  satellite,  21  =  Iranian  air  base,  22  =  Iranian  long  range  radar  site,  24 

=  Iranian  SSM  launcher,  33  =  Iranian  reconnaissance  aircraft 

Description:  Included  in  the  Appendix  for  clarity,  it  is  not  in  the  actual  input  data  file 

Op  Speed:  Normal  operating  speed  of  the  unit  in  nautical  miles  per  hour  for  ships  and  aircraft,  miles  per 

hour  for  SSM  launchers. 

Missiles  or  A/C:  Number  of  surface-to-land  missiles  on  board  a  warship,  surface-to-surface  missiles  for 
an  SSM  launcher,  or  attack  aircraft  on  board  a  carrier  or  air  base. 

Eff.  Range:  The  effective  range  of  the  missile  or  aircraft  in  miles,  used  as  a  maximum  range. 
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Unit#  Description 

#  Fighter  A/C 

Eff.  Range 

Surface  sens 

Air  sens 

ESM  sen 

X(0) 

Y(0) 

1 

U.S.  CV 

30 

200 

2 

8 

2 

695 

395 

2 

U.S.  warship 

0 

0 

2 

8 

2 

695 

393 

3 

U.S.  warship 

0 

0 

2 

8 

2 

694 

393 

4 

U.S.  warship 

0 

0 

2 

8 

2 

696 

393 

5 

U.S.  warship 

0 

0 

2 

8 

2 

694 

395 

6 

U.S.  warship 

0 

0 

2 

8 

2 

696 

395 

7 

AOR 

0 

0 

2 

3 

0 

695 

397 

8 

U.S.  warship 

0 

0 

2 

8 

2 

250 

360 

9 

U.S.  warship 

0 

0 

2 

8 

2 

252 

362 

10 

merchant 

0 

0 

2 

3 

0 

249 

362 

11 

merchant 

0 

0 

2 

3 

0 

250 

362 

12 

E-2C 

0 

0 

4 

6 

6 

695 

395 

13 

U.S.  recon  a/c 

0 

0 

2 

4 

2 

695 

395 

14 

U.S.  recon  a/c 

0 

0 

2 

4 

2 

695 

395 

15 

U.S.  recon  a/c 

0 

0 

2 

4 

2 

695 

395 

16 

Satellite 

0 

0 

4 

6 

4 

699 

200 

17 

Satellite 

0 

0 

4 

6 

4 

699 

250 

18 

Satellite 

0 

0 

4 

6 

4 

400 

399 

19 

Satellite 

0 

0 

4 

6 

4 

500 

399 

20 

Satellite 

30 

200 

2 

4 

2 

695 

396 

21 

Iranian  air  base 

30 

335 

2 

4 

4 

410 

135 

22 

Iranian  air  base 

50 

335 

4 

4 

4 

610 

250 

23 

Iranian  radar 

0 

0 

4 

4 

4 

340 

210 

24 

Iranian  radar 

0 

0 

4 

4 

4 

405 

175 

25 

Iranian  radar 

0 

0 

4 

4 

4 

490 

250 

26 

Iranian  radar 

0 

0 

4 

4 

4 

650 

295 

27 

Iranian  SSM 

0 

0 

2 

4 

2 

290 

200 

28 

Iranian  SSM 

0 

0 

2 

4 

2 

325 

205 

29 

Iranian  SSM 

0 

0 

2 

4 

2 

380 

205 

30 

Iranian  SSM 

0 

0 

2 

4 

2 

390 

202 

31 

Iranian  SSM 

0 

0 

2 

4 

2 

410 

195 

32 

Iranian  SSM 

0 

0 

2 

4 

2 

425 

195 

33 

Iranian  SSM 

0 

0 

2 

4 

2 

450 

215 

34 

Iranian  SSM 

0 

0 

2 

4 

2 

450 

230 

35 

Iranian  SSM 

0 

0 

2 

4 

2 

470 

250 

36 

Iranian  SSM 

0 

0 

2 

4 

2 

485 

265 

37 

Iranian  SSM 

0 

0 

2 

4 

2 

500 

270 

38 

Iranian  recon  a/c 

0 

0 

2 

4 

2 

610 

250 

39 

Iranian  recon  a/c 

0 

0 

2 

4 

2 

610 

250 

40 

Iranian  recon  a/c 

0 

0 

2 

4 

2 

410 

135 

Fighter  A/C:  Initial  number  of  fighter  aircraft  allocated  to  a  U.S.  CV  or  an  Iranian  air  base. 
Eff.  Range:  The  effective  range  of  the  fighter  aircraft  in  miles,  used  as  a  maximiun  range. 
Surface  sens:  Siuface  sensor  range  in  grid  units,  1  unit  =  20  miles. 

Air  sens:  Air  sensor  range  in  grid  units,  1  unit  =  20  miles. 

ESM  sens:  Electronic  supports  measures  range  in  grid  units,  1  unit  =  20  miles. 

X(0),  Y(0):  Initial  (time  =  0)  x  and  y  positions  of  a  unit  in  cartesian  coordinates. 
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Unit#  Description 

X(l) 

Y(l)  Loiter(l)  X(2) 

Y(2)  Loiter(2) 

X(3) 

Y(3)  Loiter(3) 

X(4) 

Y(4) 

Loiter(4) 

1 

U  S.  CV 

460 

395 

0 

430 

330 

720 

460 

395 

720 

595 

395 

120 

2 

U.S.  warship 

460 

393 

0 

430 

328 

720 

460 

393 

720 

595 

393 

120 

3 

U.S.  warship 

459 

393 

0  • 

429 

328 

720 

459 

393 

720 

594 

393 

120 

4 

U.S.  warship 

461 

393 

0 

431 

328 

720 

461 

393 

720 

596 

393 

120 

5 

U.S.  warship 

459 

395 

0 

429 

330 

720 

459 

395 

720 

594 

395 

120 

6 

U.S.  warship 

461 

395 

0 

431 

330 

720 

461 

395 

720 

596 

395 

120 

7 

AOR 

460 

397 

0 

430 

332 

720 

460 

397 

720 

595 

397 

120 

8 

U.S.  warship 

330 

375 

0 

330 

315 

0 

397 

237 

0 

422 

237 

0 

9 

U.S.  warship 

332 

377 

0 

332 

315 

0 

397 

239 

0 

422 

239 

0 

10 

merchant 

329 

377 

0 

329 

315 

0 

397 

239 

0 

422 

239 

0 

11 

merchant 

330 

377 

0 

330 

315 

0 

397 

239 

0 

422 

239 

0 

12 

E-2C 

630 

350 

0 

500 

350 

0 

695 

395 

0 

695 

395 

10 

13 

U.S.  recon  a/c 

680 

295 

0 

640 

295 

0 

695 

395 

0 

695 

395 

15 

14 

U.S.  recon  a/c 

640 

295 

0 

590 

295 

0 

695 

395 

0 

695 

395 

20 

15 

U.S.  recon  a/c 

590 

295 

0 

540 

295 

0 

695 

395 

0 

695 

395 

25 

16 

Satellite 

10 

200 

240 

699 

200 

240 

10 

200 

240 

699 

200 

240 

17 

Satellite 

10 

250 

120 

699 

250 

240 

10 

250 

240 

699 

250 

240 

18 

Satellite 

400 

10 

120 

400 

399 

120 

400 

10 

120 

400 

399 

120 

19 

Satellite 

500 

10 

240 

500 

399 

120 

500 

10 

120 

500 

399 

120 

20 

Satellite 

460 

396 

0 

430 

331 

0 

460 

396 

0 

595 

396 

120 

21 

Iranian  air  base 

410 

135 

9000 

410 

135 

9000 

410 

135 

9000 

410 

135 

9000 

22 

Iranian  air  base 

610 

250 

9000 

610 

250 

9000 

610 

250 

9000 

610 

250 

9000 

23 

Iranian  radar 

340 

210 

9000 

340 

210 

9000 

340 

210 

9000 

340 

210 

9000 

24 

Iranian  radar 

405 

175 

9000 

405 

175 

9000 

405 

175 

9000 

405 

175 

9000 

25 

Iranian  radar 

490 

250 

9000 

490 

250 

9000 

490 

250 

9000 

490 

250 

9000 

26 

Iranian  radar 

650 

295 

9000 

650 

295 

9000 

650 

295 

9000 

650 

295 

9000 

27 

Iranian  SSM 

300 

150 

60 

300 

150 

60 

300 

150 

60 

300 

150 

60 

28 

Iranian  SSM 

300 

150 

60 

300 

150 

60 

300 

150 

60 

300 

150 

60 

29 

Iranian  SSM 

300 

150 

60 

300 

150 

60 

300 

150 

60 

300 

150 

60 

30 

Iranian  SSM 

300 

150 

60 

300 

150 

60 

300 

.150 

60 

300 

150 

60 

31 

Iranian  SSM 

300 

150 

60 

300 

150 

60 

300 

150 

60 

300 

150 

60 

32 

Iranian  SSM 

300 

150 

60 

300 

150 

60 

300 

150 

60 

300 

150 

60 

33 

Iranian  SSM 

550 

200 

60 

550 

200 

60 

550 

200 

60 

550 

200 

60 

34 

Iranian  SSM 

550 

200 

60 

550 

200 

60 

550 

200 

60 

550 

200 

60 

35 

Iranian  SSM 

550 

200 

60 

550 

200 

60 

550 

200 

60 

550 

200 

60 

36 

Iranian  SSM 

550 

200 

60 

550 

200 

60 

550 

200 

60 

550 

200 

60 

37 

Iranian  SSM 

550 

200 

60 

550 

200 

60 

550 

200 

60 

550 

200 

60 

38 

Iranian  recon  a/c 

665 

345 

0 

565 

345 

0 

610 

250 

0 

610 

250  ■ 

30 

39 

Iranian  recon  a/c 

505 

385 

0 

555 

385 

0 

610 

250 

0 

610 

250 

15 

40 

Iranian  recon  a/c 

475 

295 

0 

430 

295 

0 

410 

135 

0 

410 

135 

30 

X(1),Y(1):  X  and  y  coordinates  of  the  first  waypoint  for  all  units  exeptstationaiy  and  SSM  units.  For 

stationary  units  (Iranian  air  base  and  radar  sites)  all  waypoints  are  equal  to  X(0)  and  Y(0).  For  an  SSM 
launcher  this  is  the  location  of  the  reload  site. 

Loiter(l):  Amount  of  time  units  remain  at  the  first  waypoint,  prior  to  moving  towards  the  2nd  waypoint. 
For  stationary  units  this  number  is  greater  than  the  total  time  of  the  battle.  For  SSM  units  this  is  the  time 
it  remains  at  the  reload  site  replenishing  its  missile  supply. 

X(2),  Y(2),  X(3),  Y(3),  X(4),  Y(4):  X  and  y  coordinates  of  the  second,  third,  and  fourth  waypoints  for 
all  units  exept  stationary  and  SSM  units.  No  meaning  for  stationary  and  SSM  units. 

Loiter(2),  Loiter(3),  Loiter(4):  Amount  of  time  units  remain  at  the  second,  thrid,  or  fourth  waypoints. 
No  meaning  for  stationary  and  SSM  units.  For  a  reconnaissance  aircraft,  Loiter(4)  it  is  the  amount  of 
time  it  remains  at  home  base  imtil  commencing  another  reconnaissance  flight 
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Unit#  Description 

X{5)  Y(5)  Loiter(5) 

X(6) 

Y(6)  Loiter(6) 

1 

U.S.  CV 

695 

395 

1440 

699 

399 

9000 

2 

U.S.  warship 

695 

393 

1440 

699 

399 

9000 

3 

U.S.  warship 

694 

393 

1440 

699 

399 

9000 

4 

U.S.  warship 

696 

393 

1440 

699 

399 

9000 

5 

U.S.  warship 

694 

395 

1440 

699 

399 

9000 

6 

U.S.  warship 

696 

395 

1440 

699 

399 

9000 

7 

AOR 

695 

397 

1440 

699 

399 

9000 

8 

U.S.  warship 

460 

370 

0 

690 

399 

9000 

9 

U.S.  warship 

460 

370 

0 

690 

399 

9000 

10 

merchant 

460 

370 

0 

690 

399 

9000 

11 

merchant 

460 

370 

0 

690 

399 

9000 

12 

E-2C 

55 

45 

1 

130 

45 

1 

13 

U.S.  recon  a/c 

15 

100 

1 

40 

100 

1 

14 

U.S.  recon  a/c 

55 

100 

1 

50 

100 

1 

15 

U.S.  recon  a/c 

105 

100 

1 

50 

100 

9000 

16 

Satellite 

10 

200 

240 

699 

200 

9000 

17 

Satellite 

10 

250 

240 

699 

250 

9000 

18 

Satellite 

400 

10 

120 

400 

399 

9000 

19 

Satellite 

500 

10 

120 

500 

399 

9000 

20 

Satellite 

695 

396 

1440 

699 

398 

9000 

21 

Iranian  air  base 

410 

135 

9000 

410 

135 

9000 

22 

Iranian  air  base 

610 

250 

9000 

610 

250 

9000 

23 

Iranian  radar 

340 

210 

9000 

340 

210 

9000 

24 

Iranian  radar 

405 

175 

9000 

405 

175 

9000 

25 

Iranian  radar 

490 

250 

9000 

490 

250 

9000 

26 

Iranian  radar 

650 

295 

9000 

650 

295 

9000 

27 

Iranian  SSM 

300 

150 

60 

300 

150 

9000 

28 

Iranian  SSM 

300 

150 

60 

300 

150 

9000 

29 

Iranian  SSM 

300 

150 

60 

300 

150 

9000 

30 

Iranian  SSM 

300 

150 

60 

300 

150 

9000 

31 

Iranian  SSM 

300 

150 

60 

300 

150 

9000 

32 

Iranian  SSM 

300 

150 

60 

300 

150 

9000 

33 

Iranian  SSM 

550 

200 

60 

550 

200 

9000 

34 

Iranian  SSM 

550 

200 

60 

550 

200 

9000 

35 

Iranian  SSM 

550 

200 

60 

550 

200 

9000 

36 

Iranian  SSM 

550 

200 

60 

550 

200 

9000 

37 

Iranian  SSM 

550 

200 

60 

550 

200 

9000 

38 

Iranian  recon  a/c 

55 

95 

22 

100 

95 

1 

39 

Iranian  recon  a/c 

105 

135 

22 

50 

135 

-1 

40 

Iranian  recon  a/c 

65 

160 

21 

45 

160 

1 

X(5),  Y(5):  X  and  y  coordinates  of  naval  units  fifth  waypoint.  No  meaning  for  stationary  and  SSM  units. 
For  reconnaissance  aircraft,  these  are  the  x  and  y  distances,  between  its  current  location  and  its  x  and  y 
destination  on  the  first  leg  of  a  reconnaissance  flight. 

Loiter(5):  Time  naval  units  remain  at  the  fifth  waypoint.  No  meaning  for  stationary  and  SSM  units.  For 
reconnaissance  aircraft,  this  is  the  Unit#  of  the  air  base  or  aircraft  carrier  that  it  is  attached. 

X(6),  Y(6):  X  and  y  coordinates  of  naval  units  sixth  waypoint.  No  meaning  for  stationary  and  SSM 
units.  For  reconnaissance  aircraft,  these  are  the  x  and  y  distances,  between  its  first  leg  endpoint  and  and 
its  X  and  y  destination  on  the  second  leg  of  a  recoimaissance  flight. 

Loiter(6);  Amount  of  time  naval  units  will  remain  at  the  sixth  waypoint.  For  stationary,  SSM,  and  U.S. 
reconnaissance  aircraft  units  it  has  no  meaning.  For  Iranian  recoimaissance  aircraft,  this  identifies  the 
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direction  from  which  the  aircraft  will  fly  away  from  the  air  base  during  its  first  leg  of  a  recoimaissance 
flight. 
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APPENDIX  C.  ALGORITHMS 


1.  NAVAL  SHIP  MOVEMENT 

IS  THE  SHIP  IN  MOTION? 

NO  >  LOITER  TIME  COMPLETE? 

NO  >  EXIT  algorithm 

YES>  Compute  course  &  speed  to  next  waypoint 
Determine  time  to  reach  next  w^oint 

{Time  =  distance  between  waypoints/speed} 

Place  ship  in  motion 
EXIT  algorithm 

YES>  IS  UNIT  AT  WAYPOINT? 

YES>  WILL  SHIP  LOITER  AT  THIS  POINT? 

YES>  Stop  ship  motion 

Compute  time  at  which  ship  will  start  movement 
NO>  Compute  course  &  speed  to  next  waypoint 
Determine  time  to  reach  next  waypoint 
Place  ship  in  motion 
EXIT  algorithm 

2.  MOBILE  SSM  LAUNCHER  MOVEMENT 

IS  THE  CURRENT  WAYPOINT  THE  RELOAD  SITE? 

YES>  IS  THE  SSM  IN  TRANSIT? 

YES>HAS  IT  REACHED  THE  SUPPLY  DEPOT  THIS  TIME  STEP? 
YES  >  Stop  SSM  motion 

Replenish  missile  supply 
Compute  time  at  which  SSM  will  start  movement 
{Time  =  current  time  +  loiter  time} 

NO>  Update  SSM  location 
EXIT  algorithm 

NO>LOITER  TIME  COMPLETE? 

NO  >  EXIT  algorithm 

YES>  Randomly  compute  next  Waypoint, 

Compute  course  &  speed  to  next  waypoint 
Determine  time  to  reach  next  waypoint 
Place  SSM  in  motion 
EXIT  algorithm 

NO>  DOES  THE  LAUNCHER  HAVE  MISSILES? 

NO  >  Compute  course  &  speed  to  reload  site 
Determine  time  to  reach  reload  site 
Place  SSM  in  motion 
Discontinue  emitting 
EXIT  algorithm 

YES>  IS  IT  ON  THE  HOUR  OR  HALF  HOUR? 
ISU(0,l)<p„<„.? 

NO  >  Use  Bernoulli  trial  to  determine  if  emitting 
EXIT  algorithm 

YES>Randomly  compute  next  waypoint 
Discontinue  emitting, 

Compute  course  &  speed  to  next  waypoint 
Determine  time  to  reach  next  waypoint 
Place  SSM  in  motion 
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EXIT  algorithm 

3.  RECONNAISSANCE  AIRCRAFT  MOVEMENT 

IS  THE  CURRENT  WAYPOINT  NUMBER  1? 

YES>  IS  THE  AIRCRAFT  IN  TRANSIT? 

NO>  Compute  course  &  speed  to  waypoint  1 
Determine  time  to  reach  waypoint  1 
YES>  HAS  THE  A/C  REACHED  WAYPOINT  1? 

YES>  Compute  course  &  speed  to  waypoint  2 
Determine  time  to  reach  waypoint  2 
GOTO  next  step 
NO>  GOTO  next  step 
Place  aircraft  in  motion 
Commence  emitting 
EXIT  algorithm 

NO>  IS  THE  CURRENT  WAYPOINT  NUMBER  2? 

YES>  HAS  THE  AIRCRAFT  REACHED  THIS  WAYPOINT? 
YES>  Compute  course  &  speed  to  waypoint  3 
Determine  time  to  reach  waypoint  3 
GOTO  next  step 
NO>  GOTO  next  step 
Place  aircraft  in  motion 
EXIT  algorithm 

NO>  IS  THE  CURRENT  WAYPOINT  NUMBER  3? 

YES>  HAS  THE  AIRCRAFT  REACHED  THIS  WP? 

YES>  Update  its  x,  y  location  to  the  CV  or  air  base 
Stop  aircraft  motion 

Determine  time  at  which  next  flight  commences 
EXIT  algorithm 

NO>  Compute  course  &  speed  to  CV  or  air  base 
Determine  time  to  reach  home  base 
Place  aircraft  in  motion 
EXIT  algorithm 

NO>  IS  THE  A/C  LOITERING  AT  ITS  HOME  BASE? 
YES>  IS  LOITER  TIME  COMPLETE? 

YES>  Compute  new  waypoint  1 
NO>  GOTO  next  step 
Update  its  x,  y  location  to  CV  or  air  base 
EXIT  algoritW 
NO>  EXIT  algorithm 

4.  ENGAGEMENT  LIST  (threat  determination  and  closest  unit) 

Cycle  through  all  U.S.  contacts  listed  on  the  Iranian  intelligence  estimate 
Beg 

Determine  the  type  of  U.S.  unit 

Cycle  through  all  Iranian  units  to  determine  which  Iranian  units  are  threatened 
by  the  U.S.  contact,  and  which  Iranian  unit  will  engage  the  U.S.  contact 
Beg 

Determine  distance  between  U.S.  &  Iranian  unit 
IS  DISTANCE  <  EFFECTIVE  WEPS  RANGE  OF  U.S.  UNIT? 
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YES>  Classify  as  high  threat 
N0>  Classify  as  low  threat 
IF  UNITS  ARE  AIR  BASE  AND  NAVAL  SHIP? 

Beg 

IS  DISTANCE  <  CURRENT  MIN  AIR  DISTANCE? 
YES>  Assigned  attack  unit  =  unit 
Min  air  distance  =  distance 
End 

F  UNITS  ARE  SSM  AND  NAVAL  SHIP? 

Beg 

IS  DISTANCE  <  CURRENT  MIN  MISSILE  DISTANCE? 
YES>  Assigned  missile  unit  =  unit 

Min  missile  distance  =  distance 

End 

F  UNITS  ARE  AIR  BASE  AND  U.S.  RECON  PLANE? 

Beg 

IS  DISTANCE  <  CURRENT  MIN  AIR  DISTANCE? 
YES>  Assigned  air  unit = unit 

Min  air  distance  =  distance 
End 

WAS  A  MISSILE  UNIT  ASSIGNED? 

YES>  GOTO  wejpjn  selection  algorithm 
NO>  WAS  AN  AIR  UNIT  ASSIGNED? 

YES>  GOTO  we^n  selection  algorithm 
NO>  GOTO  beginning  of  algorithm 
WAS  A  WEAPON  ASSIGNED  TO  STRIKE  THE  ENEMY? 

YES>  Conqjute  conflict  adjudication  time 
NO  >  GOTO  beginning  of  algorithm 

End 


End 

Cycle  through  all  Iranian  contacts  listed  on  the  U.S.  intelligence  estimate 
Beg 

Determine  the  type  of  Iranian  unit 

Cycle  through  all  U.S.  units  to  detmiine  wliich  U.S.  units  are  threatened 
by  the  Iranian  contact,  and  which  U.S.  unit  will  engage  the  Iranian  contact 
Beg 

Determine  distance  between  Iranian  &  U.S.  unit 
IS  DISTANCE  <  EFFECTIVE  WEPS  RANGE  OF  IRAN  UNIT? 
YES>  Classify  as  high  threat 
NO>  IS  THE  IRANIAN  AIR  BASE  OR  RADAR? 
YES>  Classify  as  medium  threat 
NO>  Classify  as  high  threat 
F  UNITS  ARE  CV  AND  LAND  BASED  UNIT? 

Beg 

IS  DISTANCE  <  CURRENT  MIN  AIR  DISTANCE? 
YES>  Assigned  attack  unit = unit 
Mm  air  distance  =  distance 
End 

F  UNITS  ARE  WARSHF  AND  LAND  BASED  UNIT? 

Beg 
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IS  DISTANCE  <  CURRENT  MN  MISSILE  DISTANCE? 
YES>  Assigned  missile  unit = unit 

Min  missile  distance  =  distance 

End 

IF  UNITS  ARE  C V  AND  IRANIAN  RECON  PLANE? 

Beg 

IS  DISTANCE  <  CURRENT  MIN  AIR  DISTANCE? 

YES>  Assigned  air  unit  =  unit 

Mn  air  distance  =  distance 
End 

GOTO  we^n  selection  algoritlim 

WAS  A  WEAPON  ASSIGNED  TO  STRIKE  THE  IRANIAN  UNIT? 
YES>  Conq)ute  conflict  adjudication  time 
NO  >  GOTO  beginning  of  algorithm 

End 

End 


5.  IRANIAN  WEAPON  SELECTION  AGAINST  UNITED  STATES  UNITS 

Determine  the  type  of  target 

IS  THE  TARGET  A  RECONNAISSANCE  AIRPLANE? 

YES>  Identify  the  air  base  that  is  closest  to  the  target  and  has  fighter 
aircraft  available  and  within  range  of  the  target 

NO>  1.  Identify  flie  SSM  site  that  is  closest  to  the  target  and  has  missiles 

available  and  within  range  to  strike  the  target 

2.  Identify  the  air  base  that  is  closest  to  the  target  and  has  attack 
aircraft  available  that  are  within  range  to  strike  the  target 
ARE  ANY  OF  THE  IRANIAN  UNITS  WITHIN  THE  RANGE  OF  ANY  WEAPON 
CARRIED  ONBOARD  THE  TARGET? 

YES>  Classify  the  target  as  a  high  threat 
NO>  Classify  the  target  as  a  low  threat 

IF  TARGET  IS  A/C  CARRIER? 

Beg 

IF  CLOSEST  ASSET  IS  AN  AIR  BASE 
IF  THREAT  IS  HlCai 

Assign  Ha  attack  aircraft 

ELSE 

Assign  La  attack  aircraft 
IF  CLOSEST  ASSET  IS  SSM  SITE 

Launch  all  missiles  left  on  launcher  at  A/C  carrier 

End 

IF  TARGET  IS  A  SHIP 
Beg 

IF  CLOSEST  ASSET  IS  AN  AIR  BASE 
IF  THREAT  IS  HIGH 

Assign  Hs  attack  aircraft 

ELSE 

Assign  Ls  attack  aircraft 
IF  CLOSEST  ASSET  IS  AN  SSM  SITE 
Launch  Ms  missiles  at  the  ship 
End 

IF  THE  TARGET  IS  AN  AIRPLANE 


80 


Beg 

Assign  Hf  figjiters  to  ragage  the  aircraft 

End 

6.  UNmD  STATES  WEAPON  SEUlCnONAGAmST  IRANIAN  UNITS 

Determine  the  type  of  target 

determine  the  ftiendly  combat  unit  that  is  closest  to  the  target 
Classify  the  target  as  a  high  or  low  threat 
F  THE  TARGET  IS  AN  AIR  BASE 
Beg 

F  THE  THREAT  IS  HIGH 

Assign  Ha  attack  aircraft 

ELSE 

Launch  Ma  cruise  missiles  at  the  target 

End 

F  THE  TARGET  IS  A  RADAR  SITE 
Beg 

F  THE  THREAT  IS  HIGH 

Assign  Hr  attack  aircraft 

ELSE 

Launch  Mr  cruise  missiles  at  the  target 

End 

F  THE  TARGET  IS  A  SSM  SITE 
Beg 

F  THE  THREAT  IS  HIGH 

Assign  Hs  attack  aircraft 

ELSE 

Launch  Ms  cruise  missiles  at  the  target 

End 

F  THE  TAR(FT  IS  AN  AIRPLANE 
Beg 

Assign  Hf  fighters  to  engage  the  aircraft 

End 

7.  CONFLICT  ADJUDICATION 

IS  IT  TIME  TO  ADJUDICATE  THE  CONFLICT? 

NO>  exit  algorithm 

YES>  Complete  the  following  cycle  for  total  number  of  assets  fired  at  target 
Beg.  DID  THE  WEAPON  EVADE  THE  TARGET’S  DEFENSE 
SYSTEMS? 

NO>  1.  Increment  loss  counter 

2.  Evaluate  next  asset  or  exit  cycle 
YES>  Did  the  weapon  score  a  hit? 

YES>  Increment  hit  counter 

End 

Add  the  number  of  hits  during  this  engagement  to  the  targets 
number  of  hits  counter 

IS  THE  WEAPON  A  PLANE? 

NO>  Go  to  next  step 

YES>  Return  the  number  of  aircraft  that  siuvived  the  attack  to  the 
number  of  assets  available  from  the  platform  it  was 
launched 

IS  THE  TARGET  AN  AIR  BASE  OR  AIRCRAFT  CARRIER? 
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NO>Go  to  next  step 

YES>Compute  the  number  of  aircraft  destroyed  on  the 
carrier  or  air  base  by  the  attack  and  remove  them  from  the  carrier 
or  air  base  available  assets 

HAS  THE  NUMBER  OF  HITS  TO  KILL  THE  TARGET  BEEN 
REACHED? 

yes>  1 .  Remove  the  unit  from  the  situation 

2.  Update  the  counter  to  indicate  its  destruction 
no>  Indicate  that  the  unit  is  no  longer  targeted 
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APPENDIX  D.  ALGORITHM  PARAMETERS 


1.  SENSOR  DETECTION  AND  CONTACT  LOCALIZATION 

a.  Probability  of  detection  by  a  sensor 

1.  P{Sur&ce  search  radar}  =  0 .7 

2.  P{Air  search  radar}  =  0.8 

3.  P{Electromc  Support  Measures  (ESM)}  =  0.6 

b.  Accuracy  of  sensor  (a) 

1.  Sur&ce  search  radar;  U.S.  =  4  mile,  Iranian  =  4  mile 

2.  Air  search  radar:  U.S.  =  .5  mile,  Iranian  =  .5  mile 

3.  Electronic  Support  Measures:  U.S.  =4  miles,  Iranian  =  4  miles 

2.  INTELLIGENCE  UPDATE 

If  a  unit  has  been  detected,  but  bas  not  been  detected  within  the  last  T  minutes,  the  unit  is  considered 
lost,  and  its  position  is  considered  unknown.  T=  30  minutes. 

3.  WEAPON  ALLOCATION; 

a.  Number  of  assets  to  allocate  for  a  strike 

1.  12  Iranian  attack  aircraft  against  a  Hi^  threat  carrier 

2.  6  Iranian  attack  aiiciafl:  against  a  Low  threat  carrier 

3.  4  Iranian  attack  aircraft  against  a  Hi^  threat  U.S.  warsh^ 

4.  2  Iranian  attack  aircraft  against  a  Low  threat  U.S.  warship 

5.  2  Surfac&4o-sur&ioe  missiles  against  a  U.S.  warship 

6.  All  remaining  surfece-to-surfece  missiles  against  a  carrier  at  any  threat  level 

7.  1  Fighter  aircraft  against  a  U.S.  E-2C  or  reconnaissance  airplane 

8.  12  U.S.  attack  aircr^  against  a  High  threat  air  base 

9.  4  U.S.  attack  aircraft  against  a  Low  threat  airbase 

10.  6  Tomahawk  missiles  against  a  threat  air  base 
11.6  Tomahawk  missiles  against  a  Low  threat  air  base 

12.  1 U.  S.  attack  aircraft  against  a  Medium  threat  Long  range  radar  site 
13.4  Tomahawk  missUes  against  a  Medium  threat  Long  range  radar  site 

14.  6  U.S.  attack  aircraft  against  a  Medium  threat  Long  range  radar  site 

15.  1  Tomahawk  missile  against  a  Medium  threat  Lcmg  range  radar  site 

16.  1  attack  aircraft  against  a  SSM  site  at  any  threat  level 

17.  1  Tomahawk  missile  against  a  SSM  site  at  arty  threat  level 

18.  1  U.S.  fighter  aircraft  against  a  Iranian  reconnaissance  aircraft 

b.  Minimum  weapons  inventory  (reserve  on  board) 

1.  12  U.S.  attack  aircraft  on  board  a  carrier 

2.  6  U.S.  fighter  aircraft  on  board  a  carrier 

3.  20  Iranian  attack  aircraft  at  an  air  base 

4.  8  Iranian  fighter  aircraft  at  an  air  base 

5.  6  Tomahawk  missiles  on  a  U.S.  warship 

c.  Speed  of  a  weapon 

1.  400  mph  for  a  Tomahawk  missile 

2.  600  mph  for  a  U.S.  fighter  aircraft 

3.  400  mph  for  a  U.S.  attack  aircraft 

4.  500  mph  for  an  Iranian  surface-to-surface  missile 
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5.  600  mph  for  an  Iranian  figjiter  aircraft 

6.  400  mph  for  an  Iranian  attack  aircraft 

4.  CONFLICT  ADJUDICATION 

a.  Probability  that  a  weapon  is  defeated  by  the  target’s  defense  mechanisms  (Pdefeat) 

1.  P{U.S.  attack  aircraft  against  Iranian  Long  range  radar  site}  =  .3 

2.  P{U.S.  attack  aircraft  against  Iranian  SSM  site}  = .  1 

3.  P{U.S.  attack  aircraft  against  Iranian  air  base}  =  .3 

4.  P{U.S.  fighter  aircraft  against  Iranian  reconnaissance  aircraft}  =  .05 

5.  P{Tomahawk  missile  against  Iranian  long  range  radar  site}  =  .3 

6.  P{Tomahawk  missile  against  Iranian  SSM  site}  =  .2 

7.  PfTomahawk  missile  against  Iranian  air  base}  =  .3 

8.  P{Iranian  attack  aircraft  against  U.S.  carrier}  =  .3 

9.  P{lranian  attack  aircraft  against  U.S.  wardiip}  =  .2 

10.  P{Iranian  attack  aircraft  against  U.S.  merchmit}  =  .2 

11.  P{Iranian  fighter  aircraft  against  U.S.  E-2C  or  reconnaissance  aircraft}  =  .05 

12.  P{Iranian  surfece-to-surfece  missile  against  U.S.  carrier}  = .  1 

13.  P{Iranian  surfece-to-surfece  missile  against  U.S.  warsh^}  = .  1 

14.  P{Iranian  surface-to-surfece  missile  against  U.S.  merchant}  =  .1 

b.  Probability  a  weapon  scores  a  hit  against  a  target  given  the  weapon  evaded  the  defense 
mechanisms  of  the  target  (Put) 

1.  P{U.S.  attack  aircraft  against  Iranian  Long  range  radar  she}  =  .7 

2.  P{U.S.  attack  aircraft  against  Iranian  SSM  site}  =  .9 

3.  P{U.S.  attack  aircraft  against  Iranian  air  base}  =  .4 

4.  P{U.  S.  fighter  aircraft  against  Iranian  reconnaissance  aircraft}  = .  1 

5.  P{Tomahawk  missile  against  Iranian  long  range  radar  site}  =  .7 

6.  P{Tomahawk  missile  against  Iranian  SSM  site}  =  .7 

7.  P{Tomahawk  missile  against  Iranian  air  base}  =  .7 

8.  P{Iranian  attack  aircraft  against  U.S.  carrier}  =  .8 

9.  P{Iranian  attack  aircraft  against  U.S.  warship}  =  .8 

10.  P{Iranian  attack  aircraft  against  U.S.  merchant}  =  .9 

11.  P{Iranian  fighter  aircraft  against  U.S.  E-2C  or  recxrnnaissance  aircraft}  =  .1 

12.  P{Iranian  surfece-to-sur&ce  missile  against  U.S.  carrier}  =  .9 

13.  P{Iranian  surfece-to-surface  missile  against  U.S.  warship}  =  .9 

14.  P{Iranian  surface-to-surface  missile  against  U.S.  merchant}  =  .9 

c.  Aircraft  killed  on  the  ground  at  an  air  base  or  carrier  as  a  result  of  an  enemy  strike 

1.  Ac{U.S.}  =  .3 

2.  Ac{Iran}  =  .2 
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APPENDIX  E,  C++  COMPUTER  SIMULATION  CODE 


//  Edward  R.  Martinez 
//  September  1996 
//  MSDOS  6,2 

//  Borland  C-h-h  4 , 02  for  windows 
//  THIS  IS  FILE  TerTest.cpp 

#include  <f stream. h> 

#include  <iostream.h> 

#include  <stdlib.h> 

#include  <math.h> 

# include  “ terrain. h" 

main(int  argc,  char  *argv[]) 

{  int  Hal f Hour; 

//  ALLOCATE  MEMORY  FOR  TERRAIN  MAP  AND  THEN  READ  INTO  DATA  STRUCTURE 

TerrainType  myMap; 
ifstream  infile(argv[l]); 
myMap. Terraininput (inf ile)  ; 
inf ile. close ( ) ; 
myMap . Init_Uni ts ( ) ; 

//  OPEN  FILES  TO  WRITE  OUTPUT  DATA 

ofstream  outf ile (argv[2 ] ) ; 
ofstream  display { "Motion.dat" ) ; 
ofstream  losses ( “Loss . dat" ) ; 

//  VERIFY  UNIT  INITIAL  CONDITIONS 

myMap. UnitPrint (outf ile) ; 

//  SIMULATION 

for  (int  clock  =  1;  clock  <=  2880  ;  clock++)  { 

//  IDENTIFY  HALF  HOUR  &  OUTPUT  TIME  OF  SIMULATION  ON  SCREEN 

Half Hour  =  (clock  %  30  ==  0)?1:0; 
if  (clock  ==  1)  {HalfHour  =  1;} 
cout«clock«endl ; 

//  BATTLE  TIME  STEP  CYCLE  (1  MINUTE  TIME  INTERVALS) 

myMap . Movement  (HalfHour,  clock)  ; 
myMap. Sensors (clock, display) ; 
myMap. UpdateView (clock) ; 
myMap . Target ( clock) ; 
myMap .Adjudicate ( clock, losses)  ; 
myMap. ReturnAircraft (clock)  ; 

} 

//  WRITE  RESULTS  TO  OUTPUT  FILE 

myMap.  Status  (outf  ile,  clock)  ; 
myMap. WepsStatus  (outf  ile)  ; 


} 


cout«"End  of  Program" «endl  ; 
return  1 ; 
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//  Edward  R.  Martinez 
//  September  1996 
//  THESIS 
//  MSDOS  6.22 

/ /  Borland  C-h+  4.02  for  windows 
//  THIS  IS  FILE  Terrain.}! 

// 

//  This  is  the  header  file  for  the  class  TerrainType.  This  class  models 

//  region  in  the  world.  It  utilizes  class  Grid 

//  The  source  file  Terrain,  cpp  holds  the  code  for  the  class. 

#ifndef  _ Terrain_h 

#define  _ Terrain_h 

#define  FAILOPEN  'F' 

#define  READDATA 

#include  <f stream. h> 

#include  <iostream.h> 

#include  <string.h> 

#include  <stdlib.h> 

# include  <math.h> 

# include  "RdNumGen.h" 

#include  "Grid.h” 

#include  "Unit .h“ 

# include  " Stack . h " 

class  TerrainType  { 

public  : 

TerrainType ( ) ; 

-TerrainType ( ) ; 

char  Terraininput (ifstream&) ; 

void  Init_Units ( ) ; 

void  TerrainPrint (ostream&) ; 

void  UnitPrint (ostream&) ; 

void  UnitDisplay (ostream&) ; 

void  Movement  (const  int ,  const  int )  ; 

void  Sensors  (const  int,  ostreamSc)  ; 

void  UpdateView(  const  int )  ; 

void  Target (const  int )  ; 

void  Adjudicate  (const  int,  ostream&)  ; 

void  Re turnAircraft  (const  int )  ; 

void  Status  (ostreamSc,  int )  ; 

void  WepsStatus (ostreama) ; 

private  : 

//  VARIABLES 

int  NumberOf Grids,  NumBlue,  NumRed; 

GridType*  TerrainMap[701] ; 

UnitType*  ForceMap [41] ; 

RdGen  *RandNumber; 

stackType  '^'BlueTarget ,  *RedTarget,  *RetAircraf t ; 

//  PARAMETERS  AND  COUNTERS 
doub  1  e  p_mo  ve , 

Pd_Bat_Rradar , Pd„Bat_RSSM ,  Pd_Bat_Rab , Pd_Bf _Rrecon , 
Pd_Bcm_Rradar , Pd_Bcm_RSSM, Pd_Bcm_Rab, 

Pd_Rat_Bcar ,  Pd_Rat_Bwar ,  Pd_Rat_Bmerch , 

Pd_Rf _Brecon ,  Pd_RSSM_Bcar ,  Pd_RSSM_Bwar ,  Pd_RSSM_Bmerch , 
Ph__Bat_Rradar ,  Ph__Bat__RSSM ,  Ph_Bat_Rab ,  Ph_Bf  _Rr econ , 
Ph_Bcm_Rradar,  Ph_Bcm_RSSM,  Ph_Bcm_Rab, 
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Ph_Rat_Bcar ,  Ph__Rat_Bwar ,  Ph_Rat_Binerch, 

Ph_Rf __Br econ ,  Ph_RSSM_Bcar ,  Ph__RSSM_Bwar ,  Ph_RSSM_Bmerch , 

IRMsl , USMsl , IRB, USB; 

int  RHa,RLa,RHs,RLs,RMs,RHf , 

BHa , BLa , BMa , BHr , BMr , BHs , BMs , BHf , 

BlueCM,  BlueCMLoss , BlueF,  BlueA,  IBlueF,  IBlueA, BlueFLoss,  BlueALoss, 

RedF,  RedA,  IRedF,  IRedA, RedFLoss,  RedALoss,RedSSM,  RedSSMLoss, 

SP_Bcm,  SP_Bf ,  SP_Bat,  SP_Rcm,  SP_Rf ,  SP_Rat , 

Batac , Bf ac , Ratac , Rf ac , Btom, 

ICarrier,  IWarship,  IMerchant,  IBlue_Recon,  lAirbase,  ILong_Range_Radar, 

ISSM, IRed_Recon, 

ECarrier,  EWarship,  EMerchant ,  EBlue_Recon,  EAirbase,  ELong_Range_Radar, 

ESSM, ERed_Recon; 

//  FUNCTIONS 

int  maximum  (const  int  a,  const  int  b)  {  return  a>b  ?  a:b  ;} 

int  minimum  (const  int  a,  const  int  b)  {  return  a<b  ?  a:b  ;} 

double  minimum  (const  double  a,  const  double  b)  {return  a<b  ?  a:b  ;} 
int  Getnumber  ( int  Number)  (return  ForceMap  [Number]  ->GetNumber  ( )  ; } 
float  GetXcoord(int  Number)  (return  ForceMap [Number] ->GetxCoord (); } 

float  GetYcoord(int  Number)  (return  ForceMap  [Number] ->Gety Coord  (); } 

double  DetDist (const  int,  const  int); 

double  DetDist  (const  int,  const  double,  const  double); 

void  Unitinput ( ) ; 

int  DetermineGrid( const  int); 

void  StateGrid(int  Number,  int  Grid)  (ForceMap [Number] ->SpecifyGrid (Grid)  ;  } 
void  SensorLook( const  int,  const  int,  const  int,  const  int, 

const  int ,  ostream^)  ; 

int  BlueWepSel  (const  int,  int&,  const  int.,  const  int,int&); 
int  RedWepSel  (const  int,  const  int,  const  int,  int&:)  ; 

int  BStrike  (const  int,  const  int,  const  int,  const  int,  const  int,  const  double) 
int  Fight  (const  double  ,  const  double  ,  const  int,int&); 

int  RStrike  (const  int ,  const  int ,  const  int ,  const  int ,  const  int ,  const  double  )  ; 
void  UpdateCounter (const  int); 

void  Ini tCounter  (const  int ,  const  int ,  const  int  ,  const  int )  ; 
void  ACdamage(  const  int,  const  int,  const  int )  ; 

}; 

#endif 
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//  Edward  R,  Martinez 
//  September  1996 
//  THESIS 
//  MSDOS  6.22 

/ /  Borland  C-h-h  4 .  02  for  windows 
//  THIS  IS  FILE  Terrain,  cpp 
// 

/ /  The  class  Terrain  models  the  area  of  the  World  consisting  of  the  Persian 
//  Gulf,  Gulf  of  Oman  and  Iran.  It  utilizes  class  GridType  to  hold  individual 
//  square  miles  of  area. 

finclude  <f stream. h> 

#include  <iostream-h> 
finclude  <string.h> 

#include  <stdlib.h> 

# include  <math.h> 

# include  " Terrain . h " 

# include  " RdNumGen . h " 

# include  " Grid . h " 

# include  "unit .h" 

//  CONSTRUCTOR 

TerrainType : : TerrainType ( )  { 

NumberOf  Grids  =  NumBlue  =  NuinRed  =  0; 

BlueTarget  =  new  stackType ( ) ; 

RedTarget  =  new  stackType ( ) ; 

RetAircraft  =  new  stackType () ; 

RandNumber  =  new  RdGen  ( )  ; 

//  DEFINE  CONSTANT  VARIABLES 

p_move  =  0.8;  //  Probability  a  mobile  SSM  launcher  will  move  every  1/2  hr 

//  Number  of  Assets  to  use  in  an  attack 

RHa  =  12;  //  Red  attack  aircraft  against  High  threat  CV 

RLa  =  6;  //  Red  attack  aircraft  against  Low  threat  CV 

RHs  =  4;  //  Red  attack  aircraft  against  Hight  threat  warship 

RLs  =  2 ; 

RMs  =  2;  //  Red  SSM  against  U.  S.  targets 

RHf  =  1 ; 

BHa  =6;  //  Blue  attack  aircraft  against  high  threat  air  base 

BLa  =  4;  //  low  threat  air  base 

BMa  =6;  //  Blue  cruise  missile  against  air  base 
BHr  =  1;  //  Blue  cruise  missile  against  high  threat  radar  site 
BMr  =  1;  //  against  medium  threat  radar  site 

BHs  =  1;  // 

BMs  =  1;  // 

BHf  =  1;  //  Blue  fighter  against  Red  recon  plane 

//  MINIMIM  PERMISSIABLE  LEVELS  TO  STRIKE 
Batac  =  12; 

Bfac  =  10; 

Ratac  =20; 

Rf ac  =  8 ; 

Btom  =  6 ; 

//  AIRCRAFT  DAMAGED  DURING  A  STRIKE  ON  CV  OR  AIR  BASE 

IRMsl  =  .01; 

USMsl  =  .02; 

IRB  =  .02; 

USB  =  .02; 

//  COUNTER  TO  KEEP  TRACK  OF  TOTAL  FORCE  WEAPONS  USED/LOST 
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BlueCM  =  0 
Blue?  =  0 
BlueA  =  0 
IBlueF  =  0 
IBlueA  =  0 
BlueCMLoss  =  0 
BlueFLoss  =  0 
BlueALoss  =  0 
RedF  =  0 
RedA  =  0 
IRedF  =  0 
IRedA  =  0 
RedSSMLoss  =  0 
RedFLoss  =  0 
RedALoss  =  0 
RedSSM  =  0 


//  COUNTER  TO  KEEP  TRACK  OF  UNITS  LOST 

ICarrier  =  0; 

IWarship  =  0; 

IMerchant  =  0 ; 

IBlue_Recon  =  0; 
lAirbase  =  0; 

ILong_Range_Radar  =  0 ; 

ISSM  =  0; 

IRed_Recon  =  0; 

ECarrier  =  0; 

EWarship  =  0; 

EMerchant  =  0 ; 

EBlue_Recon  =  0 ; 

EAirbase  =  0;. 

ELong__Range_Radar  =  0 ; 

ESSM  =  0; 

ERed_Recon  =  0; 


SP_Bcm  =  400; 

SP_Bf  =  600; 

SP_Bat  =  400; 

SP_Rcin  =  500; 

SP_Rf  =  600; 

SP__Rat  =  400; 

//  PROBABILITY  A  STRIKE  WEAPON  IS  DEFEATED  BY  TARGET  DEFENSIVE  MECHANISMS 


Pd_Bat_Rradar  =  . 3 ; 
Pd_Bat_RSSM  =  .  1  ; 
Pd_Bat_Rab  =  .  3  ; 
Pd_Bf _Rr Gcon  =  .05; 

Pd__Bcin_Rradar  =  .  3  ; 
Pd_Bcm_:RSSM  =  ,  2  ; 
Pd_Bcin_Rab  =  .  3  ; 

Pd_Rat_Bcar  =  . 3 ; 
Pd_Rat_Bwar  =  .  2  ; 
Pd__Rat_Bmerch  =  .  2  ; 
Pd_Rf ^Brecon  =  .05; 

Pd_RSSM_Bcar  =  .  1  ; 
Pd_RSSM_Bwar  =  .  1  ; 


Pd_RSSM_Binerch  =  .1; 

//  PROBABILTY  A  STRIKE  WEAPON  SCORES  A  DIRECT  HIT  AGAINST  TARGET  AFTER  IT 
//  AVOIDED  DEFENSE  MECHANISMS 
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Ph_Bat_Rradar  =  .  7  ; 

Ph_Bat_RSSM  =  .  9  ; 

Ph_Bat_Rab  =  .  4  ; 

Ph_Bf_Rrecon  =  .1; 

Ph_Bcin_Rradar  =  .  7  ; 

Ph_Bcin_RSSM  =  .7; 

Ph_Bcm_Rab  =  . 7 ; 

Ph_Rat_Bcar  =  .  8 ; 

Ph_Rat_Bwar  =  .  8  ; 

Ph_Rat_Binerch  =  .  9  ; 

Ph_Rf_Brecon  =  .  1  ; 

Ph_RSSM_Bcar  =  .  9  ; 

Ph_RSSM__Bwar  =  .  9  ; 

Ph_RSSM_Binerch  =  .  9 ; 

for  (int  i  =  0;  i<=  700;  i++)  { 

TerrainMap [ i ]  =  new  GridType ( ) ; 

} 

for  (int  j  =  0;  j<=:  40;  j++)  { 

ForceMap [ j ]  =  new  UnitType ( ) ; 

} 

} 

//  DESTRUCTOR 

TerrainType : : -TerrainType ( )  { 

NumberO  f Gr i ds  =  0  ; 
delete  BlueTarget ; 
delete  RedTarget; 
delete  RetAircraft; 
delete  RandNnmber; 

for  (int  i  =  0;  i<=  700;  i++)  { 
delete  TerrainMap [i ] ; 

} 

for  (int  j  =  0;  j<=  40;  j++)  { 

delete  ForceMap [ j ] ; 

} 

} 

//  INPUT  THE  TERRAIN  AND  INITIAL  UNIT  DATA  BASES 

char  TerrainType: iTerrainInput (if stream^  inputfile)  { 

//  DEFINE  VARIABLES 

char  *data,  Lineof Data [160 ] ; 

int  xl,x2,yl,y2,  t,w,U,H,T,M,B,MS,MR,MSR,BR,SRR,ARR,ESMR,XP0,YP0,XPl,YPl; 
int  LPl , XP2 , YP2 , LP2 , XP3 , YP3 , LP3 , XP4 , YP4 , LP4 , XP5 , YP5 ,  LP5 , XP6 , YP6 , LP6  ; 

//  READ  IN  NUMBER  OF  GRIDS  FOR  THE  TERRAIN  AND  NUMBER  OF  UNITS  FOR  EACH  FORCE 

inputfile. getline (LineofData, 160, ' \n' ) ; 
data  =  strtok (Lineof Data, ; 

NumberOfGrids  =  atoi(data); 
data  =  strtok(NULL, ; 

NumBlue  =  atoi(data); 
data  =  strtok  (NULL,  )  ; 

NumRed  =  atoi  (data)  ; 

//  READ  IN  X  AND  Y  COORDINATES  OF  FOUR  CORNERS  OF  TERRAIN  GRID  AND  STORE  AS 
//  A  TERRAIN  MAP 

for  (int  counter  =  1;  counter  <=  NumberOfGrids;  counter++)  { 
input file. getline (LineofData, 160,  ' \n' ) ; 
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data  =  strtok(LineofData, " ,  " )  ; 

xl  =  atoi(data); 

data  =  strtok(NULL,  " , “ ) ; 

x2  =  atoi(data); 

data  =  strtok(NULL, “ , “ ) ; 

yl  =  atoi(data); 

data  =  strtok(NULL, ; 

y2  =  atoi(data); 

data  =  strtok(NULL, “ , " ) ; 

t  =  atoi(data); 

data  =  strtok(NULL,  ; 

w  =  atoi (data) ; 

//  STORE  AN  INDIVIDUAL  GRID 

TerrainMap [counter] ->SetGridPoints (xl,x2,yl,y2, t,w) ; 

} 


//  READ  IN  THE  UNIT  INITIAL  FORCE  STRUCTURE  AND  STORE  AS  A  UNIT  MAP 

for  (counter  =  1;  counter  <=  (NiomBlue  +  NumRed)  ;  counter++)  { 
inputfile.getline(LineofData,160, *\n' ) ; 
data  =  strtok(LineofData, " , " ) ; 

U  =  atoi (data); 

data  =  strtok(NULL, ", ■) ; 

H  =  atoi (data); 

data  =  strtok(NULL, ", ") ; 

T  =  atoi (data); 

data  =  strtok(NULL, ",  ")  ; 

MS  =  atoi (data); 

data  =  strtok(NULL, ", ") ; 

MR  =  atoi (data); 

data  =  strtok(NULL, ", ") ; 

M  =  atoi (data); 

data  =  strtok(NULL, ", ") ; 

MSR  =  atoi (data); 
data  =  strtok(NULL, ", ") ; 

B  =  atoi (data) ; 

data  =  strtok(NULL, ", ") ; 

BR  =  atoi (data); 

data  =  strtok(NULL, ", ") ; 

SRR  =  atoi (data); 
data  =  strtok (NULL, " , " ) ; 

ARR  =  atoi (data); 
data  =  strtok (NULL, ",") ; 

ESMR  =  atoi (data); 
data  =  strtok(NULL, ", ") ; 

XPO  =  atoi (data); 
data  =  strtok(NULL, ",  ") ; 

YPO  =  atoi (data); 
data  =  strtok (NULL, ", “); 

XPl  =  atoi (data); 
data  =  strtok (NULL, ",") ; 
yPl  =  atoi (data); 
data  =  strtok(NULL, ", ") ; 

LPl  =  atoi (data); 
data  =  strtok (NULL, “,") ; 

XP2  =  atoi (data); 
data  =  strtok (NULL, " , “ ) ; 

YP2  =  atoi (data); 

data  =  strtok (NULL, " , ") ; 

LP2  =  atoi (data); 
data  =  strtok(NULL, ", ") ; 

XP3  =  atoi (data); 
data  =  strtok (NULL, “ , " ) ; 
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YP3  =  atoi(data); 
data  =  strtok (NULL, " , “ ) ; 

LP3  =  atoi(data); 
data  =  strtok  (NULL,  ** ,  “ )  ; 

XP4  =  atoi(data); 
data  =  strtok(NULL, ; 

YP4  =  atoi (data) ; 
data  =  strtok(NULL, ; 

LP4  =  atoi (data); 
data  =  strtok (NULL, ",") ; 

XP5  =  atoi (data) ; 
data  =  strtok{NULL, ", ") ; 

YP5  =  atoi (data) ; 
data  =  strtok (NULL, ” , " ) ; 

LP5  =  atoi (data); 
data  =  strtok (NULL, ",") ; 

XP6  =  atoi (data); 
data  =  strtok (NULL, ",") ; 

YP6  =  atoi (data); 
data  =  strtok (NULL, " , ” ) ; 

LP6  =  atoi(data); 

InitCounter (T,H,M,B) ; 

//  STORE  AN  INDIVIDUAL  UNIT  IN  THE  UNIT  MAP 

ForceMap [ counter] ~>SetUnit (U, H, T, SRR, ARR, ESMR, M, B, LPl , LP2 , LP3 , LP4 , LP5 , LP6 , 
MS , MR, MSR, BR, XPO , YPO , XPl , YPl , XP2 , YP2 , XP3 , YP3 , XP4 , YP4 , XP5 , YP5 , XP6 , YP6 ) ; 

} 

return  READDATA  ; 


//  READ  IN  INITIAL  UNIT  LOCATIONS  TO  THE  GRID  MAP 

void  TerrainType : : Init_Units { )  { 

int  Gridnuit±)er; 

for  (int  i  =  1 ; i<= (NumBlue  +  NumRed) ; i++)  { 

if  (ForceMap [i] ->GetType ( )  i=  -1)  { 

Gridnuinber  =  DetennineGrid(i)  ; 

TerrainMap [Gridnuinber]  ->AddUnit  (i ,  1 ,  ForceMap  [i ]  ->GetType  ( )  , 
ForceMap [ i ] ~>GetxCoord ( ) , ForceMap [ i ] ~>GetyCoord { ) ) ; 
ForceMap [i] ->SpecifyGrid (Gridniimber)  ; 

} 

} 


//  OUTPUT  TerrainMap  TO  AN  OUTPUT  FILE 

void  TerrainType:  :TerrainPrint  (ostreamSc  method)  { 
for (int  i  =  1;  i  <=  Number Of Grids;  i++)  { 

if  (TerrainMap [i] ->AnyUnits 0 )  { 

method« "  [  “«i«  ”  ]  “ ; 

TerrainMap [i]  “>GridPrint  (method)  ; 
method«endl  ; 

} 

} 


//  OUTPUT  UnitMap  TO  AN  OUTPUT  FILE 

void  TerrainType:  :  Unit  Print  (ostream  Scmethod)  { 

for  (int  i  =  1;  i  <=  40;  i++)  (ForceMap [i ] ->PrintUnit (method) ; } 

} 

//  PRINT  OUT  X  AND  Y  COORDINATES  OF  ALL  UNITS  TO  AN'  OUTPUT  FILE 

void  TerrainType:  : Uni tDisplay  (ostream  Scmethod)  { 
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for  (int  i  =  1;  i  <=  40;  i++)  { 

ForceMap[i]->PrintLoc  (method) ; 

} 


//  PERFORM  MOVEMENT  UPDATE  STEP  FOR  ALL  UNITS  THAT  ARE  FUNCTIONING 

void  TerrainType:  : Movement  (const  int  Flag,  const  int  TheTime)  { 
int  Update  =  0,  NewGrid,  Unit , direct, AB; 
double  cs ; 


for  (int  i  =  1;  i<(NumBlue  +  NumRed  +  1);  i++)  { 

Unit  =  ForceMap [ i ] ->GetType ( ) ; 
if  (Unit  >0)  { 

if  (Unit  <  11  1  I  Unit  ==  19) {Update  =  ForceMap [i] ->MoveShip (TheTime) ; } 
else  if  (Unit  <  20)  { 

AB  =  ForceMap [i ]  *->GetLP5  0  ; 
cs  =  ForceMap [AB] ->GetCourse ( ) ; 
cs  =  (cs  >  0)?cs:~cs; 
direct  =  (cs  >1.57)?-1:1; 

Update  =  ForceMap [i] ->MoveAir (TheTime, ForceMap [AB] “>GetxCoord( ) , 

ForceMap [AB] ->GetyCoord ( ) , direct ,-1,0) 

} 

else  if  (Unit  ==24)  { 

Update  =  ForceMap[i] ->MoveSSM(TheTime,Flag,p_move)  ; 
if  (Update  ==  1)  { 

NewGrid  =  DetermineGrid ( i ) ; 

if  (TerrainMap [NewGrid] ->GetType()  ==  1)  { 

ForceMap [ i ] ->ChangeSSM ( ) ; 

} 

} 

} 

else  if  (Unit  >  30)  { 

direct  =  ForceMap [i] ->GetLP6 () ; 

AB  =  ForceMap [i] ->GetLP5 0 ; 

Update  =  ForceMap [i] ->MoveAir (TheTime, ForceMap [AB] ->GetxCoord( ) , 

ForceMap [AB] “>GetyCoord ( ) , direct ,1,1); 

} 

else  {Update  =  0;} 


} 


if  (Update  ==  1)  {ForceMap [i ] ->UpdateMot ion (. 01667 ); } 
if  (Update  >0)  { 

NewGrid  =  DetermineGrid(i) ; 

TerrainMap [ForceMap [i]->GetGr id  0 ] ->RemoveUnit ( i ) ; 
TerrainMap [NewGrid] ->AddUnit ( i , ForceMap [i ] ->GetEMC0N ( ) , 
ForceMap [i ] ->GetType ( ) , ForceMap [i ] ->GetxCoord ( ) , 
ForceMap [ i ] ->GetyCoord ( ) ) ; 
if  (NewGrid  !=  ForceMap [i ] ->GetGrid () )  { 

ForceMap [i] ->SpecifyGr id (NewGrid) ; 

} 

} 


//  LOOK  AT  THE  WORLD  WITH  EACH  OF  THREE  TYPES  OF  SENSORS 


void  TerrainType:  : Sensors  (const  int  TheTime,  ostream  £cmethod)  { 
int  Gridnumber; 
int  p  =  1; 

method«endl«"Time  =  "«TheTime«endl ; 

//  ALLOW  EACH  BLUE  UNIT  TO  USE  HIS  THREE  SENSORS  TO  LOOK  FOR  THE  ENEMY 


for  (p  =  1;  p  <  (NumBlue  +  1);  p++)  { 

method«p« "  *'  ; 
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If  (ForceMap [p] “>GetType ( )  >0)  { 

Gridnumber  =  ForceMap [p] ->GetGrid () ; 
inethod«Gridnumber« "  “ ; 

SensorLook (Gridnumber, 0, ForceMap [p] ->GetSRR( ) , TheTime, 1, method)  ; 
me  t hod« "  Sur  f  " ; 

SensorLook  (Gridnumber,  0 ,  ForceMap  [p]  ->GetARR  ( )  ,  TheTime,  2 , method)  ; 
me t hod« " Air  " ; 

SensorLook (Gridnumber, 0, ForceMap [p] ->GetESMR ( ) , TheTime, 3 , method) ; 
method«"ESM  "«endl; 

} 

} 

//  ALLOW  EACH  RED  UNIT  TO  USE  HIS  THREE  SENSORS  TO  LOOK  FOR  THE  ENEMY 

for  (p  =  (NumBlue  +  1);  p  <  (NumBlue  +  NumRed  +.  1)  ;  p++)  { 

if  (ForceMap [p] ->GetType( )  >0)  { 

method«p«  *'  " ; 

Gridnumber  =  ForceMap[p]~>GetGrid(); 
method«Gridnuinber« "  '*  ; 

SensorLook(Gridnumber,l,  ForceMap [p]->GetSRR()  , TheTime,  1 , method)  ; 
method«  ”  Surf  "  ; 

SensorLook(Gridnumber,l, ForceMap [p]->GetARR() , TheTime, 2 , method)  ; 
method«  "Air  " ; 

SensorLook  (Gridnumber,  1,  ForceMap  [p]  ->GetESMR( )  ,  TheTime,  3 ,  method)  ; 
method« "  ESM  "  «endl  ; 

} 

} 

} 


//  SEARCH  FOR  TARGETS  WITH  ALL  A  SENSOR 

void  TerrainType:  :  SensorLook  (const  int  GridNum,  const  int  BlueRed, 

const  int  sensor_range,  const  int  Timer,  const  int  SensType,  ostream  &method)  { 
int  row  =  0; 
int  col  =  0; 


} 


if  ( sensor_range  ==  0)  { 

method«"sr  =  "«sensor__range; 

if  (SensType  ==  1)  {TerrainMap [GridNum] ->LookSurf (BlueRed, Timer, method) ;} 

else  if  (SensType  ==  2)  {TerrainMap [GridNum] ->LookAir (BlueRed, Timer , method) ; } 
else  if  (SensType  ==  3)  (TerrainMap [GridNum] ->LookESM (BlueRed,  Timer , method)  ;  } 

} 


else  { 

method« "  srm  =  "«sensor_range; 

row  =  GridNum/35; 

col  =  GridNum  -  row*35; 

if  (col  ==  0)  (row — ;  col  =  35;} 

int  col_low  =  maximum ( 1,  col- sensor_range)  ; 

int  col_high  =  minimum  (35 ,  col +sensor_range)  ; 

int  row_low  =  maximum  ( 0 ,  row- sensor_range)  ; 

int  row_high  =  minimum (19 , row+ sens or_range) ; 

int  deltacol  =  col_high  -  col_low; 

int  deltarow  =  row__high  -  row_low; 

int  start_point  =  35*row_low  +  col_low; 

for  (int  r  =  0;  r  <=  deltarow;  r++)  { 

for  (int  c  =  0;  c  <=  deltacol;  C++)  { 

if  (SensType  ==  1) 

{TerrainMap[ ( start_point+35^r  +  c) ]->LookSurf (BlueRed, Timer, method) ; } 
else  if  (SensType  ==  2) 

{TerrainMap[ ( start_point+35*r  +  c) ] ->LookAir (BlueRed, Timer , method) ; } 
else  if  (SensType  ==  3) 

{TerrainMap[ (start_point+35*r  +  c) ] ->LookESM (BlueRed, Timer , method) ; } 

} 

} 

} 
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//  ACCOUNT  FOR  LOST  CONTACT 


void  TerrainType:  :UpdateView( const  int  TheTime)  { 

for(int  counter  =  1;  counter  <=  NumberOf Grids;  counter++)  { 

TerrainMap [counter ] ”>UpdateMap (TheTime) ; 

} 

} 

//  DEVELOP  TARGET  LIST 

void  TerrainType :: Target  (const  int  TheTime)  { 

int  UN,  UTfr,  UTtg,AU,  AUat,  AUms,WepType,  Threat ,  Time,  Conf=0 ,  NumAs sets; 
double  MinDistat,  MinDistms, TheDist ,  Xunit,  Yunit,  AimptError; 

//method«endl ; 

for  (int  i  =  1;  i<=  NumberOfGrids; i++)  { 

TerrainMap [ i ] ->Ini tQueue ( ) ; 

while  (TerrainMap[i] -*>LookAtView{UN,  0, Xunit,  Yunit, Conf)  )  { 

//  me thod« “BLUE [Time  =  “«TheTime«“  ]  Grid  =  "«i«"  Unit  =  “«UN«endl 
UTtg  =  ForceMap  [UN]  -->GetType  ( )  ; 

if  (ForceMap [UN] ->GetTarget ( )  ==  0  &&  UTtg  >0)  { 

Threat  =  0;  MinDistat  =  807;  MinDistms  =  807;  AUat  =  0;  AUms  =  0; 
for  (int  j  =  1;  j<=  NuinBlue;  j++)  { 

UTfr  =  ForceMap[j]~>GetType() ; 

TheDist  =  DetDist (j ,Xunit, Yunit) ; 

//  THREAT  DETERMINATION 
if  (UTfr  >  0  ScSc  UTfr  <  10)  { 

if  (TheDist<ForceMap [UN] ->GetMR( )) {Threat  =  1;} 

) 

if  (UTtg  ==  22)  (Threat  =2;} 

if  (  Threat  I =  1  &&  UTtg  ==  24)  (Threat  =  2;} 

//  DETERMINE  CLOSEST  UNITS  THAT  HAVE  WEAPONS  TO  SHOOT 
if  (UTfr  ==  1  &&  UTtg  <  30  &&  ForceMap [j ] ->GetNumMS{ ) >Batac 
&&  ForceMap [j ] ->GetMR( )  >=  TheDist)  ( 

AUat  =  (MinDistat  >  TheDist ) ?j : AUat ; 

MinDistat  =  minimum (MinDistat, TheDist ) ; 

) 

if  (UTfr  ==  2  &&  UTtg  <  30  &&  ForceMap [j ] ->GetNumMS () >Btom 
&&  ForceMap [j ] ->GetMR( )  >=  TheDist)  ( 

AUms  =  (MinDistms  >  TheDist ) ?j : AUms; 

MinDistms  =  minimum  (MinDistms,  TheDist)  ; 

) 

if  (UTfr  ==  1  ScSc  UTtg  >  30  &&  ForceMap  [j  ] ->GetNumB  ( )  >  Bfac 
ScSc  ForceMap[  j]“->GetBR( )  >=  TheDist)  ( 

AUat  =  (MinDistat  >  TheDist ) ?j :AUat; 

MinDistat  =  minimum (MinDistat , TheDist ) ; 

) 

) 

//  SELECT  WEAPON 

WepType  =  BlueWepSel  (UN, AUat, AUms, Threat, NumAs sets)  ; 

//  IF  A  WEAPON  HAS  BEEN  SELECTED,  COMPUTE  IMPACT  TIME,  AIMPOINT 
//  ERROR  AND  PLACE  ON  THE  TARGET  LIST 
if  (WepType  >0)  ( 

if  (WepType  <=  3)  (Time  =  TheTime  + 

floor (RandNumber->Norm{10, 2 ) + ( 60*DetDist (UN, 1) /SP_Bcm) ) ; } 
if  (WepType  ==  4)  (Time  =  TheTime  + 

floor  (RandN'umber->Norm(5 , 1)  +  ( 60*DetDist  (UN,  1) /SP_Bf )  )  ;  } 
if  (WepType  >=  5)  (Time  =  TheTime  + 

floor ( 60*DetDist (UN, 1) /SP_Bat ) ; } 


AimptError  =  DetDist (UN, Xunit , Yunit) ; 
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// 

// 


AimptError  =  (AimptError  <=  10)  ?exp(-‘AimptError/10.0)  :0.2; 
BlueTarget->PUSH  (Time , UN, AUat , WepType , NiomAssets , AimptError,  1.0)  ; 
method«”  "«AUafc«"  Shot  ^«NuinAssets«"  of  type  "«WepType 

«"  will  arrive  at  ”«Time«endl; 

ForceMap [UN] ->SetTarget ( 1 ) ; 

} 

} 

} 

while  (TerrainMap[i]“>LookAtView(UN,l,Xunit, Yunit,Conf ) )  {  ■ 

//  method«"  RED  [Time  -  "«TheTime«”  ]  Grid  =  "«i«"  Unit  -  "«UN«endl ; 

UTtg  =  ForceMap [UN] ~>GetType 0 ; 

if  (ForceMap [UN] ->GetTarget ( )  ==  0  &&  UTtg  >0)  { 

Threat  =  0;  MinDistat  =  807;  MinDistms  =  807;  AUat  =  0;  AUms  =  0; 

AU  =  0; 

for  (int  j  =  (NumBlue  +  1);  j<=  (NumBlue  +  NiimRed)  ;  j++)  { 

UTf  r  =  ForceMap  [j  ]  -'>GetType  ( )  ; 

TheDist  =  DetDist ( j ,Xunit , Yunit ) ; 

//  THREAT  DETERMINATION 
If  (UTfr  >  0  ScSc  UTfr  <  30)  { 

if  (TheDist  <=  ForceMap [UN] ->GetMR( ) )  {Threat  =  1;} 

} 

//  DETERMINE  CLOSEST  UNITS  THAT  HAVE  WEAPONS  TO  SHOOT 
if  (UTfr  ==  21  ScSc  UTtg  <10  &&  ForceMap  [j  ]  “>GetN\imMS  ()  >Ratac 
&&  ForceMap[ j ] ->GetMR( )  >=  TheDist)  { 

AUat  =  (MinDistat  >  TheDist ) ?j : AUat ; 

MinDistat  =  minimum (MinDistat , TheDist ) ; 

} 

if  (UTfr  ==  24  &&  UTtg  <  10  &&  ForceMap [ j ] ->GetNumMS () >  0 
ScSc  ForceMap[  j  ] ->GetMR( )  >=  TheDist)  { 

AUms  =  (MinDistms  >  TheDist ) ?j : AUms ; 

MinDistms  =  rainimiam (MinDistms, TheDist )  ; 

} 

if  (UTfr  ==  21  ScSc  UTtg  >  10  ScSc  ForceMap  [j  ] ->GetNuinB  ()  >  Rfac 
ScSc  ForceMap[  j  ] ->GetBR( )  >=  TheDist)  { 

AUat  =  (MinDistat  >  TheDist) ?j : AUat ; 

MinDistat  =  minimum (MinDistat , TheDist ) ; 

} 

} 

//  SELECT  WEAPON 
if  (AUms  >  0) 

{AU  =  AUms;WepType  =  RedWepSel  (UN,  AUms,  Threat , NiimAs sets)  ;  } 
else  if  (AUat  >  0) 

{AU  =  AUat;WepType  =  RedWepSel  (UN, AUat , Threat , NiomAssets)  ;  } 
else  {WepType  =  0 ; } 

//  IF  A  WEAPON  HAS  BEEN  SELECTED,  COMPUTE  IMPACT  TIME,  AIMPOINT 
//  ERROR  AND  PLACE  ON  THE  TARGET  LIST 
if  (WepType  >0)  { 

if (WepType  <=  13)  {Time  =  TheTime  + 

floor (RandNumber“>Norm(5, 1) + ( 60*DetDist (UN, 1) /SP_Rcm) ) ; } 
if  (WepType  ==  14)  {Time  =  TheTime  + 

floor (60*DetDist (UN, 1) /SP_Rf) ;} 
if  (WepType  >=  15)  {Time  =  TheTime  + 

floor (RandNumber->No rm( 10,2) +  (60*DetDist (UN, 1) /SP_Rat) ) ; } 
AimptError  =  DetDist (UN,Xunit, Yunit ) ; 

AimptError  =  (AimptError  <=  10 ) ?exp ( -AimptError/10 . 0) : 0 . 2 ; 
RedTarget->PUSH  (Time, UN,  AU,  WepType,  NumAssets,  AimptError,  1.0)  ; 

//  method«"  "«AU«"  Shot  "«NimiAssets«"  a  "«WepType 

//  «"  will  arrive  at  ‘'«Time«endl; 

ForceMap [UN] ->SetTarget ( 1 ) ; 


} 


} 
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} 

} 

} 

//  SELECT  A  WEAPON  TO  ATTACK  THE  RED  FORCES 

Int  TerrainType : : BlueWepSel (  const  int  UN, int  &AUatt,  const  int  AUmsl, 

const  Int  TH,  Int  &assets)  { 
int  TheType  =  ForceMap [UN] ->GetType ( ) ; 

//  HIGH  PRIORITY  THREAT 

If  (TH  ==  1)  { 

if  (AUatt  >  0  &&  TheType  ==21)  { 

assets  =  BHa  <  ForceMap [AUatt ] ->GetNuinMS () ?BHa:BLa; 
if  (assets  !=  BHa)  { 

assets  =  BLa  <  ForceMap [AUatt] ->GetNuinMS () ?BLa: ForceMap [AUatt] ->GetNuiriMS (); ) 
BlueA  =  BlueA  +  assets; 

ForceMap [AUatt] ->UseMS (assets) ; 
return  5 ; 

} 

else  if  (AUmsl  >  0  &&  TheType  ==21)  { 

assets  =  BMa  <  ForceMap [AUmsl ] ->GetNumMS ( ) ?BMa : ForceMap [AUmsl ] ->GetNtmiMS ( ) ; 
BlueCM  =  BlueCM  +  assets; 

ForceMap [AUmsl] ->UseMS (assets) ;AUatt  =  AUmsl; 
return  1 ; 

} 

else  if  (AUatt  >  0  &&  TheType  ==22)  { 

assets  =  6  <  ForceMap [AUatt ] ->GetNumMS () ?4: ForceMap [AUatt] ->GetNumMS () ; 

BlueA  =  BlueA  +  assets; 

ForceMap [AUatt] ->UseMS (assets) ; 
return  6 ; 

) 

else  if  (AUmsl  >  0  &&  TheType  ==22)  { 

assets  =  6  <  ForceMap [AUmsl ] ->GetNumMS () ?6 : ForceMap [AUmsl] ->GetNumMS () ; 

BlueCM  =  BlueCM  +  assets; 

ForceMap [AUmsl ]->UseMS (assets) ; 

AUatt  =  AUmsl; 
return  6 ; 

} 

else  if  (AUatt  >  0  &&  TheType  ==  24)  { 

assets  =  BHs  <  ForceMap [AUatt ] ->GetNumMS () ?BHs : ForceMap [AUatt] ->GetNumMS () ; 
BlueA  =  BlueA  +  assets; 

ForceMap [AUatt ] ->UseMS ( assets ) ; 

return  7 ; 

} 

else  if  (AUmsl  >  0  &&  TheType  ==24)  { 

assets  =  BMs  <  ForceMap [AUmsl ] ->GetNumMS () ?BMs : ForceMap [AUmsl] ->GetNumMS () ; 
BlueCM  =  BlueCM  +  assets; 

ForceMap [AUmsl] ->UseMS (assets) ; 

AUatt  =  AUmsl; 
return  2  ; 

} 

else  return  0 ; 


if  (TH  ==  2)  { 

if  (AUatt  >  0  &&  TheType  ==21)  { 

assets  =  BLa  <  ForceMap  [AUatt  ] ->GetNvimMS  ()  ?BLa:  0; 
if  (assets  >0)  { 

BlueA  =  BlueA  +  assets;  ForceMap [AUatt] ->UseMS (assets) ;  return  5; 

} 

} 

if  (AUmsl  >  0  &&  TheType  ==  21)  { 

assets  =  BMa  <  ForceMap [AUmsl] ->GetNumMS () ?BMa: 0 ; 
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if  (assets  >0)  { 

BlueCM  =  BlueCM  +  assets; 

ForceMap [AUmsl ] ->UseMS (assets) ; 

AUatt  =  AUmsl; 
return  1 ; 

} 

else  {return  0;} 

} 

if  (AUatt  >  0  &&  TheType  ==  24)  { 

assets  =  BHs  <  ForceMap [AUatt ] “>GetNumMS () ?BHs : 0; 
if  (assets  >0)  { 

BlueA  =  BlueA  +  assets;  ForceMap  [AUatt] ->UseMS  (assets)  ;  return  7; 

) 

} 

if  (AUmsl  >  0  &&  TheType  ==24)  { 

assets  =  BMs  <  ForceMap  [AUmsl  ]  ->GetNuinMS  ( )  ?BMs :  0  ; 
if  (assets  >0)  { 

BlueCM  =  BlueCM  +  assets; 

ForceMap [AUmsl ] ->UseMS (assets) ; 

AUatt  =  AUmsl; 
return  2 ; 

} 

else  (return  0;} 

} 

if  (AUatt  >  0  ScSc  TheType  ==22)  { 

assets  =  BHr  <  ForceMap [AUatt ] ->GetNumMS () ?BHr: 0 ; 

if  (assets  >0)  { 

BlueA  =  BlueA  +  assets;  ForceMap [AUatt] ->UseMS (assets) ;  return  6; 

} 

} 

if  (AUmsl  >  0  TheType  ==24)  { 

assets  =  BMr  <  ForceMap [AUmsl ] ->GetNumMS ( ) ?BMr : 0 ; 
if  (assets  >0)  { 

BlueCM  =  BlueCM  +  assets; 

ForceMap [AUmsl] “>UseMS (assets) ; 

AUatt  =  AUmsl; 
return  2 ; 

} 

else  (return  0;} 

} 

return  0 ; 


if  (AUatt  >  0  ScSc  TheType  ==33)  ( 

assets  =  BHf  <  ForceMap [AUatt ] ->GetNumB ( ) ?BHf : 0 ; • 
if  (assets  >0)  { 

BlueF  =  BlueF  +  assets;  ForceMap [AUatt ] ->UseB (assets) ;  return  4; 

} 

else  (return  0;} 

} 

return  0  ; 


//  SELECT  A  WEAPON  TO  ATTACK  THE  BLUE  FORCES 

int  TerrainType:  rRedWepSel  (const  int  UN,  const  int  AU,  const  int  TH,int  bassets)  ( 
int  TheType  =  ForceMap [UN] ->GetType ( ) ; 
int  AvailAssets  =  ForceMap [AU] ->GetNumMS () ; 

if  (TheType  ==  1  &&  ForceMap [AU] “>Get Type ()  ==21)  { 

if  (TH  ==  1)  (assets  =  RHa  <  Aval lAssets?RHa : AvailAssets ; 

RedA  =  RedA  +  assets;  ForceMap [AU] ->UseMS (assets) ;  return  15;} 
else  (assets  =  RLa  <  AvailAssets?RLa : AvailAssets ; 

RedA  =  RedA  +  assets ; ForceMap [AU] ”>UseMS (assets) ;  return  15;} 

} 
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if  (TheType  ==  1  &&  ForceMap [AU] ->GetType ( )  ==  24)  { 

assets  =  AvailAssets; 

ForceMap [AU] ">UseMS( assets) ; 

RedSSM  =  RedSSM  +  assets; 
return  1 1  ; 


if  (TheType  ==  2  &&  ForceMap [AU] ->Get Type ( )  ==21)  { 

if  (TH  ==  1)  {assets  =  RHs  <  Avai 1 As sets?RHs: AvailAssets; 

RedA  =  RedA  +  assets;  ForceMap  [AU] ->UseMS  (assets)  ;  return  16;} 
else  (assets  =  RLs  <  Avai lAs set s ?RLs : AvailAssets; 

RedA  =  RedA  +  assets; ForceMap [AU] ->UseMS (assets) ;  return  16;} 


if  (TheType  ==  2  ScSc  ForceMap  [AU] ->Get Type  ( )  ==  24)  { 

assets  =  RMs  <  AvallAssets?RMs:AvailAssets; 

ForceMap  [AU] ->UseMS  (a&c;ets)  ; 

RedSSM  =  RedSSM  +  RMs; 
return  12; 


if  (TheType  ==  3  &&  ForceMap [AU] ->Get Type ( )  ==  21)  { 

assets  =  RLs  <  AvailAssets?RLs : AvailAssets ; 

RedA  =  RedA  +  assets;  ForceMap  [AU] ->UseMS  (assets)  ;  return  17; 


if  (TheType  ==  3  S:&  ForceMap  [AU] ’'>GetType  ( )  ==  24)  { 

assets  =  RMs  <  Avai 1 As se t s ?RMs: AvailAssets ; 

ForceMap [AU] ->UseMS { assets )  ; 

RedSSM  =  RedSSM  +  assets; 
return  13; 


if  (ForceMap [AU] ->GetType ( )  ==  21  &&  (TheType  ==12  II  TheType  ==  13))  { 
assets  =  RHf  <  ForceMap  [AU]  ->GetNiimB  ( )  ?RHf :  ForceMap  [AU]  ->GetNuinB  ( )  ; 
ForceMap [AU]“>UseB( assets)  ; 

RedF  =  RedF  +  assets; 
return  14; 

} 

return  0  ; 


//  ADJUDICATE  ANY  CONFLICTS 

void  TerrainType:  :  Adjudicate  (const  int  TheClock,ostream  ^method)  { 
int  UN,  WepType,  Time,  Assigned_Unit  ,NumAssets; 
double  AimptErr,  dummy,  Hits_On_Target ; 
int  Blue  =  BlueTarget->stackSize ( ) ; 
int  Red  =  RedTarget->stackSize ( ) ; 

for  (int  i  =  1;  !<=  Blue;  i++)  { 

BlueTarget“>POP (Time, UN, Assigned_Unit , WepType, NumAssets,  AimptErr, dummy) ; 
if  (Time  >  TheClock)  { 

BlueTarget “>PUSH( Time, UN, Assigned_Unit, WepType, NumAssets, AimptErr, dummy) ; 

} 

else  { 

Hits_On_Target  =  BStrike  (TheClock,  WepType,  UN,  NumAssets,  Ass igned_Unit,  AimptErr) 
ForceMap  [UN]  ■“>Attack  (Hits__On_Target )  ; 
if  (WepType  ==  1)  (ACdamage (UN, Hits_On_Target , 1) ; } 
if  (WepType  ==  5)  (ACdamage (UN, Hits_On_Target , 2 ) ; } 

} 

if  (WepType  ==  4)  {ERed_Recon++; } 
if  (ForceMap [UN] ->Dead()  ==  1)  { 

UpdateCounter ( ForceMap [UN] ->GetType ()); 

ForceMap [UN] “>SetType (-1) ; 
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method«TheClock«  ”  "«UN«endl  ; 

} 

if  (ForceMap  [UN] ->Dead( )  ==  0)  {ForceMap  [UN] -*>SetTarget  (0)  ; } 


for  (i  =  1;  !<=  Red;  i++)  { 

RedTar ge t  -  >PO P  { T ime ,  UN ,  As  s  igned_Uni t ,  WepType ,  NumAs  sets,  Aimpt Er r ,  diiinmy )  ; 
if  (Time  >  TheClock)  { 

RedTarget->PUSH  ( Time ,  UN,  Assigned_Unit ,  WepType ,  NumAsset s ,  AimptErr,  dummy )  ; 

} 

else  { 

Hits_On__Target  =  RStrike  {TheClock,  WepType,  UN,  NumAs  sets,  Assigned_Unit,^^inpt  Err)  ; 
ForceMap[UN]  -  >At  tack  (Hit  s_On__Target)  ; 
if  (WepType  ==  11)  {ACdamage  (UN,  Hits_On_Target ,  1 )  ;  } 

if  (WepType  ==  15)  {ACdamage  (UN,  Hits_On_Target ,  2 )  ;  } 

} 

if  (WepType  ==  14)  {EBlue_Recon++; } 
if  ( ForceMap [UN] ->Dead( )  ==1)  { 

UpdateCoun ter (ForceMap [UN] “>GetType { ) ) ; 

ForceMap [UN] ~>SetType ( -1 ) ; 
method«TheClock«  “  “ «UN«endl«^ndl ; 

} 

if  (ForceMap [UN] ->Dead( )  ==  0)  {ForceMap [UN] ->SetTarget«))  ; } 

} 

} 


//  DETERMINE  ASSETS  LOST  AND  HITS  RECORDED  FOR  A  BLUE  ATTACK  AGAINST  RED 

int  TerrainType:  :BStr ike  (const  int  PresTime,  const  int  WT,  const  int  UN,  const  int  NA,  cons 
t  int  AU,  const  double  Err)  { 
int  RetTime, Strikes  =  0; 
int  Lost  =  0; 

if  (WT  ==  1)  { 

if  ( ForceMap  [UN] ->GetType( )  >0)  { 

Strikes  =  Fight  ( Ph_Bcm_Rab* Err ,  Pd_Bcm_Rab ,  NA,  Lost )  ; 

BlueCMLoss  =  BlueCMLoss  +  Lost; 

} 

return  Strikes; 

} 

if  (WT  ==  2)  { 

if  (ForceMap [UN] ^>GetType { )  >  0)  { 

Strikes  =Fight  (Ph_Bcm_Rradar*Err,  Pd_Bcm_Rradar,NA,  Lost)  ; 

BlueCMLoss  =  BlueCMLoss  +  Lost; 

} 

return  Strikes; 

} 

if  (WT  ==3)  { 

if  (ForceMap  [UN] -‘>GetType()  >  0)  { 

Strikes  =Fight  (Ph_Bcm_RSSM*Err ,  Pd_Bcm_RSSM,  NA,  Lost )  ; 

BlueCMLoss  =  BlueCMLoss  +  Lost; 

.  } 

return  Strikes; 

} 

if  (WT  ==  4)  { 

if  (ForceMap [UN] ->GetType ( )  >0)  { 

Strikes  =  Fight (Ph_Bf_Rrecon*Err, Pd_Bf_Rrecon,NA, Lost) ; 

BlueFLoss  =  BlueFLoss  +  Lost; 

} 

RetTime  =  PresTime  +  floor (RandNumber->Norm{90, 30)  +  { 60*DetDist (UN, AU) /SP__Bf ) ) ; 
RetAircraf t->PUSH (RetTime, AU, UN, 2 , NA-Lost , 1 . 0 , 1 . 0 ) ; 
return  (Strikes  >  0)?1:0; 

(WT  ==  5)  { 

if  (ForceMap [UN] ->GetType()  >0)  { 


} 

if 
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Strikes  =  Fight (Ph_Bat_Rab*Err, Pd_Bat_Rab,NA, Lost) ; 

BlueALoss  =  BlueALoss  +  Lost; 

} 

RetTime  =  PresTirae  +  floor  {RandNuitiber->Nonn(90, 30)  +  {60*DetDist  (UN,  AU) /SP_Bf )  )  ; 
RetAircraf t->PUSH ( RetTime , AU, UN, 1 , NA-Lost , 1 . 0 , 1 . 0 ) ; 
return  Strikes; 

} 

if  (WT  ==  6)  { 

if  (ForceMap[UN] ->GetType ( )  >0)  { 

Strikes  =  Fight (Ph_Bat_Rradar*Err, Pd_Bat_Rradar,NA, Lost) ; 

BlueALoss  =  BlueALoss  +  Lost; 

} 

RetTime  =  PresTime  +  floor (RandNumber->Norm(90,30)  +  ( 60*DetDist (UN, AU) /SP_Bf ) ) ; 
RetAircraf t - >PUSH ( RetTime , AU, UN, 1 , NA-Lost , 1 . 0 , 1 . 0 )  ; 
return  Strikes; 

} 

if  (WT  ==  7)  { 

if  (ForceMap[UN]->GetType()  >0)  { 

Strikes  =  Fight ( Ph_Bat_RSSM*Err , Pd_Bat_RSSM, NA, Lost ) ; 

BlueALoss  =  BlueALoss  +  Lost; 

} 

RetTime  =  PresTime  +  floor (RandNumber->Norm(90, 30)  +  (60*DetDist (UN, AU) /SP_Bf ) ) ; 
RetAircraf t->PUSH (RetTime, AU, UN, 1 , NA-Lost , 1 . 0 , 1 . 0)  ; 
return  Strikes; 

} 

return  0 ; 


//  DETERMINE  ASSETS  LOST  AND  HITS  RECORDED  FOR  A  RED  ATTACK  AGAINST  BLUE 

int  TerrainType:  :RStrike( const  int  PresTime, const  int  WT,  const  int  UN,  const  int  NA,  con 
St  int  AU,  const  double  Err)  { 
int  Strikes=  0; 
int  RetTime, Lost  =  0; 
if  (WT  ==  11)  { 

if  (ForceMap[UN]->GetType()  >0)  { 

Strikes  =  Fight ( Ph_RSSM_Bcar*Err , Pd_RSSM_Bcar , NA, Lost ) ; 

RedSSMLoss  =  RedSSMLoss  +  Lost; 

} 

return  Strikes; 

} 

if  (WT  ==  12)  { 

if  (ForceMap [UN] ->GetType ( )  >0)  { 

Strikes  =Fight ( Ph_RSSM_Bwar*Err, Pd_RSSM_Bwar,  NA, Lost) ; 

RedSSMLoss  =  RedSSMLoss  +  Lost; 

} 

return  Strikes; 

} 

if  (WT  ==13)  { 

if  (ForceMap [UN] ->GetType ( )  >0)  { 

Strikes  =Fight ( Ph_RSSM_Bmerch*Err , Pd_RSSM_Bmerch, NA, Lost) ; 

RedSSMLoss  =  RedSSMLoss  +  Lost; 

}  ■ 

return  Strikes; 

} 

if  (WT  ==  14)  { 

if  (ForceMap [UN] ->GetType ( )  >0)  { 

Strikes  =  Fight (Ph_Rf_Brecon*Err, Pd_Rf_Brecon,NA, Lost) ; 

RedFLoss  =  RedFLoss  +  Lost; 

} 

RetTime  =  PresTime  +  floor(RandNumber->Norm(90,30)  +  (60*DetDist(UN,AU)/SP_Bf)); 
RetAircraf t->PUSH (RetTime, AU, UN, 2 , NA-Lost , 1 . 0 , 1 . 0 ) ; 
return  (Strikes  >  0)?1:0; 

) 

if  (WT  ==  15)  { 

if  (ForceMap [UN] ->GetType()  >0)  { 
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Strikes  =  Fight  (Ph_Rat_Bcar* Err,  Pd__Rat_B car,  NA,  Lost)  ; 

RedALoss  =  RedALoss  +  Lost; 

} 

RetTime  =  PresTime  +  floor (RandNumber->Nonn  (9 0,30)  +  (60*DetDist(UN,AU)/SP_Bf)); 
RetAircraf t->PUSH (RetTime , AU, UN, 1 , NA-Lost , 1 . 0 , 1 . 0 ) ; 
return  Strikes; 

} 

if  (WT  ==  16)  { 

if  (ForceMap [UN] ->GetType ( )  >0)  { 

Strikes  =  Fight  (Ph__Rat_Bwar*Err,  Pd__Rat_Bwar,NA,  Lost)  ; 

RedALoss  =  RedALoss  +  Lost; 

} 

RetTime  =  PresTime  +  f  loor  (RandNuinber->Norm  (90 , 30 )  +  ( 60*DetDist  (UN,  AU) /SP_Bf )  )  ; 
RetAircraf t->PUSH  (RetTime ,  AU, UN,  1 , NA-Lost ,  1 , 0 , 1 . 0 )  ; 
return  Strikes; 

} 

if  (WT  ==17)  { 

if  ( ForceMap [UN] ->GetType ( )  >  0 )  { 

Strikes  =  Fight (Ph_Rat__Bmerch*Err, Pd_Rat_Bmerch, NA, Lost) ; 

RedALoss  =  RedALoss  +  Lost; 

} 

RetTime  =  PresTime  +  floor (RandNuinber->Norm  (90 , 30 )  +  ( 60*DetDist (UN, AU) /SP_Bf ) ) ; 
RetAircraf  t->PUSH  ( RetTime ,  AU,  UN,  1 , NA-Lost ,  1 . 0 , 1 . 0 )  ; 
return  Strikes; 

} 

return  0 ; 


//  RETURN  AIRCRAFT  TO  SERVICE  THAT  HAVE  LANDED  AND  RELOADED 

void  TerrainType:  :ReturnAircraft  (const  int  TheClock)  { 
int  UN,  WepType,Time,duinmyl,NumAssets; 
double  dummy 2  ,  dummy 3  ; 
int  Size  =  RetAircraf t->stackSi2e () ; 

for  (int  i  =  1;  i<=  Size;  i++)  { 

RetAircraf  t->POP  (Time,  UN,  dummyl,WepType,NumAs sets,  dummy 2 ,  dummy3  )  ; 
if  (Time  >  TheClock)  { 

RetAircraf  t->PUSH  (Time,  UN,  dummy  1  ,WepType,  NumAs sets ,  dummy2  ,  dummy3  )  ; 

} 

else  { 

if  (WepType  ==  1)  {ForceMap  [UN] ->UseMS  (-NumAs sets)  ;  } 
if  (WepType  ==  2)  (ForceMap  [UN] ->UseB (-NumAssets)  ; } 

}  •  . 

} 

} 


//  PRINT  OUT  THE  STATUS  OF  ALL  UNITS  AT  A  GIVEN  TIME 

void  TerrainType:  : Status  (ost ream  ^method,  int  TheTime)  { 
method«TheTime«Gndl ; 

UnitPrint (method) ; 

} 

//  OUPUT  WEAPONS  USED  DURING  SIMULATION 

void  TerrainType: iWepsStatus (os tream  ^method)  { 
method«endl«  "Blue"«endl 

«"CM:  "«BlueCM«”  ’•«BlueCMLoss«endl 

«“A  A/C:  •'«BlueA«”  ‘•«BlueALoss«endl 

«’‘F  A/C:  "«BlueF«''  ‘'«BlueFLoss«endl 
«”  Red"«endl 
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« '•  SSM :  “«RedSSM«  “  "«RedSSMLoss«endl 

«"A  A/C:  “«RedA«“  ''«RedALoss«endl 
«“F  A/C:  “«RedF«"  "«RedFLoss«endl ; 


//  DETERMINE  DISTANCE  BETWEEN  TWO  UNITS 

double  TerrainType:  :DetDist  (const  int  FirstUnit,  const  int  SecondUnit)  { 
double  xcomp,  ycomp; 

xcomp  =  ForceMap [FirstUnit] ->GetxCoord{ )  -  ForceMap [SecondUnit] ->GetxCoord( ) ; 
ycomp  =  ForceMap [FirstUnit] ->GetyCoord { )  -  ForceMap [SecondUnit ] ->GetyCoord () ; 
return  sqrt {xcomp*xcomp  +  ycomp*ycomp) ; 


double  TerrainType:  :DetDist  (const  Int  FirstUnit ,  const  double  X,  const  double  Y)  { 
double  xcomp,  ycomp; 

xcomp  =  ForceMap [ FirstUnit ] ->GetxCoord( )  -  X; 
ycomp  =  ForceMap [ FirstUnit ] ->GetyCoord( )  -  Y; 
return  sqrt  (xcomp*xcomp  +  ycomp*ycomp)  ; 


//  DETERMINE  THE  GRID  THAT  A  UNIT  IS  LOCATED  ON 

int  TerrainType:  :DetermineGrid( const  int  Unit)  { 
double  x,y; 
int  counter  =1; 

X  =  ForceMap [Unit ] ->GetxCoord ( ) ; 
y  =  ForceMap [Unit] ->GetyCoord ( ) ; 

while  (y>TerrainMap [counter]  -->GetyUp ( )  &&  y  >TerrainMap [counter] ->GetyDown ( ) ) 
{counter  =  counter+5;} 

while  (x>TerrainMap [counter] ~>GetxLeft ( )  &:&x>TerrainMap [counter] ->GetxRight ( ) ) 
{counter++; } 
return  counter; 


//  ADJUDICATE  A  CONFLICT  BETWEEN  TWO  UNITS 

int  TerrainType:  : Fight  (const  double  Ph,  const  double  Pd,  const  int  NC,int  ^Losses) 

{ 

int  Strike  =  0; 

Losses  =0; 

for  (int  i  =  1;  i  <=  NC;  i++)  { 

if  (RandNumber~>Unif (0,1)  >  Pd)  { 

if (RandNumber“>Unif (0, 1)  <  Ph)  {Strike++; } 

} 

else  {Losses++;} 

} 

return  Strike; 


//  KEEP  TRACK  OF  TOTAL  FORCE  UNIT  LOSSES  AFTER  A  CONFLICT  ADJUDICATION 
void  TerrainType:  :UpdateCounter( const  int  UNtype)  { 
if  (UNtype  ==  1)  {ECarrier++; } 

else  if  (UNtype  ==  2)  {EWarship++ ; } 

else  if  (UNtype  ==  3)  {EMerchant++ ; ) 

else  if  (UNtype  ==  21)  {EAirbase++; ) 

else  if  (UNtype  ==  22)  {ELong_Range_Radar++ ; } 

else  if  (UNtype  ==  24)  {ESSM++;} 


//  RECORD  THE  TYPE  AND  NUMBER  OF  EACH  TYPE  OF  UNIT  PRIOR  TO  BATTLE  COMMENCEMENT 


void  TerrainType:  :  Ini tCounter( const  int  UNtype,  const  int  HTK,  const  int  Att, 

const  int  F)  { 

if  (UNtype  ==  1)  { 
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ICarriGr++; 

IBlueA  =  IBlueA  +  Att; 

IBlueF  =  IBlueF  +  F; 

} 

else  if  (UNtype  ==  2)  {IWarship++; } 

else  if  (UNtype  ==  3)  {IMerchant++ ; } 

else  if  (UNtype  ==12  i I  UNtype  ==  13)  {IBlue_Recon  =  IBlue_Recon  +  HTK;  } 
else  if  (UNtype  ==21)  { 

IAirbase++; 

IRedA  =  IRedA  +  Att; 

IRedF  =  IRedF  +  F; 

} 

else  if  (UNtype  ==  22)  {ILong_Range_Radar++ ; } 
else  if  (UNtype  ==  24)  {ISSM++;} 

else  if  (UNtype  ==  33)  {IRed_Recon  =  IRed_Recon  +  HTK;} 


void  TerrainType:  rACdamage  (const  int  Tget,  const  int  Hits,  const  int  MorB)  { 
int  Tot Ac , AC down , At t ; 
double  Percent,  Effectiveness; 

TotAc  =  ForceMap[Tget]  “>GetNuinMS  ( )  +  ForceMap  [Tget  ]  “>GetNuinB  ( )  ; 

Effectiveness  =  RandNuinber*->Unif  (0, 2 )  ; 

if  (Tget  <  20)  {  Percent  =  (MorB  ==  1) ?USMsl :USB; } 
else  {  Percent  =  (MorB  ==  1 ) ?IRMsl : IRB; } 

ACdown  =  floor  (double  (Hits)  *Percent*Ef  f  ectiveness*double  (TotAc)); 

Att  =  floor  (double  (ACdown)  ^double  (ForceMap  [Tget  ] ->GetNiainMS  ()) /double  (TotAc)) 
Att  =  (ForceMap  [Tget] ->GetNuinMS  { )  >  Att )  ?Att :  ForceMap  [Tget] '~>GetNuitiMS  ()  ; 
ForceMap  [Tget]  ~>UseMS  (Att)  ; 

Att  =  ACdown  -  Att; 

Att  =  (ForceMap  [Tget  ] ->GetNuinB  ( )  >  Att )  ?Att :  ForceMap  [Tget  ] -“>GetNuinB  ()  ; 
ForceMap [Tget ]“>UseB( Att); 
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//  Edward  R.  Martinez 
//  September  1996 
//  MSDOS  6.22 

//  Borland  C++  4 . 02  for  windows 
//  THIS  IS  FILE  Grid.h 

//  This  class  models  a  20  square  mile  grid  of  terrain 
//  The  file  Grid.cpp  holds  the  code  for  the  class. 

#ifndef  _ ^Grid_h 

#define  _ Grid_h 

#include  <iostreain.h> 

#include  <f stream. h> 

# include  " Grid . h " 

#include  "queue. h" 

# include  " stack. h" 

#include  "RdNuinGen.h" 

//DEFINE  CLASS  GridType 

class  GridType  { 
public  : 

GridType ( ) ; 

-GridType ( ) ; 

void  UpdateMap  (const  int )  ; 

void  SetGridPoints (int ,  int ,  int ,  int ,  int ,  int ) ; 

void  AddUnit  (const  int,  const  int,  const  int,  const  double,  const  double); 

void  RemoveUni  t  ( int )  ; 

void  GridPrint  (ostream&)  ; 

void  LookSur f  ( int ,  const  int ,  ostream& )  ; 

void  LookAir  ( int ,  const  int ,  ostreamSc )  ; 

void  LookESM  ( int ,  const  int ,  ostream& )  ; 

void  Sensor_Process  (const  int,  const  int,  const  int ,  const  int,  const  double) 

int  LookAtView(int&,  const  int,  doubled,  doubled,  int&)  ; 

void  SetWeather  ( int  NewValue)  {Weather  =  NewValue;} 

double  GetWeather  ( )  (return  Weather; } 

float  GetxLeft  ( )  (return  xLeft ; } 

float  GetxRight  ( )  (return  xRight;  } 

void  InitQueue()  (BlueView->InitLookPtr  ( )  ;  RedView->InitLookPtr  ( )  ;  } 

float  GetyUp  ( )  (return  yUp;  } 

int  GetTypeO  (return  Type_of_Grid; } 

float  GetyDownO  (return  yDown; } 

int  AnyUnits ( ) ; 

private  : 

int  xLeft,  xRight,  yUp,  yDown; 
int  Weather,  Type_of_Grid; 

RdGen  *RandNum; 

queueType  *Blue,  *Red,  *BlueView,  *RedView; 

//  STATIC  VARIABLES 

static  double  ASurfB,  AAirB,  AESMB,  ASurfR,  AAirR,  AESMR; 

}; 

#endif 
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//  Edward  R,  Martinez 
//  THESIS 
//  September  1996 
//  MSDOS  6.2 

//  Borland  C++  4 . 02  for  windows 
//  THIS  IS  FILE  Grid,  cpp 
// 

//  This  class  models  a  Grid  on  a  Terrain  Map 

//  This  class  is  utilized  by  class  TerfainType  in  the  file  Terrain,  cpp 

#include  <iostreain.h> 

#include  <f stream. h> 

# i nc 1 ude  " RdNumGen . h ” 

# include  " Gr i d . h “ 

# include  ” queue. h" 

//  DEFINE  STATIC  VARIABLES 

double  GridType: :ASurfB  =  0; 
double  GridType: :AAirB  =  0; 
double  GridType: :AESMB  =  0; 
double  GridType: lASurfR  =  4; 
double  GridType:  :AAirR  =  .5; 
double  GridType:  :AESMR  =  4; 

//  CONSTRUCTOR 

GridType : : GridType ( )  { 

Blue  =  new  queueType ( ) ; 

Red  =  new  queueType  ( )  ; 

BlueView  =  new  queueType ( ) ; 

RedView  =  new  queueType () ; 

RandNum  =  new  RdGen ( ) ; 

} 

//  DESTRUCTOR 

GridType : : -GridType ( )  { 

delete  Blue; 
delete  Red; 
delete  BlueView; 
delete  RedView; 
delete  RandNum; 

} 

//  UPDATE  INTELLIGENCE  PERCEPTION 

void  GridType:  rUpdateMap  (const  int  PresentTime)  { 
int  UN,  Qsize; 

BlueView->InitLookPtr ( ) ; 

if  (! BlueView- >queueEmpty () )  { 

Qsize  =  BlueView->queueSize ( ) ; 
for  (int  i  =  0;  i  <  Qsize;  i++)  { 

if  ((PresentTime  -  BlueView->GetLookTime ( ) )  >15)  { 

UN  =  BlueView- >LookUnit 0 ; 

BlueView->IncLookPtr ( ) ; 

BlueView->RemoveMember  (UN)  ; 

} 

else  {BlueView->IncLookPtr  ( )  ;  } 

} 

} 

RedView->InitLookPtr ( ) ; 

if  ( ! Re dView->queue Empty ( ) )  { 

Qsize  =  RedView- >queueSize 0 ; 
for  (int  j  =  0;  j  <  Qsize;  j++)  { 

if  ((PresentTime  -  RedView->GetLookTime ( ) )  >15)  { 
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UN  =  RedView->l40okUnit  ( )  ; 

RedView->IncLookPtr ( ) ; 

RedView->ReinoveMember  (UN)  ; 

} 

else  {RedView->IncLookPtr ( ) ; } 

) 

} 

} 

//  INITIALIZE  A  GRID 

void  GridType:  : SetGridPoints  (int  xL,  int  xR,  Int  yU,  int  yD,  int  t,  int  w)  { 
xLeft  =  xL; 

xRight  =  xR; 

yUp  =  yU; 

yDown  =  yD ; 

Type_of_Grid  =  t; 

Weather  =  w; 

} 

//  ADD  A  UNIT  TO  A  GRID 

void  GridType :  :AddUnit  (const  int  Unit^Ntomber,  const  int  Number, 

const  int  TheTime,  const  double  XNew,  const  double  YNew) 
if  (Unit_Number  <  21)  { 

Blue- >AddFront (Uni t^Number, Number ,  TheTime,  XNew,  YNew); 

} 

if  (Unit_Number  >20)  { 

Red- >AddFront (Uni t_Number, Number,  TheTime, XNew, YNew) ; 

} 

} 

//  REMOVE  A  UNIT  FROM  A  GRID 

void  GridType:  :RemoveUnit  (int  Unit^Number)  { 
if  (Unit_Number  <  21)  { 

Blue ->RemoveMeinber  (Uni t_Number)  ; 

//  BlueView->ReinoveMeinber  (Unit^Number)  ; 

} 

if  (Unit^Number  >20)  { 

Red->RemoveMeinber  (Unit_Number)  ; 

//  RedVi  ew->RemoveMember  ( Uni  t_,Number )  ; 

} 


//  PRINT  OUT  A  LISTING  OF  GRIDS 

void  GridType: iGridPrint (ostream  &method)  { 

if  (Blue->queueSize ( )  >  0  ||  Red->gueueSize ( )  >0)  { 

method«xLef  t« "  "«xRight«‘'  ”«yUp«"  "«yDown 

«H  "«Type_of_Grid«''  Num  Blue  =  "«Blue->queueSize  ( ) 
«“  Red  Saw  =  ''«RedView->queueSize  ( ) 

«”  Num  Red  =  “«Red->queueSize  ( ) 

«"  Blue  Saw  =  ’'«BlueView->queueSize  ( )  ;  } 


int  GridType: :AnyUnits ( )  { 

if  (Blue“>queueSize( )  >  0  II  Red“>queueSize ( )  >0)  {return  1;} 
return  0  ; 


//  SEARCH  USING  SURFACE  SEARCH  RADAR 

void  GridType:  :LookSurf  (int  Color,  const  int  DeTime,  ostream  Scmeth)  { 
int  UnitNumb,  Emit; 
if  (Color  ==  0),  { 
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Red-‘>InitLookPtr  ( )  ; 

if  (  !  Red->queueEmpty  ( )  )  { 

for  (int  i  =  0;  i  <  Red->queueSize ( ) ;  i++)  { 

UnitNiimb  =  Red- >LookAt Grid  (Emit )  ; 
if  (UnitNumb  <  30  RandNiim->Unif  { 0, 1)  >.3  )  { 

Sensor_Process  (0,  DeTime,  UnitNvimb,  Emit^ASurfB)  ; 
meth« "  [  ‘•«UnitNuinb«  "  ]  "  ; 

} 

Red->IncLookPtr ( ) ; 

} 

} 

} 

else  { 

Blue->InitLookPtr ( ) ; 

if  ( !Blue“>queueEmpty { ) )  { 

for  (int  i  =  0;  i  <  Blue->queueSize ( ) ;  i++)  { 

UnitNumb  =  Blue->LookAtGrid(Emit ) ; 
if  (UnitNumb  <  10  &&  RandNum->Unif (0, 1)  >  .3)  { 

Sensor_Process (1^  DeTime,  UnitNumb,  Emit,ASurfR) ; 
meth«  “  [  '‘«UnitNumb« "  ]  "  ; 

} 

Blue->IncLookPtr ( ) ; 

} 

} 

} 

} 

//  SEARCH  USING  AIR  SEARCH  RADAR 

void  GridType:  :LookAir  (int  Color,  const  int  DeTime,  ostream  &  meth)  { 
int  UnitNumb,  Emit; 

if  (Color  ==  0)  { 

Red“>InitLookPtr ( ) ; 

if  ( i Red->queueEmpty ( ) )  { 

for  (int  i  =  0;  i  <  Red->queueSize ( ) ;  i++)  { 

UnitNumb  =  Red->LookAt Grid (Emit ) ; 
if  (UnitNumb  >  30  &&  RandNum“>Unif (0 , 1)  >.2)  { 

Sensor_Process (0,  DeTime,  UnitNumb,  Emit,AAirB) ; 
me th« "  [ "  «Uni tNumb«  “  ]  "  ; 

} 

Red->IncLookPtr ( ) ; 

} 

} 

} 

else  { 

Blue->InitLookPtr ( ) ; 

if  ( 1 Blue“>gueueEmpty ( ) )  { 

for  (int  i  =  0;  i  <  Blue->queueSize ( ) ;  i++)  { 

UnitNumb  =  Blue->LookAtGrid (Emit )  ; 

if  (UnitNumb  >  10  &&  UnitNumb  i=  16  &&  UnitNumb  !:=!?&& 

UnitNumb  1=18  &&  UnitNumb  1=  19  &&  RandNum->Unif  ( 0 , 1 )  >.2){ 
Sensor_Process (1,  DeTime,  UnitNumb,  Emit,AAirR) ; 
meth«  “  [  ”  «Uni tNumb«  ”  ]  ” ; 

} 

Blue->IncLookPtr ( ) ; 

} 

} 

} 

} 

//  SEARCH  USING  ELINT  RECEIVERS 

void  GridType ::  LookESM  (int  Color,  const  int  DeTime,  ostream  Scmeth)  { 
int  UnitNumb,  Emit; 
if  (Color  ==  0)  { 
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Red->InitLookPtr ( ) ; 

if  ( ! Red->queueEmpty ( ) )  { 

for  (int  i  =  0;  i  <  Red->queueSize ( ) ;  i++)  { 

UnitNumb  =  Red->LookAtGrid(Emit)  ; 

if  (  Emit  ==  1  &&  RandNum->Uni f ( 0 , 1 )  >  .4)  { 

Sensor_Process (0,  DeTime,  UnitNumb,  Emit,AESMB) ; 
meth« "  [  “  «Uni  tNumb«  “  ]  "  ; 

} 

Red->IncLookPtr ( ) ; 

} 

} 

} 

else  { 

Blue->InitLookPtr ( ) ; 

if  ( !Blue->queueEmpty{) )  { 

for  (int  i  =  0;  i  <  Blue->queueSize ( ) ;  i++)  { 

UnitNumb  =  Blue->LookAtGrid(Emit) ; 

if  (Emit  ==  1  Sc&  UnitNumb  !=  16  &&  UnitNumb  !=  17  && 

UnitNumb  !=  18  &&  UnitNumb  !=  19  &&  RandNum->Unif (0, 1)  >  .4){ 
Sensor_Process ( 1 ,  DeTime,  UnitNumb,  Emit,AESMR) ; 
meth«  “  [ "  «UnitNumb«  “  ]  "  ; 

} 

Blue->IncLookPtr ( )  ; 

} 

} 

} 

} 

//  SENSOR  PROCESS 


void  GridType : : Sensor_Process ( const 

const 

double  X,Y; 


int  CL,  const  int  DT,  const 
int  NU,  const  double  ACC)  { 


int  UN, 


if  (CL  ==  0)  {Red->LookUpXY(X,Y) ; } 

else  {Blue->LookUpXY(X,Y) ; } 

X  =  RandNum->Norm(X, ACC) ; 
if  (X  <  0)  {X  =  0;  } 
if  (X  >  700)  {X  =  700;) 

Y  =  RandNum->Norm(Y, ACC) ; 
if  (Y  <  0)  {Y  =  0;) 
if  (Y  >  400)  {Y  =  400;} 

if  (CL  ==  0)  C 

if  (BlueView->FindElement (UN) )  { 

BlueView->SecondSight (DT,X, Y) ; 

} 

else  {BlueView->AddFront (UN,NU,DT,X, Y) ;} 

} 

else  { 

if  (RedView->FindElement (UN) )  { 
RedView->SecondSight (DT,X, Y) ; 

} 

else  {RedView->AddFront (UN, NU, DT, X, Y) ; } 

} 

} 


//  LOOK  AT  VIEW 

int  GridType:  :LookAtView( int  &Un,  const  int  Cl,  double  ScX,  double  &Y,int  &.C)  { 
if  (Cl  ==  0)  { 

if  (BlueView->queueEmpty  ( )  )  {return  0;) 
if  (BlueView->NoLookAvail  ( )  )  (return  0;) 

Un  =  BlueView->LookUnit ( ) ; 
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C  =  BlueView->LookConf ( ) ; 
BlueView~>LookUpXY{X,Y) ; 
BlueVieW“>IncLookPtr ( )  ; 
return  1  ; 

} 

if  (Cl  ==  1)  { 

if  (RedView->queueEmpty ( ) )  {return  0;} 
if  (RedView->NoLookAvail ( ) )  (return  0;} 
Un  =  RedView~>LookUnit ( ) ; 

C  =  RedView->LookConf ( ) ; 

RedView->LookUpXY (X, Y) ; 
RedView->IncLookPtr ( )  ; 
return  1  ; 

} 

return  0 ; 
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//  Edward  R.  Martinez 
//  September  1996 
//  THESIS 
//  MSDOS  6,22 

//  Borland  C++  4 , 02  for  windows 
//  THIS  IS  FILE  Unit,h 

//  This  class  models  a  naval  or  ground  combat  unit 
//  The  file  Unit.cpp  holds  the  code  for  the  class. 
//  This  class  is  used  by  terrain,  cpp 


fifndef  Unit  h 
#define  _ ^Unit_h 

#include  <iostreain.h> 
#include  <f stream. h> 
#include  <math.h> 
#include  <stdlib.h> 

# include  " RdNumGen . h " 


//DEFINE  CLASS  UnitType 

class  UnitType  { 
public  : 

UnitType ( ) ; 

void  SetUnit  (const  int ,  const  int ,  const  int ,  const  int ,  const  int ,  const  int 

const  int ,  const  int ,  const  int ,  const  int ,  const  int ,  const  int ,  const  int , 

const  int,  const  double  , const  double  , const  double  , const  double  , 

const  double  ,  const  double  ,  const  double  ,  const  double  ,  const  doxible  , 

const  doiible  ,  const  double  ,  const  double  ,  const  double  ,  const  double  , 

const  double  , const  double  , const  double  , const  double  )  ; 

int  MoveSSM  (const  int  ,  const  int,  const  double); 

int  MoveShip(  const  int )  ; 

int  MoveAir  (const  int,  const  double,  const  double,  const  int,  const  int, 

const  int ) 

void  UpdateMot ion  (const  dotible  )  ; 
void  PrintUnit  (ostreamSc)  ; 
void  PrintLoc  (ostreamfic)  ; 
void  ChangeSSM  ( )  ; 
void  NewSSMLoc  ( )  ; 

double  GetxCoordO  {return  xLocation; } 

double  GetyCoordO  (return  yLocation; } 

int  GetNumberO  (return  Number;} 

double  GetCourse  0  (return  Course;} 

int  GetLPSO  (return  LP[5];} 

int  GetLP6  ( )  (return  LP[6];} 

void  SetType  (const  int  val)  (Type  =  val;} 

int  GetType  ( )  (return  Type;} 

int  GetGridO  (return  GridLocation; } 

,  int  GetSRRO  (return  Surf__Range; } 

int  Ge tARR  ( )  ( return  Air__Range ;  } 

int  GetESMR  ( )  (return  ESM_Range ;  } 

void  SetSRR( const  int  range)  (Surf_Range  =  range;} 

void  SetARR(const  int  range)  (Air_Range  =  range;} 

double  GetMRO  (return  Missile_Range;  } 

double  GetBR ( )  { return  Bomb_Range ; } 

int  GetNumMS ( )  (return  Missiles;} 

int  GetNumB  ( )  (return  Bombs;} 

void  UseMS  (const  int  val)  (Missiles  =  Missiles  -  val;} 
void  UseB(const  int  val)  (Bombs  =  Bombs  -  val;} 
void  ReloadMS  (const  int  val)  (Missiles  =  Missiles  +  val;} 
void  Reloads  (const  int  val)  (Bombs  =  Bombs  +  val;} 
void  Sped fyGrid( int  GridLoc)  (GridLocation  =  GridLoc;  } 
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int  GetEMCON 0  {return  Emitting;} 

void  SetTarget  (const  int  val)  (Target  =  val;} 

int  GetTarget ( )  (return  Target;} 

int  DeadO  (return  Hits  <  Hits_To_Kill?0 : 1; } 

void  Attack(const  int  Niim)  (Hits  =  Hits  +  Num; } 

private  : 

//  VARIABLES 

int  UnitNumber, Number, Type, GridLocation,  WP,  Motion; 
double  Course,  Speed,  MaxSpeed, MaxRange,  xLocation,  yLocation; 
double  XP[7],  YP[7]; 
int  LP[7]; 

int  Missiles,  Bombs,  Surf_Range,  Air_Range,  ESM_Range; 
double  Missile__Range,  Bomb_Range; 

int  Emitting,  Target,  Hits,  Hits„To_Kill , StopTime; 

RdGen  RandNumb; 

//  FUNCTIONS 

double  CheckBoundsX  ( int  & ,  const  double  )  ; 
double  CheckBoundsY(int &,  const  double); 
double  CalcCourse  (double  ,  double  )  ; 
double  Distance  ( )  ; 
void  ChartCourse  (const  int )  ; 

}; 

#endif 
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//  Edward  R,  Martinez 
//  THESIS 
//  September  1996 
//  MSDOS  6.2 

/ /  Borland  C-i-h  4 .  02  for  windows 
//  THIS  IS  FILE  Unit.cpp 
// 

//  This  class  models  a  unit  in  a  combat  scenario 


#include  <iostreain.h> 
#include  <f stream. h> 

# include  <math.h> 
#include  <stdlib.h> 
#include  <string.h> 
#include  "Unit .h" 

# i nc 1 ude  " RdNumGen . h “ 


//  CONSTRUCTOR 


UnitType: :UnitType ( )  { 
UnitNumber  =  0; 
Number  =  1 ; 

Type  =  -1; 
GridLocation  =  0; 

WP  =  1; 


Motion  =  0; 

Course  =  0.0; 

Speed  =  0.0; 

MaxSpeed  =  0.0; 
xLocation  =  0.0; 
yLocation  =0.0; 

XP[0]  =  0.0;XP[1]  =  0.0;XP[2]  =  0.0;XP[3]  =  0.0;XP[4]  = 
XP[6]  =  0.0; 

YP[0]  =  0.0;yP[l]  =  0.0;YP[2]  =  0.0, -YPO]  =  0.0;YP[4]  = 
YP[6]  =  0.0; 

LPtO]  =  0;LP[1]  =  0;LP[2]  =  0;LP[3]  =  0;LP[4]  =  0;LP[5] 
Missiles  =  0; 

Bombs  =  0; 

Surf_Range  =  0; 

Air_Range  =  0; 

ESM_Range  =  0; 

Missile_Range  =  0.0; 

Bomb^Range  =  0.0; 

Emitting  =  1;  Target  =  0;  Hits  =  0;  Hits_To_Kill  =  0; 
StopTime  =  1; 


} 


0.0;XP[5] 
0.0;yP[5] 
=  0;LP[6] 


=  0,0; 
=  0.0; 
=  0; 


void  UnitType:  :SetUnit  (const  int  UU, const  Int  UH, const  int  UT^const  int  USRR, 
const  int  UARR,  const  int  UESMR,  const  int  UM,  const  int  UB,  const  int  ULPl, 
const  int  ULP2 ,  const  int  ULP3,  const  int  ULP4,  const  int  ULP5,  const  int  ULP6, 
const  double  UMS,  const  double  UMR,  const  double  UMSR,  const  double  UBR, 
const  double  UXPO,  const  double  UYPO,  const  double  UXPl,  const  double  UYPl, 
const  double  UXP2,  const  double  UYP2,  const  double  UXP3,  const  double  UYP3 , 
const  doiible  UXP4,  const  double  UYP4,  const  double  UXP5,  const,  double  UYP5, 
const  double  UXP6,  const  double  UYP6)  { 

UnitNumber  =  UU; 

Hits_To_Kill  =  UH; 

Type  =  UT; 

xLocation  =  UXPO; 

yLocation  =  UYPO; 

MaxSpeed  =  UMS; 

MaxRange  =  UMR; 

Missiles  =  UM; 

Missile_Range  =  UMSR; 

Bombs  =  UB; 
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} 


Bomb_Range  =  UBR; 
Surf__Range  =  USRR; 
Air_Range  =  UARR; 
ESM_Range  =  UESMR; 

xp[i]  =  uxpi;  yp[i] 

XP[2]  =  UXP2;  YP[2] 
XP[3]  =:  UXP3;  YP[3] 
XP[4]  =  UXP4;  YP[4] 
XP[5]  =  UXP5;  YP[5] 
XP[6]  =  UXP6;  YP[6] 


= 

UYPl; 

LP[1] 

= 

ULPl 

=r 

UYP2; 

LP[2] 

= 

ULP2 

= 

UYP3; 

LP[3] 

= 

ULP3 

r= 

UYP4; 

LP[4] 

= 

ULP4 

= 

UYP5; 

LP[5] 

= 

ULP5 

= 

UYP6; 

LP[6] 

= 

ULP6 

int  Uni tType :  :MovgSSM (const  int  TheTime, const  int  Flag,  const  double  pmove) 


//  Determine  if  a  SSM  has  reached  the  next  waypoint,  if  it  has  determine 
//  if  it  will  loiter  at  this  point 


{ 


if  (WP  ==  2)  { 

if  (Motion  ==  1)  { 

if  (TheTime  ==  StopTime)  { 

StopTime  =  TheTime  +  LP[2]; 
Missiles  =  16; 
return  0 ; 

} 

else  return  1 ; 

} 

else  ( 

if  (TheTime  ==  StopTime)  { 

NewSSMLoc ( )  ; 

Chart Course (TheTime) ; 

WP  =  1; 
return  1  ; 

} 

return  0  ; 

} 

} 


if  (Missiles  ==  0)  { 

WP  =  2; 

Chart Course (TheTime) ; 
return  1 ; 


if  (WP  ==  1)  .{ 

if  (Motion  ==  1)  { 

if  (TheTime  ==  StopTime)  { 

Motion  =  0; 

Course  =  0 ; 

Speed  =  0 ; 
return  0  ; 

} 

else  {return  1;} 

} 

} 

if  (Flag  ==  1  ScSc  RandNumb . Uni f  ( 0 , 1 )  <  pmove)  { 

NewSSMLoc  ( )  ; 

ChartCourse (TheTime) ; 

Emitting  =  0; 
return  1  ; 

} 

if  (Flag  ==  1) (Emitting  =  (RandNumb.Unif ( 0 , 1 )  <  0.5)?0:1;} 
return  0 ; 


//  ALL  NAVAL  VESSEL  MOVEMENT 
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int  UnitType:  :MoveShip( const  Int  TheTime)  { 

//  SHIP  IS  NOT  MOVING  1)  Done  Loitering^  pull  anchor  2)  Continue  Loitering 

if  (Motion  ==  0)  { 

if  (TheTime  ==  StopTime)  { 

ChartCourse (TheTime)  ; 
return  1 ; 

} 

else  {return  0;} 


//  SHIP  IS  MOVING  1)  Reached  destination,  loiter  or  new  course 
//  2)  Maintain  present  course  and  speed 

if  (Motion  ==  1)  { 

if  (TheTime  ==  StopTime)  { 

WP++; 

if  (LP[WP-1]  >  0)  { 

StopTime  =  TheTime  +  LP[WP-1]; 

Motion  =  0; 
return  0 ; 

} 

else  { 

ChartCourse (TheTime) ; 
return  1 ; 

} 

) 

else  (return  1;} 

} 

return  1 ; 

} 

//  RECONAISSANCE  AIRCRAFT  MOVEMENT 

int  UnitType:  :MoveAir( const  int  TheTime,  const  double  X,  const  double  Y, 

const  int  dirX,  const  int  dirY,  const  int  Color)  { 


if  (WP  ==  1)  { 

if  (Motion  ==  0)  (ChartCourse (TheTime) /Emitting  =  1;} 
else  { 

if  (TheTime  ==  StopTime)  ( 

WP  =  2; 

ChartCourse (TheTime) ; 

} 

}  • 

Motion  =  1; 
return  1  ; 


if  (WP  ==  2)  ( 

if  (TheTime  ==  StopTime)  ( 
WP  =  3; 

XP[3]  =  X; 

YP[3]  =  Y; 

ChartCourse (TheTime) ; 

} 

Motion  =  1; 
return  1 ; 


if  (WP  ==  3)  { 

if  (TheTime  ==  StopTime)  { 
WP  =  4; 

xLocation  =  X; 
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yLocation  =  Y; 

Motion  =  0; 

Emitting  =  0; 

StopTime  =  TheTime  +  LP[4]; 
return  0  ; 

} 

XP[3]  =  X; 
yP[3]  =  Y; 

Chart Course (TheTime)  ; 

Motion  =  1; 
return  1  ; 

} 

if  (WP  ==  4)  { 

if  (TheTime  ==  StopTime)  { 

WP  =  1; 

YP[1]  =  Y  +(dirY*YP[5] ) ; 

XPtl]  =  X  +  (dirX*XP[5]) ; 

YP[2]  =  YP[1] ; 

if  (Color  ==  0)  {XP[2]  =  XP[1]  +  (dirX*XP[6]);} 

else  {XP[2]  =XP[1]  -  (dirX*XP [6] ) ; } 

Motion  =  0; 

} 

xLocation  =  X; 
yLocation  =  Y; 
return  2  ; 

} 

return  2 ; 


//  UPDATE  LOCATION  OF  A  UNIT  THAT  HAS  MOVED 

void  UnitType :  :UpdateMo t ion  (const  double  deltaT)  { 
int  dummy  =  0; 

xLocation  =  CheckBoundsX(diimmy,  cos  (Course)  *Speed* deltaT  +  xLocation) 
yLocation  =  CheckBoundsY (dummy, sin (Course) *Speed*deltaT  +  yLocation) 

} 

//  PRINT  OUT  INFORMATION  DESCRIBING  A  UNIT 

void  UnitType:  : PrintUnit  (ostream  Scmethod)  { 

method«UnitNumber«"  "«Type«"  HITS/HTK  =  ("«Hits«”,” 
«Hits_To_Kill«  ” )  " 

«“#  Msl  =  ”«Missiles«"  " 

"Location  =  ( "«xLocation« " ,  "«yLocation« " )  " 

«"Grid  =:  [  "«GridLocation«"  ]  " 

«endl ; 


//  PRINT  OUT  COORDINATES  TO  DATA  FILE  FOR  DISPLAY  PURPOSES 

void  UnitType:  :  PrintLoc  (ostream  Scmethod)  { 

method«UnitNumber« "  “«xLocation«"  "«yLocation«endl ; 

} 

//  ENSURE  SSM  STAYS  ON  LAND 

void  UnitType : : ChangeSSM ( )  { 

int  flag  =  0; 

XP[1]  =  CheckBoundsX(flag,  10  +  xLocation); 
if  (flag  ==  1)  { 

if  (XP[1]  ==  0)  {XP[1]  =  XP[1]  +  5;} 

else  {XP[1]  =  XP[1]  -  5; } 

} 

flag  =  0; 

YP[1]  =  CheckBoundsY (flag,  10  -  yLocation); 
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if  (flag  ==  1)  { 

if  (YP[1]  ==  0)  {YP[1]  =  YP[1]  +  5;} 

else  {YP[1]  =  YP[1]  -  5;} 

} 

) 

//  ENSURE  A  UNIT  DOES  NOT  EXCEED  X  LIMITS  OF  GRID  MAP 

double  UnitType: :CheckBoundsX(lnt  &flag,  const  double  Location)  { 
double  temp  =  Location; 

if  (Location  <0)  {  temp  =  0;  flag  =  1;} 
if  (Location  >  700)  {  temp  =  700;  flag  =1;} 
return  temp; 


//  ENSURE  A  UNIT  DOES  NOT  EXCEED  Y  LIMITS  OF  GRID  MAP 

double  UnitType: :CheckBoundsY(int  &flag,  const  double  Location)  { 
double  temp  =  Location; 

if  (Location  <0)  {  temp  =  0;  flag  =  1;} 
if  (Location  >  400)  {  temp  =  400;  flag  =1;} 
return  temp; 


//  DETERMINE  A  COURSE  TO  REACH  A  GIVEN  POINT 

doiible  UnitType:  :CalcCourse (double  X,  double  Y)  { 
double  Xval  =  (X  -  xLocation)  ; 
if  (Xval  !=  0)  { 

return  atan2 (Y-yLocation, Xval ) ; 

) 

else  { 

return  Y  >  yLocation?! . 570796 : -1 . 570796 ; 

) 


double  UnitType: : Distance ( )  { 

double  xcomp,  ycomp; 
xcomp  =  xLocation  -  XP[WP]; 
ycomp  =  yLocation  -  YP[WP] ; 
return  sqrt ( (xcomp*xcomp) + (ycomp*ycomp) ) ; 

) 

//  DETERMINE  COURSE  AND  SPEED  TO  REACH  NEXT  WAYPOINT  AND  THE  TIME  OF  ARRIVAL 

void  UnitType:  iChartCourse  (const  int  T)  { 

Course  =  CalcCourse (XP[WP] , YP[WP] ) ; 

Speed  =  -1; 

while  (Speed  <  0)  { 

Speed  =  (Type  ==  19 ) ?MaxSpeed:RandNumb.Norm (MaxSpeed, 1) ; 

} 

StopTime  =  T  +  floor (60*Distance () /Speed) ; 

Motion  =  1; 


// 

void  UnitType: :NewSSMLoc ( )  { 
int  dummy  =  0; 

XP[1]  =  CheckBoundsX(dumit^,RandNumb.  Norm  (0,10)  +  xLocation); 
yp[l]  =  Che  ckBoundsY(dxammy,RandNuirLb.  Norm  (0,10)  +  yLocation); 
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//  Edward  R,  Martinez 
//  Thesis 
//  September  1996 
//  MSDOS  6,22 

/ /  Borland  C-h-f  4,02  for  windows 
//  THIS  IS  FILE  Queue, h 

//  This  file  models  a  Queue,  It  uses  the  class  QueueElement  to  model  an 
//  element  of  the  queue.  It  is  utilized  by  the  class  GridType  in  the  file 
//  Grid.cpp 

#ifndef  _ queue_h 

#define  cmeue  h 

#include  <iostreain.h> 

#include  <f stream. h> 

# include  ’*  Queuel e .  h  “ 

//DEFINE  CLASS  queueType 

class  queueType  { 
public  : 

queueType 0  {nodeCount  =  0;  Head  =  Tail  =  NULL;} 

-queueType ( )  {Reset queue ( ) ; } 

void  AddFront  (const  int ,  const  int ,  const  int ,  const  double,  const  double) 
void  Display  (os treamSc)  ; 

int  queueEmptyO  (return  (nodeCount  >  0  ?  0  :  1);} 
int  queueSizeO  {return  nodeCount;} 
int  Remo veMember  (const  int )  ; 
void  Resetqueue ( ) ; 

int  LookConf ( )  {return  LookPtr“>GetConf idence ( ) ; } 
int  FindElement  (const  int )  ; 

void  SecondSight  (const  int,  const  double,  const  double); 
void  InitLookPtr ( )  {LookPtr  =  Head;} 
void  LookUpXY  (double  &  ,  double  &)  ; 
int  LookAtGrid ( int  & ) ; 

int  LookUnit ( )  {return  LookPtr->GetIdentif ier ( ) ; } 
int  GetLookTime  ( )  {return  LookPtr«>GetDetectTime  ( )  ;  } 
void  IncLookPtr()  {LookPtr  =  LookPtr~>GetNext ( ) ; } 
int  NoLookAvail ( )  {return  LookPtr  ==  NULL?1:0;} 
private  : 

QueueElement  *Head,  *Tail,  *LookPtr; 
int  nodeCount ; 

}; 


#endif 
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//  Edward  R.  Martinez 
//  Thesis 
//  September  1996 
//  MSDOS  6,22 

//  Borland  C-h-h  4.02  for  windows 
//  THIS  IS  FILE  queueType .  cpp 

#include  <iostreain.h> 

#include  <fstreain.h> 

#include  " queue. h" 

#include  "Queuele.h” 

void  queueType :  :  LookUpXY  (double  &Xp,  double  &yp)  { 
Xp  =  LookPtr->GetX { ) ; 

Yp  =  LookPtr-->GetY  ( )  ; 


void  queueType:  rSecondSight  (const  int  DTime,  const  double  Xn,  const  double  Yn)  { 
double  Xad  j  ,  Yad  j  ; 

Xadj  =  (LookPtr->GetX ( ) *LookPtr->GetConfidence ( )  +  Xn) / (LookPtr->GetConf idence ( ) +1) ; 
Yadj  =  (LookPtr->GetY() *LookPtr->GetConfidence()  +  Yn) / (LookPtr->GetConfidence ( ) +1)  ; 
LookPtr-'>BoostConf  idence  (DTime)  ; 

LookPtr->SetXY (Xadj , Yadj ) ; 


int  queueType:  : FindElement  (const  int  Number)  { 

InitLookPtr ( ) ; 

for  (int  i  =  0;  i  <  queueSizeO;  i++)  { 

if  (Number  ==  LookPtr->GetIdentif  ier  ( )  )  {return  1;} 
IncLookPtr ( )  ; 

} 

return  0 ; 


int  queueType :  :  LookAtGr id  ( int  &Niainber )  { 
int  Temp ; 

Temp  =  LookPtr->GetElement (Number) ; 
return  Temp; 

) 


//  Function  Resetqueue 

void  queueType :  :  Resetqueue  ( )  { 

if  (  IqueueEmpty  ( )  )  { 

QueueElement  *NEXTptr,  *CurrentPtr  =  Head; 
while  ( Current Ptr)  { 

NEXTptr  =  CurrentPtr->GetNext ( ) ; 
delete  Current Ptr; 

Current Ptr  =  NEXTptr; 

} 

nodeCount  =  0 ; 

Head  =  Tail  =  NULL; 

} 

} 

/ /  Functi  on  AddFron  t 

void  queueType:  :AddFront  (const  int  Newnode,  const  int  NewValue,  const  int  Timer, 

const  doaible  Xn,  const  double  Yn)  { 
if  ( queueEmpty  ( )  ) 

Head  =  Tail  =  new  QueueElement (Newnode, NewValue , Timer , Xn, Yn) ; 
else  (QueueElement  *temp; 

Head~>SetPrevious (new  QueueElement (Newnode, NewValue, Timer , Xn, Yn) ) ; 
temp  =  Head- >Get Previous ( ) ; 
temp->SetNext (Head) ; 

Head  =  temp; 
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} 

nodeCount++ ; 


//  Function  Display 

void  queueType:  :  Display  {os treamSc  method)  { 

if  (  !  queueEmpty  ()  )  { 

QueueElement  *CurrentNode  =  Head; 
do  { 

CurrentNode~>Display  (method)  ; 

CurrentNode  =  CurrentNode->GetNext ( ) ; 

}  while  (CurrentNode); 

} 

) 

//  Function  RemoveMewber 

int  queueType: :  Remo veMember( const  int  ElementToFind)  { 
if  ( queueEmpty  0  )  return  0; 

QueueElement  *CurrentNode  =  Head; 
do  { 

if  (ElementToFind  ==  CurrentNode->GetIdentif ier ( ) )  {. 

if  (nodeCount  ==1)  { 

delete  CurrentNode ; 

Head  =  Tail  =  NULL; 
nodeCount  =  0; 
return  1 ;  } 

if  (CurrentNode  ==  Head)  { 

Head  =  Head->GetNext ( ) ; 

Head->SetPrevious (NULL) ; 
delete  CurrentNode; 
nodeCount — ; 
return  1 ;  } 

if  (CurrentNode  ==  Tail)  { 

Tail  =  Tail->GetPrevious ( ) ; 

Tail->SetNext (NULL) ; 
delete  CurrentNode; 
nodeCount--; 
return  1 ; } 

QueueElement  *TempBef ore  =  CurrentNode ->Get Previous ( ) ; 

QueueElement  *TempAf ter  =  CurrentNode->GetNext ( ) ; 

TempBef ore->SetNext (TempAfter) ; 

TempAfter->Set Previous (TempBef ore) ; 

delete  CurrentNode ; 

nodeCount--; 

return  1 ; 

} 

CurrentNode  =  CurrentNode- >GetNext ( ) ; 

}  while  (CurrentNode) ; 
return  0 ; 
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//  Edward  R,  Martinez 
//  Thesis 
//  September  1996 
//  MSDOS  6,2 

//  Borland  C+-h  4,02  for  windows 
//  THIS  IS  FILE  QueueElem,cpp 

//  This  file  models  an  element  of  a  Queue,  It  is  utilized  by  the  class 
//  dequeType  in  the  file  queuetype ,  cpp 

#ifndef  _ QueueEle_h 

#define  _ QueueEle_h 

iinclude  <iostreain.h> 

♦include  <f stream. h> 


//DEFINE  CLASS  QueueElement 


class  QueueElement  { 
public  : 

QueueElement  (const  int  Newintl,  const  int  Newint2 , const  Int  DetTime, 
const  double  Xnew,  const  double  Ynew) 


{ 

inti  =  Newintl;  int2 

int3  =  1;  int4 

Reall  =  Xnew;  Real 2 

Real3  =  0.0;  Real4 

NEXT  =  PREVIOUS  =  NULL; 


Newint2; 
DetTime; 
=  Ynew; 

=  0.0; 

} 


-QueueElement 0  {NEXT  =  PREVIOUS  =  NULL;} 

void  Display (ostreamS:  method)  {method«intl«”  ";} 

void  SetNext (QueueElement  *nextptr)  (NEXT  =  nextptr; } 

void  Set Previous (QueueElement  *previousptr)  (PREVIOUS  =  previousptr; } 

int  Get  Element  (int  ScEleValue)  (EleValue  =  int2;  return  inti;} 

QueueElement*  GetNext ( )  (return  NEXT;} 

QueueElement*  GetPrevious ( )  (return  PREVIOUS;} 

void  Boos tConfidence (const  int  DetTime)  (int3++;  int4  =  DetTime;} 

void  SetXY(const  double  Xnew,  const  double  Ynew)  (Reall  =  Xnew;Real2=Ynew; } 

int  GetDetectTime  ( )  (return  int 4; } 

int  Getidentif ier ( )  (return  inti;} 

int  GetConf idence ( )  (return  int3 ; } 

double  GetXO  (return  Reall;} 

double  GetY  ( )  (return  Real2;} 

private  : 

QueueElement  *NEXT,  *PREVIOUS; 
int  inti,  int2,  int3,  int4; 
double  Reall ,Real2,Real3,Real4; 


#endif 
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//  Edward  R,  Martinez 
//  Thesis 
//  September  1996 
//  MSDOS  6.22 

//  Borland  4.02  for  windows 
//  THIS  IS  FILE  stack,  h 

//  This  file  models  a  stack.  It  uses  the  class  stackElement  to  model  an 
//  element  of  the  stack.  It  is  utilized  by  the  class  GridType  in  the  file 
//  Grid. cpp 


#ifndef  _ stack_h 

# define  stack_h 

#include  <iostreain.h> 

#include  <fstreain,h> 

#include  " stackele . h” 

//DEFINE  CLASS  stackType 

class  StackType  { 
public  : 

StackType  0  {nodeCount  s=  0;  Head  =  NULL;} 

-StackType ( )  {Reset stack ( ) ; } 

void  PUSH  (const  int ,  const  int ,  const  int ,  const  Int ,  const  Int , 
const  double  ,  const  double  )  ; 
void  Display  (ostreain&)  ; 

Int  stackEmpty  ( )  (return  (nodeCount  >  0  ?  0  :  1);} 
int  stackSize()  (return  nodeCount;} 

int  POP(int&,  int&,  int&,  intS:,  int double  &, double  S^)  ; 
void  Resetstack ( ) ; 

private  : 

StackElement  *Head,  *Tail; 
int  nodeCount ; 


#endif 
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//  Edward  R,  Martinez 
//  Thesis 
//  September  1996 
//  MSDOS  6.22 

/ /  Borland  C-i-+  4.02  for  windows 
//  THIS  IS  FILE  stackType.cpp 

//  This  class  models  a  stack.  It  uses  the  class  stackElement .  It  is 
//  utilized  by  the  class  Gridiype  in  file  Grid.cpp 

#include  <iostreain.h> 

#include  <f stream. h> 

#include  " stack. h" 

#include  "stackele.h" 

//  Function  Resetstack 

void  stackType : : Resetstack ( )  { 

if { ! stackEmpty ( ) )  { 

StackElement  *NEXTptr,  *CurrentPtr  =  Head; 
while  (CurrentPtr)  { 

NEXTptr  =  CurrentPtr->GetNext ( ) ; 
delete  CurrentPtr ; 

CurrentPtr  =  NEXTptr; 

} 

nodeCount  =  0; 

Head  =  NULL; 

) 

} 

//  Function  PUSH 

void  StackType:  :  PUSH  (const  int  Newl,  const  int  New2,  const  int  New3, 

const  int  New4,  const  int  New5, const  double  New6,  const  double  New?)  { 
if  ( StackEmpty ( ) ) 

Head  =  Tail  =  new  stackElement  (Newl,New2,New3, New4,New5,New6, New?)  ; 
else  { 

Tail->SetNext  (new  stackElement  (Newl ,  New2  ,  New3 , New4 ,  New5 ,  New6 ,  New? )  ) 
Tail  =  Tail~>GetNext ( ) ; 

} 

nodeCount ++; 

) 

//  Function  Display 

void  StackType: : Display (ostream&  method)  { 

if  ( I stackEmpty 0 )  { 

StackElement  *CurrentNode  =  Head; 
do  { 

CurrentNode->Display (method) ; 

CurrentNode  =  CurrentNode->GetNext ( ) ; 

}  while  (CurrentNode) ; 

} 

} 

//  Function  POP 

int  StackType:  :  POP  (int  &Eleml,  int  &Elem2,  int  ScElem3,int  &:Elem4,  int  &Elem5 
double  &Elem6,  double  &Elem?)  { 
if  ( StackEmpty  ( )  )  return  0 ; 

Eleml  =  Head“>GetElementl ( ) ; 

Elem2  =  Head“>GetElement2 ( ) ; 

Elem3  =  Head->GetElement3 ( )  ; 

Elem4  =  Head->GetElement4 ( ) ; 
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EleitiB  =  Head->GetElement5  ( ) 
Elem6  =  Head->GetEleinent6  ( ) 
Elem?  =  Head->GetElGment7 ( ) 


if  (nodeCount  ==  1)  { 

delete  Head; 

Head  =  Tail  =  NULL; 
nodeCount  =  0 ; 
return  1 ;  } 

s tackEl ement  *  temp ; 
temp  =  Head; 

Head  =  Head->GetNext ( ) ; 
delete  temp; 
nodeCount — ; 
return  1  ; 


) 


//  Edward  R.  Martinez 
//  Thesis 
//  September  1996 
//  MSDOS  6,2 

//  Borland  C’h-h  4 . 02  for  windows 
//  THIS  IS  FILE  StackEle.cpp 

//  This  file  models  an  element  of  a  Stack,  It  is  utilized  by  the  class 
//  stackType  in  the  file  stack,  cpp 

#ifndef  _ stackEle_h  ' 

#define  _ stackEle_h 

#include  <iostreain.h> 

#include  <fstreain.h> 

//DEFINE  CLASS  stackElement 

class  StackElement  { 
public  : 

StackElement  (const  Int  Elel,  const  int  Ele2,  const  Int  Ele3, 

const  int  Ele4,  const  int  Ele5,  const  double  Ele6,  const  double  Ele7)  { 
Elementl  =  Elel ; Element2  =  Ele2 ; Elements  =  Ele3 ; Element4  =  Ele4; 
Elements  =  Ele5 ; Elements  =  Ele6; Element?  =  Ele7;NEXT  =  NULL;} 

•^StackElement  ( )  {NEXT  =  NULL;} 
void  Display (os tream&  method)  { 

method«Elementl«“  “«Element2«“  “«Element3«”  ”«Element4 
«"  ••«Element5«”  "«Element6«“  “«Element7«endl; } 
void  SetNext (stackElement  *nextptr)  (NEXT  =  nextptr; } 
int  GetElementl  ( )  (return  Elementl ; } 
int  GetElement2  ( )  (return  Elements  ;  } 
int  GetElementS  ( )  (return  Elements,-} 
int  GetElement4  ( )  (return  Element4 ; } 
int  GetElementS  ( )  (return  Elements;} 
double  GetElementS  ( )  (return  Elements;  } 
double  GetElement?  ( )  (return  Element?;} 

StackElement*  GetNext ( )  (return  NEXT;} 

private  : 

StackElement  *NEXT; 

int  Elementl,  Elements,  Elements,  Element 4,  Elements,- 
double  Elements,  Element?  ; 

} ;  . 


#endif 
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//  Edward  R,  Martinez 
//  September  1996 
//  THESIS 
//  MSDOS  6.22 

//  Borland  4 . 02  for  windows 
//  THIS  IS  FILE  RdNumGen.h 

//  This  class  was  written  by  Amie  Buss^  Naval  Postgraduate  School and  this 
//  code  is  a  specific  section  of  his  Random  Number  Generation  Class 


#ifndef  _ RdNumGen_h 

# define  _ ^RdNuinGen_h 

# include  <inath.h> 

#include  <stdlib.h> 

#include  <time.h> 

const  long  MODLUS  =  2147483  6471; 
const  long  MULTI  =  241121; 
const  long  MULT2  =  261431; 

class  RdGen  { 
public  : 

RdGen ( ) ; 

RdGen  (const  long  )  ; 
double  Uniform  ( )  ; 

double  Unif  (const  double,  const  double); 
long  Unif  I  (const  long  a,  const  long  b)  ; 
double  Normal  ( )  ; 

double  Norm (const  double,  const  double); 
private  : 

long  StartingSeed, CurrentSeed; 
int  BoxFlop; 
double  BoxSave; 


}; 


#endif 
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//  Edward  R,  Martinez 
//  September  1996 
//  THESIS 
//  MSDOS  6,22 

/ /  Borland  C-h-h  4 ,  02  for  windows 
//  THIS  IS  FILE  RdNumGen.h 

//  This  class  was  written  by  Amie  Buss,  Naval  Postgraduate  School,  and  this 
//  code  is  a  specific  section  of  his  Random  Number  Generation  Class, 

//  The  one  major  modification  to  his  code  is  to  allow  the  random  number 
//  generator  seed  to  be  the  present  computer  clock  time 

# include  <math.h> 

#include  <stdlib.h> 

#include  <time.h> 

# include  " RdNumGen . h " 


RdGen : : RdGen ( )  { 

CurrentSeed  =  time (NULL); 
StartingSeed  =  CurrentSeed; 
BoxFlop  =  1; 


RdGen:  : RdGen  (const  long  initSeed)  { 
StartingSeed  =  InitSeed; 
CurrentSeed  =  initSeed; 

BoxFlop  =  1; 


double  RdGen : : Uni form ( )  { 

long  zi,  lowprd,  hi31; 


zi  =  CurrentSeed; 

lowprd  =  (zi  Sc  655351)  *  MULTI; 

hi31  =  (zi  »  16)  *  MULTI  +  (lowprd  »  16); 

zi  =  ((lowprd  Sc  655351)  -  MODLUS)  +  (  (hi31  Sc  32767)  «  16)  +  (hi31  »  15) 


if  (zi  <  0)  {  zi  +=  MODLUS;} 

lowprd  =  (zi  Sc  655351)  *  MULT 2  ; 

hi31  =  (zi  »  16)  *  MULT2  +  (lowprd  »  16) ; 

zi  =  ((lowprd  Sc  655351)  -  MODLUS)  +  (  (hi31  Sc  32767)  «  16)  +  (hi31  »  15) 


if  (zi  <0)  {  zi  +=  MODLUS;} 


CurrentSeed  =  zi; 

return  double  (((zi  »  7)  I  1)  +  1 ) /16777216 . 0  ; 


double  RdGen::  Unif  (const  double  a,  const  double  b)  { 
if  (  a  <=  b)  {return  a  +  (b  ~  a)  *  UniformO;} 
else  {  return  (a  +  b)/2-0;} 


long  RdGen:  :  Unif  I  (const  long  a,  const  long  b)  { 

if  (  a  <=  b)  {return  a  +  (long)  floor  (  (b  -  a  +  1)  *  UniformO);} 
else  (return  0;} 


double  RdGen :: Normal ( )  { 

double  vl,  v2,  w,  y; 

if  (BoxFlop)  { 
w  =  2.0; 

while  (w  >  1.0  I  I  w  <  .0001)  { 
vl  =  2.0  *  UniformO  -  1.0; 
v2  =  2.0  *  UniformO  -  1.0; 
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return  BoxSave ; 

} 

) 


double  RdGen: :  Norm  (const  double  mean,  const  double  std)  { 
If  (std  <=  0.0)  {return  mean;} 
else  (return  mean  +  std  *  Normal  (); } 

) 
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APPENDKF.  SIMULATION  RUN  DATA 
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98 
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MOE  1  vs.  Trial#  [10  replications  of  each  trial] 


MOE  2  vs.  Trial#  [1 0  replications  of  each  trial] 
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multinational  forca  ships  survived  battle  Fraction  of  CVBG  ships  survived  battle 


MOE  3  vs.  Trial#  [1 0  replications  of  each  trial] 
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MOE  5  vs.  That#  [1 0  replications  of  each  trial] 
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Iranian  attack  a/c  launched  against  U.  S.  ships 


Trial 
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